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SHOCK  INDUCED  DETONATION  ON  PROJECTILES  IN  HYPERSONIC  FLOWS  OF  DETONABLE 

GASMDCTURES 

This  Final  Report  on  Contract  N6817194C9065  presents  the  devel(q)ments  in  the  research  program  in  the 
period  1  May  1994  to  ^ril  30, 1995.  The  research  effort  and  results  which  were  obtained  in  the  period 
1994  to  January  1995  were  presented  in  the  Interim  Reports  1,  2  and  3.  During  the  period  of  this  research 
contract  the  investigations  included  the  following  programs:  A  CFD  results  on  the  External  Propulsion 
Accelerator  (EPA)  projectile  configurations.  B.  An  analytical  study  on  the  stability  of  Itypersonic  reacting 
flow  at  the  stagnation  region  of  a  blunt  bo(fy  using  Dynamical  System  Analysis.  C.  The  Use  of  the  EPA  for 
Sciamjet  Combustion  ResearcL  D.  The  Use  of  the  EPA  for  Hypersonic  Aeroctynamic  Test  Facility.  E. 
Analysis  of  the  Imtiation  of  Detonation  on  a  Hypervelocity  Prcyectile  and  It’s  Maximum  Velocity  in  the  EPA. 
F.  Prq)arations  for  testing  at  the  Amty  Research  Laboratories  Ram  Accelerator  fecility  at  Aberdeen,  MD.  The 
results  of  these  investigations  are  summarized  in  this  report  and  presented  in  its  Annexes.  A  summaiy  of  the 
research  results  will  be  presented  at  the  31st  AIAA  Joint  Propulsion  Conference,  July  1995  (Ref.  1)  and  is 
included  in  Annex  1 

The  research  performed  under  this  contract  provides  a  numerical  proof  of  concqit  for  the  principles  of  the 
External  Propulsion  Accelerator.  The  numerical  results  indicate  that  considerable  thrust  can  be  obtained  on  a 
properly  designed  projectile  enabling  accelerations  of  tens  of  thousands  of  g’s.  Such  accelerator  can  be 
utilized  to  increase  the  muzzle  velocity  of  conventional  gun  projectiles  as  well  as  for  new  artillery  tystems. 
The  EPA  can  be  used  for  launching  hypervelocity  missiles  for  theater  missile  defense,  either  for  terminal 
intercept  or  for  boost  phase  intercept.  A  very  intriguing  application  can  be  for  single  stage  to  orbit  as  the 
Piaximum  velocity  achievable  in  the  EPA  is  estimated  to  be  6  times  the  detonation  velocity,  well  b^ond 
escape  velocity  fi’om  the  earth  gravitational  force.  It  is  now  necessary  to  carry  out  experimental  tests  to  verify 
the  concept.  A  single  firing  is  planed  to  be  performed  at  ARL  in  the  Ml  of  1995,  as  a  result  of  the  present 
program.  However,  a  much  more  extensive  test  program  is  required  in  order  to  have  a  reasonable  proof  of 
concept  experimental  verification.  This  remains  to  be  done  in  the  next  phase  of  the  research  program. 

1.  SCIENTMC  WORK 

a.  Introduction 

The  question  of  the  establishment  and  stabilization  of  a  combustion  or  detonation  from  on  a  projectile  flying 
at  hypersonic  speeds  in  a  detonable  mixture  is  of  interest  both  as  a  fundamental  combustion  problem  and  for 
its  practical  applications  in  the  External  Propulsion  Accelerator  and  it  is  also  relevant  to  the  Oblique 
Detonation  Wave  engine,  the  Ram  Accelerator  and  various  scramjet  engines.  The  External  Propulsion 
Accelerator  was  proposed  by  J.  Rom  (Ref.  2)  and  certain  aspect  of  its  characteristics  were  presented  in  Refs.  3 
and  4.  Some  analytical  investigations,  based  on  many  simplifying  assumptions,  for  the  ftigtahlishmpnf  of 
combustion  firont  ahead  of  the  forward  facing  step  and  on  ^herical  nosed  blunt  bodies  in  hypersonic  flows  of 
detonable  mixtures,  including  the  studies  of  the  oscillations  that  may  occur  were  investigated  by  Tivanov  and 
Rom  (Ref.  5, 6  and  7). 

The  External  Propulsion  Accelerator  has  been  developed  following  the  development  of  the  Ram  Accelerator 
by  A.  Hertzberg  (Ref.  8)  at  the  University  of  Washington.  The  Ram  Accelerator  method  utilizes  the  idea  of 
using  a  premixed  fiiel/oxidizer  mixture  in  the  launcher  barrel  and  injecting  a  projectile  at  supersonic  speeds 
into  this  mixture.  Thrust  is  generated  by  the  ramjet  combustion  obtained  between  the  tube  wall  and  the 
projectile,  where  the  projectile  acts  as  a  ramjet  centerbody  and  the  barrel  is  basically  a  long  cowling,  and 
combustion  is  initiated  and  stabUized  by  the  interaction  of  the  shock  waves  fi’om  the  projectile  nose  with  the 
barrel  walls.  However,  in  deference  fi’om  the  Ram  Accelerator,  in  the  External  Propulsion  Accelerator  the 
thrust  is  generated  purely  by  the  interactions  between  the  shock  waves  fiom  the  projectile  nose  with 
combustion  generated  by  normal  shock  waves  generated  on  the  projectile  body  by  a  ramp  or  a  forward  faring 
step  or  Ity  a  ring  wing,  completely  independent  of  the  barrel  walls.  Since  the  projectile  flies  in  a  premixed 
fiiel/oxidizer  atmosphere,  the  difficult  problem  of  mixing  the  fuel  and  oxidizer,  which  plagues  all 
airbr^thing  propulsion  methods  proposed  for  hypersonic  flight,  is  eliminated.  However,  there  are  the 
remaining  difficult  issues  of  the  ability  to  initiate  and  stabilize  the  combustion  front  on  the  projectile  in  such 
a  way  that  positive  thrust  can  be  generated.  The  previous  studies  (Refs.  3  and  4)  indicated  that,  if  a  relatively 
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small  forward  feeing  stq>  on  the  projectile  shoulder  will  be  able  to  induce  a  detonation  wave,  then  reasonable 
thrust  levels,  about  the  same  as  those  obtained  in  the  Ram  Accelerator,  can  be  achieved. 

In  order  to  further  investigate  the  flow  over  the  configuration  with  a  forward  feeing  step  a  cooperative 
program  was  established  with  the  Ram  Accelerator  project  group  at  the  Army  Research  Laboratory  in 
Aberdeen,  MD.  During  the  period  from  March  1994  to  September  1994  the  Principal  Investigator,  J.  Rom, 
was  at  the  University  of  Maryland  and  established  this  cooperation  with  the  Ram  Accelerator  Project  group  at 
the  Aberdeen  Proving  Ground  and  with  the  Hypersonic  Research  Center  of  the  department  of  Aerospace 
Engineering  at  the  University  of  Maryland,  College  Park,  MD.  During  this  period  a  very  significant  progress 
was  achieved  in  this  research  program.  The  main  results  are  presented  under  the  following  subjects: 

A.  CFD  results  on  the  External  Propulsion  Accelerator  (EPA)  projectile  configurations. 

The  application  of  the  solution  of  the  Navier-Stokes  equations  for  flows  with  chemical  processes  to  the  flow 
around  the  EPA  projectile  was  achieved  by  the  adaptation  of  the  Ram  Accelerator  calculations  done  at  ARL 
by  Nusca  (Refs.  9,  10)  to  the  External  Propulsion  projectile  geometry.  This  calculation  uses  the  Rockwell 
Science  Center  USA-RG  (Unified  Solution  Algorithm  Real  Gas)  code  written  by  Chakravarthy  et  al  (Refe.  11, 
12).  This  CFD  code  includes  the  full  3D  unsteaefy  Ronolds-Averaged  Navier-Stokes  (RANS)  equations 
including  equations  for  chemical  kinetics  (finite-rate  and  equilibrium).  These  equations  are  cast  in 
conservation  form  and  converted  to  algebraic  equations  using  upwind  finite-difference  and  finite-volume 
formulations.  The  equations  are  solved  using  a  second-order  TVD  (total  variation  diminishing)  scheme 
which  is  used  to  insure  non-oscillatory  numerical  behavior.  The  results  of  these  calculations  indicated  that 
combustion-detonation  is  established  on  forward  feeing  steps  of  height  greater  than  0.5  mm.  Based  on  these 
calculation  a  projectile  configuration  giving  reasonably  high  thrust  levels  was  selected  for  an  experimental 
shot  in  the  Ram  Accelerator  fecility  at  ARL.  The  experiment  will  be  performed  using  a  conical  nosed  model 
with  a  forward  feeing  step  which  will  be  shot  into  methane-oxygen-nitrogen  mixture  in  the  ARL  Ram 
Accelerator  fecility  (Refe.  13,  14)  using  subcaliber  diameter  model.  Part  of  the  results  of  these  calculations 
are  presented  in  AIAA  Paper  95-0259  (Ref.  15)  and  is  included  in  Annex  2. 

The  first  projectile  design  is  composed  of  a  shallow  sharp  nosed  cone  (10°  half  angle)  followed  by  a  forward 
feeing  step,  a  short  length  of  a  cylindrical  body  and  a  conical  afteibody  with  a  blunt  base.  In  the  present 
calculations  effects  of  varying  step  height  on  generating  the  combustion  front  are  investigated.  It  is  shown 
that  for  steps  of  height  above  0.5  mm  at  flow  Mach  number  of  6  there  is  a  sizable  combustion  front  followed 
by  combustion  over  the  afterbody  and  in  the  wake  behind  the  blunt  base.  The  combustion  results  in  a  region 
of  high  pressure,  high  temperature  gas  engulfing  the  back  parts  of  the  projectile  resulting  in  significant  net 
positive  thrust.  This  region  of  combustion  is  confined  by  the  hypersonic  flow  structure  generated  the 
interactions  of  the  shock  wave  from  the  conical  nose  with  the  detached  shock  wave  ahead  of  the  step  which  is 
followed  by  the  combustion  -  detonation  front.  This  is  an  excellent  demonstration  of  external  propulsion  in 
which  supersonic  combustion  is  confined  by  the  external  flow  without  any  structural  confinement. 

The  calculations  program  includes  the  parametric  study  of  the  effects  of  step  height,  flow  pressure  and  Mach 
number  and  projectile  geometry  on  the  combustion  process  and  their  effects  on  the  thrust  generated  on  the 
projectile.  Particular  emphasis  is  placed  on  the  investigation  of  the  flow  ahead  of  the  forward  facing  step, 
including  the  formation  of  the  combustion-detonation  front.  Some  of  these  results  are  presented  in  Ref.  15 
and  Annex  2. 

B.  Studies  of  the  Stability  of  Detonations  in  the  Hypersonic  Flow  using  Dynamical  System  Analysis. 

The  analytical  studies  of  the  combustion-detonation  in  the  stagnation  region  of  blunt  nosed  bodies  at 
hypersonic  speeds  indicated  a  veiy  interesting  analogy  between  the  combustion-detonation  phenomena  in  the 
hypersonic  stagnation  flow  and  the  mechanical  -  electrical  vibrating  systems.  Therefore,  the  methods  of 
Dynamical  systems  analysis  can  be  applied  to  the  combustion-detonation  process  resulting  in  very  good 
results.  Some  of  this  work  is  presented  in  a  paper  entitled:  “Stability  of  Hypersonic  Reacting  Stagnation  Flow 
of  a  Detonable  Gas  Mixture  by  Dynamical  Systems  Analysis”  by  Tivanov  and  Rom.  This  paper  is  accepted 
for  publication  in  Combustion  and  Flame  (Ref.  16)  and  the  paper  is  presented  in  Annex  3. 
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C.  The  Use  of  the  EPA  for  Scramjet  Combustion  Research. 

The  possibility  of  using  a  ring  wing  model  for  scramjet  combustion  research  is  being  studied  first  by  CFD 
solutions  on  the  ring  wing  configuration.  The  projectile  in  this  case  will  consist  of  a  nose  cone,  cylindrical 
mid-section  (no  step!)  and  aft-cone  shaped  to  conform  to  the  scramjet  engine  desigrt  A  ring,  which  will 
simulate  the  scramjet  engine  cowling  will  be  attached  so  that  the  projectile  will  act  as  the  scramjet  anginp 
centerbody.  Firing  such  projectiles  in  the  EPA  will  enable  studies  of  combustion  characteristics  in  scramjet 
engine  configurations  using  various  premixed  fuels.  In  particular  this  fecility  can  be  used  for  studies  of 
initiation  of  supersonic  combustion  and  dynamics  of  combustion  of  premixed  hydrocarbon  fuels.  This  concept 
is  presented  in  a  paper  at  the  35th  Israel  Annual  Conference  on  Aeroqrace  Sciences,  held  in  Tel-Aviv  and 
Hai&  in  February  15-16,  1995  (Ref.  16)  and  is  included  in  Armex  4. 

D.  The  Use  of  the  EPA  for  Hypersonic  Aerodynamic  Test  Facility 

It  is  proposed  to  use  the  EPA  as  the  launcher  of  aeroctynaniic  models  to  hypersonic  speeds  in  a  ground  test 
facility,  pie  models,  encased  in  a  properly  designed  sabot,  are  accelerated  in  the  EPA  to  the  desired 
hypersonic  speed  and  then  fired  into  an  instrumented  ballistic  range.  Since  the  EPA  is  capable  theoretically  to 
launch  projectiles  to  velocity  of  up  to  about  6  times  the  detonation  velocity  of  the  fuel-oxidizer  mixture  which 
can  be  in  the  order  of  12,000  m/sec  to  18,000  m/sec,  it  is  well  suited  to  be  a  launcher  for  hypersonic  ballistic 
ranges.  Using  enclosed  baUistic  ranges  will  enable  simulation  of  various  altitudes  in  air  or  in  any  other  gas 
mixture,  which  can  be  used  to  simulate  atmospheres  of  other  planets.  The  concept  is  presented  at  the  AIAA 
Hypersonic  Technology  Conference  held  in  Chattanoga,  Term,  on  April  3-7,  1995  (Ref  18)  and  is  attached  as 
Annex  5. 

E.  Analysis  of  the  Initiation  of  Detonation  on  a  Hypervelodty  Projectile  and  It’s  Maximum  Velocity  in  the 
EPA. 

The  initiation  and  the  subsequent  acceleration  of  the  projectile  in  the  in-tube  chemical  launchers  for 
accelerating  projectiles  to  hypervelodty  are  investigated  using  an  energy  balance  analysis.  It  is  shown  that  the 
energy  balance  analysis  can  be  used  to  evaluate  the  velocity  required  to  initiate  the  detonation  process  on  the 
rear  part  of  the  projectile.  This  analysis  indicated  that  the  initiation  of  detonation  by  a  step  on  the  projectile  in 
the  EPA  as  well  as  in  the  Ram  Accelerator  can  be  achieved  at  velodties  lower  then  the  detonation  velodty. 

The  energy  balance  analysis  can  be  used  to  evaluate  the  maximum  velodty  that  can  be  achieved  in  the  in-tube 
chemical  accelerators  by  balancing  the  drag  work  of  the  projectile  with  the  utilization  of  the  total  available 
chemical  reaction  energy.  It  is  shown  that  the  EPA  is  capable  to  reach  maximum  velodties  of  about  6  times 
the  detonation  velodty  while  the  Ram  Accelerator  is  limited  to  about  1.3  times  the  detonation  velodty  of  the 
mixture.  This  analysis  is  presented  in  Technion  Aerospace  Engineering  (TAE)  Report  729,  December  1994 
(Ref  19)  which  is  enclosed  as  Armex  6.  The  paper  is  accepted  for  presentation  at  the  20th  International 
Symposium  on  Shock  Waves,  July  23-28  1995  in  Pasadena,  CA 

F.  Preparations  for  Test  Firing  at  the  ARL  Fadlity  at  the  Aberdeen  Proving  Grounds,  MD.  USA. 

In  order  to  prove  the  concept  of  the  external  propulsion  on  the  step  model  a  test  program  was  initiated  at  the 
ARL  Ram  Accelerator  facility.  The  first  experiment  was  performed  using  the  120mm  gun  launcher  using  a 
conical  nosed  model  vrith  a  forward  facing  step.  The  first  test  of  a  subcaliber  model  of  32  mm  diameter  was 
tested  in  the  120  mm  Ram  Accelerator  &cility  at  ARL  on  July  26,  1994.  The  projectile  was  not  accelerated 
correctly  and  disintegrated  in  the  launching  gun  probably  due  to  the  sabot  failure. 

For  the  experimental  validation,  a  new  test  is  planned,  probably  in  the  fall  of  1995,  which  will  be  done 
utilizing  the  lessons  learned  from  the  failure  of  the  sabot  system  in  the  previous  test  of  July  1994.  The 
experiment  will  be  performed  using  a  conical  nosed  model  with  a  forward  facing  step  which  will  be  shot  into 
methane-oxygen-carbon  dioxide  mixture  in  the  modified  ARL  Ram  Accelerator  facility  (Refs.  13,  14).  Since 
it  was  decided  to  examine  a  more  suitable  mixture  for  the  EPA  operation  using  the  mixture  composed  of 
methane-oxygen  diluted  with  carbon  dioxide. 
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The  projectile  is  designed  to  be  aenxfynamically  stable  at  the  flight  Mach  numbers.  Preliminary  aerorfynamic 
characteristics  are  measured  in  the  supersonic  wind  tuimel  at  the  Technion  and  the  results  will  be  used  to 
finalize  the  fins  size  required  for  stable  flight  trajectory. 
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ANNEXl 


RECENT  RESULTS  WITH  THE  EXTERNAL  PROPULSION  ACCELERATOR 
J.  Rom’ ,  M.  J.  Nusca^ ,  D.  Kraczynski^  MJ.  Lewis^ ,  A.K.  Gupta^  and  J.  Sabean^ 


ABSTRACT 

The  present  p^r  presents  the  results  of  recent 
developments  in  die  research  on  the  External  Propulsion 
Accelerator  (EPA).  This  research  is  motivated  the  feet 
that  the  chemical  propellants  are  about  three  order  of 
magnitude  more  compact  in  weight  and  size  than 
electromagnetic  energy  production  and  storage  systems. 
Therefore,  there  is  great  interest  in  developing  accelerators 
of  hypervelocity  projectiles  using  chemical  propellants. 
There  are  two  methods  for  in-tube  chemical  accelerators 
utilizing  premixed  gaseous  detonative  mixtures;  the  Ram 
Accelerator  (RA)  and  the  External  Propulsion  Accelerator 
(EPA).  These  in-tube  chemical  launchers  for  accelerating 
projectiles  to  hypervelocity  utilize,  in  the  superdetonative 
mode  of  operation,  the  possibilities  of  generating 
continuous  thrust  by  initiating  detonation  in  the  premixed 
fuel/oxidizer  mixture  by  shock  wave  interactions.  The  first 
method  proposed  for  an  in-tube  chemical  launcher  was  the 
Ram  Accelerator,  originated  and  developed  by  A. 
Hert^rg  and  his  colleagues  at  the  University  of 
Washington.  The  second  method  for  operating  the 
chemical  in-tube  accelerator  is  based  on  the  utilization  of 
the  external  propulsion  cycle,  proposed  by  J.  Rom  at  the 
Technion-Israel  Institute  of  Technology.  In  the  External 
Propulsion  Accelerator  the  projectile  is  fired  into  the 
launcher  tube  which  is  filled  with  premixed  fuel/oxidizer 
mixture,  however,  here  the  projectile  diameter  is  much 
smaller  than  the  tube  diameter  (about  25%)  so  that  there 
is  no  interaction  between  the  flow  over  the  projectile  and 
the  tube  wall  over  the  complete  length  of  the  projectile.  In 
this  case  the  detonation  is  established  by  aerocfynamic 
means  on  the  projectile,  such  as  a  forward  facing  step  on 
the  projectile  shoulder  or  by  the  blunt  leading  edge  of  a 


ring  wing  positioned  on  the  center/rear  part  of  the 
projectile.  By  the  interactions  of  the  detonation  wave  with 
the  nose  shock  wave  an  “external  combustion  chamber'’  is 
produced.  This  aerodynamically  confined  region  which  is 
filled  with  the  hot  chemical  reaction  products  is  then 
e^qianded  on  the  rear  part  and  into  the  base  region  of  the 
projectile,  producing  thrust  on  the  projectile.  In  this  paper 
we  present  some  investigations  of  the  operational 
performance  limits  of  the  Ram  and  the  External 
Propulsion  Accelerators  using  an  energy  balance  analysis. 
The  energy  balance  analysis  is  applied  to  evaluate  the 
maximum  velocity  which  can  be  achieved  in  the  in-tube 
accelerators  when  the  available  chemical  reaction  energy  is 
utilized.  It  is  shown  that  the  maximum  velocity  of  the 
projectile  in  the  RA  is  limited,  due  to  the  high  drag  caused 
by  the  choking,  to  about  1.3  the  detonation  velocity  while 
the  maximum  velocity  of  the  projectile  in  the  EPA  can 
reach  up  to  6  times  the  detonation  velocity.  Then  the 
results  of  numerical  simulation  of  the  flow  over  projectile 
configurations  and  some  performance  parameters  of  the 
External  Propulsion  Accelerator  are  presented.  The 
applications  of  the  EPA  to  groimd  test  fecilities  such  as  a 
launcher  for  a  hypersonic  range  and  as  a  fecility  for 
combustion  research  in  scramjet  engines  are  discussed. 


INTRODUCTION 

The  present  paper  presents  the  results  of  recent 
developments  in  the  research  on  the  External  Propulsion 
Accelerator  (EPA).  This  research  is  motivated  ly  foe  fact 
that  foe  chemical  propellants  are  about  three  order  of 
magnitude  more  compact  in  weight  and  size  than 
electromagnetic  energy  production  and  storage  systems. 
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Therefore,  there  is  great  interest  in  developing  accelerators 
of  hypervelodty  projectiles  using  chemical  propellants. 
There  are  two  methods  for  in-tube  chemical  accelerators 
utilizing  premixed  gaseous  detonative  mixtures;  the  Ram 
Accelerator  and  the  External  Propulsion  Accelerator. 
These  in-tube  chemical  launchers  for  accelerating 
projectiles  to  hypervelodty  utilize,  in  the  superdetonathre 
mode  of  operation,  the  possibilities  of  generating 
continuous  thrust  Ity  initiating  detonation  in  the  premixed 
hiel/oxidizer  mixture  by  shock  wave  interactions.  The  first 
method  proposed  for  an  in-tube  chemical  launcher  was  the 
Ram  Accelerator,  originated  and  developed  Ity  A. 
Hertzberg  and  his  colleagues  at  the  University  of 
Washington  (Ref.  1).  The  concept  of  the  Ram  Accelerator, 
operating  in  the  superdetonative  mode,  is  based  on 
utilization  of  the  scramjet  (tycle,  where  the  projectile  acts 
as  a  fiee  centerbocty  and  the  tube  as  an  extended  cowling 
(Fig.  1).  The  sharp  nosed  projectile  diameter  is  slightly 
less  (typically  70%  to  80%  )  than  the  tube  diameter, 
therefore,  the  nose  shock  wave  is  reflected  fi'om  the  tube 
wall  into  the  projectile  centeibody.  Under  proper 
conditions  this  reflected  shock  wave  initiates  a  detonation 
process  so  that  when  the  products  of  the  chemical  reactions 
are  expanded  on  the  rear  part  of  the  projectile,  thrust  is 
generated. 

Another  method  for  operating  the  chemical  in-tube 
accelerator  is  based  on  the  utilization  of  the  external 
propulsion  cycle,  proposed  by  Rom  (Ref.  2).  In  the 
External  Propulsion  Accelerator  the  projectile  is  fiired  into 
the  launcher  tube  which  is  filled  with  premixed 
fuel/oxidizer  mixture,  however,  here  the  projectile 
diameter  is  much  smaller  than  the  tube  diameter  (about 
25%)  so  that  there  is  no  interaction  between  the  flow  over 
the  projectile  and  the  tube  wall  over  the  complete  length  of 
the  projectile  (Fig.  2).  In  this  case  the  detonation  is 
established  by  aero<tynamic  means  on  the  projectile,  such 
as  a  forward  &cing  step  on  the  projectile  shoulder  or  by  the 
blunt  leading  edge  of  a  ring  wing  positioned  on  the 
center/rear  part  of  the  projectile.  By  the  interactions  of  the 
detonation  wave  with  the  nose  shock  wave  an  “external 
combustion  chamber”  is  produced.  This  aerodynamically 
confined  region  which  is  filled  with  the  hot  chemical 
reaction  products  is  then  expanded  on  the  rear  part  and 
into  the  base  region  of  the  projectile,  producing  thrust  on 
the  projectile. 

Some  characteristics  of  the  EPA  and  various  applications 
are  discussed  in  Refs.  3,  4,  and  5.  The  CFD  calculations 
indicate  that  it  is  possible  to  establish  and  stabilize  a 
combustion-detonation  front  on  projectiles  of  various 
geometry  flying  at  hypersonic  speeds  in  detonable  gas 
mixtures.  The  characteristics  of  the  hypersonic 
combustion-detonation  on  the  projectiles  in  the  External 


Propulsion  Accelerator  are  also  relevant  to  the  Oblique 
Detonation  Wave  engine,  the  Ram  Accelerator  and  various 
scramjet  engines.  Certain  characteristics  of  the  EPA  were 
presented  in  Re&.  6,  and  7.  Some  analytical  investigations, 
based  on  many  simplifying  assunq>tions,  for  the 
establishment  of  combustion  front  ahead  of  the  forward 
facing  step  and  on  spherical  nosed  blunt  bodies  in 
hypersonic  flows  of  detonable  mixtures,  including  the 
studies  of  the  oscillations  that  mty  occur  were  investigated 
by  Tivanov  and  Rom  (Ref.  8, 9  and  10). 

In  this  paper  we  present  some  investigations  of  the 
operational  performance  limits  of  the  Ram  and  the 
External  Propulsion  Accelerators  using  an  energy  balance 
analysis.  The  energy  balance  analysis  is  applied  to  evaluate 
the  maximum  velocity  which  can  be  achieved  in  the  in¬ 
tube  accelerators  when  the  available  chemical  reaction 
energy  is  utilized.  It  is  shown  that  the  maximum  velocity 
of  the  projectile  in  the  RA  is  limited  due  to  the  high  drag 
caused  Ity  the  choking,  to  about  1.3  the  detonation  velocity 
while  the  maximum  velocity  of  the  projectile  in  the  EPA 
can  reach  up  to  6  times  the  detonation  velocity.  Then  the 
results  of  numerical  simulation  of  the  flow  over  projectile 
configurations  and  some  performance  parameters  of  the 
External  Propulsion  Accelerator  are  presented.  The 
applications  of  the  EPA  to  ground  test  facilities  such  as  a 
laimcher  for  a  hypersonic  range  and  as  a  facility  for 
combustion  research  in  scramjet  engines  are  discussed. 


AN  ENERGY  BALANCE  ANALYSIS  FOR  THE 
DETONATION  DRIVEN  PROJECTILE  IN  THE  IN¬ 
TUBE  CHEMICAL  ACCELERATORS 

The  performance  limits  of  the  in-tube  chemical 
accelerators,  the  RA  and  EPA,  can  be  evaluated  using  an 
energy  balance  analysis.  In  this  analysis  we  equate  the 
energy  of  the  chemical  reaction  with  the  drag  work  of  the 
flying  projectile.  Using  the  energy  required  to  initiate 
detonation  is  used  to  determine  the  minimtun  projectile 
velocity  which  is  needed  to  initiate  the  detonation.  While 
utilizing  the  total  energy  available  in  the  mixture  will 
determine  the  maximum  velocity  of  the  projectile. 

A.  The  Minimum  Projectile  Velocity  Required  to  Initiate 
Detonation 

The  minimum  velocity  of  the  projectile  flying  in  the  Ram 
Accelerator  (RA)  or  in  the  External  Propulsion  Accelerator 
(EPA)  needed  to  initiate  detonation  can  be  evaluated  Ity 
equating  the  energy  required  to  initiate  detonation  in  the 
fuel/oxidizer  mixture  with  the  drag  work  of  the  flying 
projectile,  as  presented  in  Ref.  12.  The  energy  for  the 
initiation  of  detonation  by  the  hypervelocity  projectile  is 
evaluated  by  the  use  of  the  blast  wave  analogy  (Ref.  12).  It 
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is  shown  in  Ref.  12  that  for  the  flow  of  the  fuel/oxidizer 
mixture  at  hypersonic  velocity  over  the  projectile  in  the 
EPA,  as  illustrated  in  Fig.  2,  the  minimiun  Mach  number 
for  the  initiation  of  detonation  on  a  projectile  with  a 
forward  feeing  step  can  be  evaluated 


H 

-^  =  2.504— 
d 


(1) 


Where  d  is  the  projectile  diameter,  the  step  height  is  H 
(where  H/d  «  1),  D  is  the  Chapman-Jouguet  detonation 
speed  and  Xr  is  the  reaction  time.  Since  ^ical  H/d  values 
can  be  1/30  and  smaller  ,  Dxr  /H  should  be  less  than  i 
(probably  in  the  range  of  1/3  to  1/5)  for  detonation  to  occur 
3°d  Conep  >  the  drag  coefficient  for  a  forward  feeing  step  is 
about  1.  Therefore,  introducing  these  values  into  Eq.  1, 
the  minimum  projectile  velocity  that  produces  energy  to 
balance  the  critical  initiation  energy  for  detonation  is 
evaluated  to  be  much  below  the  detonation  velocity.  Sinryi 
the  projectile  in  the  EPA  is  expected  to  fly  well  above  the 
detonation  q)eed  then  there  is  no  problem  of  securing  the 
energy  required  to  initiate  the  detonation  both  in  the  step 
and  also  the  ring  wing  configurations. 

This  energy  balance,  presented  in  Eq.  1,  is  also  applicable 
to  the  projectile  in  the  RA.  In  this  case,  the  drag 
coefficient  is  equated  with  that  of  the  drag  of  the  projectile 
including  the  effects  of  the  interaction  between  the 
projectile  and  the  tube  wall  and  the  reference  area  is  the 
cross  section  of  the  projectile  body.  As  the  projectile  speed 
increases  the  drag  force  on  the  projectile  increases  so  as  to 
provide  the  energy  for  initiation  of  the  detonation.  At  the 
lower  initial  spee^  the  required  drag  is  provided  by  the 
perforated  piston  orbutrator  which  is  used  in  the  RA  for 
initiation  of  the  process. 


B.  The  Maximum  Projectile  Velocity  in  the  In-Tube 
Chemical  Accelerators. 

It  is  obvious  that  the  flow  as  well  as  the  combustion 
process  are  very  different  in  the  cases  of  the  RA  and  in  the 
EPA,  however,  in  both  cases  the  work  done  by  the  drag 
force  acting  on  the  hypervelocity  projectile  is  balanced  by 
the  energy  released  by  the  detonating  gas  mixture. 


and  following  Lee  (Ref  13),  we  can  assume  a  value  of  2/3 
for  Tic  (assuming  high  combustion  temperature  of  3000°  K 
and  low  temperature  of  1000°  K).  In  2,  p  is  rfefinpd  as 
the  ratio  of  the  outside  diameter  of  the  external  combustion 
zone  to  the  projectile  diameter.  In  the  EPA  the  thickness  of 
the  external  combustion  layer  varies  as  a  function  of  the 
flow  Mach  number.  It  was  found  fiom  the  niiTnerical 
calculations,  which  are  reported  in  part  in  Refe.  3,  11  that 
for  the  32  mm  diameter  projectile  with  a  1  mm  step  the 
combustion  layer  thickness  varies  fiom  3  stq)  heights  at 
Mach  number  5  to  5.5  step  heights  at  Mach  number  6  and 
to  17  step  heights  at  Mach  number  10.  The  value  of  Cd  for 
the  cone  with  the  1  mm  step  is  about  0.26-0.28  for  Mach 
numbers  above  5  (for  cone-cylinder  without  the  step  the 
drag  coefficient  is  about  0.07-0.08).  Then,  using  Eq.  2,  the 
maximum  projectile  velocity  can  be  about  1.9  times  the 
detonation  velocity  at  Mach  munber  5,  about  2.1  times  the 
detonation  velocity  at  Mach  number  6  and  about  3.3  times 
the  detonation  velocity  at  Mach  number  10.  At  higher 
flight  Mach  numbers  the  combustion  may  extend  to  the 
tube  wall,  then  for  the  120  mm  tube  with  the  32  mm 
projectile  the  maximum  velocity  is  about  6  times  the 
detonation  velocity,  which  is  well  above  the  escape  velocity 
to  orbit  into  space.  Therefore,  when  the  projectile  is 
injected  at  initial  velocity  between  Mach  5  to  6,  only  about 
21%  to  28%  of  the  energy  of  combustion  is  used  to 
overcome  the  drag,  so  that  about  72%  to  79%  of  the 
combustion  energy  is  available  for  accelerating  the 
projectile.  Actually  the  available  combustion  energy 
increases  as  the  Mach  number  increases,  reaching  to  about 
90%  at  Mach  number  10  and  will  be  even  higher  as  the 
Mach  number  is  increased  further.  In  this  case  of 
accelerating  the  projectile  in  the  EPA  we  can  conclude  that 
only  a  small  fiaction  of  the  available  combustion  energy 
which  can  be  released  in  the  detonation  process  is  needed 
to  overcome  the  drag  of  the  hypervelocity  projectile  and 
this  available  energy  can  be  used  for  accelerating  the 
projectile  to  higher  hypervelocity.  Of  course  as  the 
projectile  velocity  increases  the  problems  of  heating  and 
ablation  of  the  projectile  surfaces  and  the  aero(tynamic  and 
acceleration  loads  on  the  projectile  structure  become  more 
critical. 


2.  Maximiun  Velocity  of  the  Projectile  in  the  RA 


1.  Maximum  velocity  of  the  projectile  in  the  EPA. 

Following  the  discussion  in  Ref.  12,  the  maximum  velocity 
is 


The  maximum  velocity  of  the  projectile  in  the  RA  is 
evaluated  using  a  similar  analysis  (Ref.  12).  Such 
evaluation  was  first  done  by  Lee  in  Ref.  13.  Lee  acRumag  a 
value  of  the  drag  coefficient  for  the  projectile  flying  in  the 
RA  tube  to  be  of  order  1,  Cd  =  0(1),  which  is  equal  to  the 
drag  coefficient  of  a  blunt  body  at  hypersonic  speeds.  So, 
using  Eq.  2  with  y  =  1.4  and  p  =  1.5,  he  estimates  the 
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value  of  the  maximum  projectile  velocity  to  be  about  1.3 
tiiTiRs  the  detonation  velod^.  He  then  concludes  that  when 
the  projectile  velocity  is  equal  to  the  detonation  velocity 
about  75%  of  the  av^able  energy  is  used  to  overcome  the 
drag  so  that  only  a  small  fraction  of  the  energy  is  available 
for  additional  acceleration  of  the  projectile.  This  result  is 
in  agreement  with  the  experimental  data  on  projectiles 
obtained  in  the  RA  of  the  University  of  Washington.  The 
projectiles  tested  at  the  University  of  Washington  and 
ARL,  as  presented  in  Refr.  1  and  11,  have  a  drag 
coefficient  of  about  0.2  -0.25  at  Mach  numbers  above  5  in 
free  flight  without  the  tube  wall  interference.  However,  the 
multiple  shock  wave  reflections  from  the  tube  wall  to  the 
projectile  center-section  and  the  blunt  leading  edges  of  the 
fins  which  may  initiate  detonation  increase  the  drag 
coefficient  of  the  projectile  considerably  in  its  flight  in  the 
accelerator  tube.  At  these  conditions  ffie  assumption  of  a 
drag  coefficient  of  about  1,  which  is  the  value  for  a  blunt 
body,  seems  reasonable.  As  the  projectile  velocity  increases 
the  strength  of  the  reflected  oblique  shock  waves  increase 
and  the  total  pressure  loss  in  the  shock  interference  region 
on  one  hand  provides  energy  for  the  initiation  of  the 
detonation  but  on  the  other  hand  due  to  the  extremely  high 
drag  consumes  an  increasing  portion  of  the  available 
reaction  energy  until  the  energy  required  to  overcome  the 
drag  force  will  be  larger  than  the  av^able  energy  from  the 
chemical  reaction  and  we  will  face  the  conditions  of 
“unstart”  if  the  projectile  will  deform  and  disintegrate  or  it 
will  exit  the  tube  at  this  maximum  velocity.  The  fact  that 
the  “unstart”  is  obtained  at  about  1.15  to  1.2  the  detonation 
velocity  is  shown  in  Figs.  3a  and  3b,  based  on  the  data  of 
Ref  14,  for  aluminum  projectiles.  Using  a  titanium 
projectile,  which  did  not  deform,  the  maximum  velocity  is 
obtained  and  the  projectile  flies  at  this  maximum  velocity 
for  the  last  3  meters  of  the  tube,  as  shown  in  Fig.  4.  This 
maximum  velocity  is  1.2  the  detonation  velocity  and  is  in 
very  good  agreement  with  the  energy  analysis  value. 


THE  PERFORMANCE  OF  PROJECTILES  FLYING  AT 
HYPERSONIC  SPEED  IN  THE  EXTERNAL 
PROPULSION  ACCELERATOR 

A.  Description  of  the  External  Propulsion  Accelerator 
(EPA) 

An  outline  of  an  EPA  facility  is  presented  in  Fig.  5. 

The  first  stage  of  the  EPA  fecility  can  be  a  large  caliber 
smooth  bore  conventional  powder  gun  for  the  initial 
acceleration  of  the  sabot-projectile/missile  payload.  Since 
the  projectile/missile  must  be  accelerated  to  velocities 
higher  than  the  detonation  velocity  of  the  fiiel/oxidizer 
mixture  in  the  accelerator  tube  (detonation  speeds  of  1000 
m/sec  to  1500  m/sec),  we  require  muzzle  exit  velocities  of 


1300  m/sec  to  1800  m/sec  in  order  to  accommodate 
various  mixture  con:qx)sitions.  These  muzzle  velocities  can 
be  achieved  in  large  caliber  powder  guns  with  barrel 
lengths  offimto  10  m 

The  second  stage  most  be  a  sabot  separation  section.  There 
are  various  methods  used  for  sabot  separation,  mechanical 
and/or  aerodynamic  methods,  which  can  be  utilized  in  this 
design.  In  Fig.  5  a  sabot  stripper  plate  is  used.  The  “clean” 
piojectile/missile  flies  throng  the  sabot  stripper  and  enters 
at  a  velocity  which  is  higher  than  the  detonation  velocity 
into  the  fiiel/oxidizer  mixture  in  the  accelerator  tube. 

The  third  stage  is  the  accelerator  section.  This  section  is  a 
large  tube,  its  diameter  should  be  in  the  order  of  4  times 
the  projectile/missile  maximum  diameter,  which  is  filled 
with  various  piemixed  gaseous  fuel  and  oxidizer  mixtures. 
Both  ends  of  the  accelerator  tube  must  be  sealed,  either  by 
diaphragms  which  are  pierced  by  the  flying  missile  or  by 
quick  acting  valves  which  are  opened  in  synchronization 
with  the  flying  missile.  This  tube  can  be  divided  into 
sections  closed  by  separating  diaphragms,  so  that  different 
mixtures  can  be  filled  into  each  of  these  subsections.  In 
this  way  mixtures  of  increasing  detonation  speeds  can  be 
matched  to  the  projectile/missile  velocity  as  this  velocity 
increases  due  to  the  acceleration  of  the  projectile/missile  as 
it  flies  down  the  tube  through  these  sub-sections. 

The  tube,  in  this  case,  acts  only  as  a  containing  vessel  for 
the  piemixed  fiiel/oxidizer  atmosphere  through  which  the 
projectile/missile  flies.  The  tube  structure  must  be  designed 
to  withstand  the  initial  mixture  pressure  which  can  be  of 
the  order  of  20  atm.  to  200  atm.,  if  we  wish  to  use  the  tube 
for  a  single  shot  operation,  or  it  can  be  made  to  withstand 
the  peak  pressures  experienced  after  the  passage  of  the  gun 
blast  and  the  explosion  of  the  remaining  mixture  after  the 
passage  of  the  projectile/missile. 

At  the  exit  of  the  accelerator  section  the  projectile/missile 
emerges,  after  piercing  the  exit  diaphragm,  at  its  highest 
velocity.  The  exit  velocity  is  determined  by  the  thrust 
which  is  generated  on  the  projectile/missile  and  the  length 
of  the  accelerator  section.  The  calculations  reported  in  Ref. 
3  indicate  that  thrust  levels  of  F/pA  (where  F  is  the  thrust 
force,  p  is  the  initial  pressure  in  the  tube  and  A  is  the 
maximum  cross  section  of  the  projectile/missile)  of  about  3 
can  be  achieved.  This  thrust  level  can  be  increased  by 
using  more  energetic  fuels  and  by  optimizing  the 
projectile/missile  design.  This  thrust  level  is  sufficient  to 
accelerate  reasonable  size  projectile/missile  to  above 
50,000  g’s.  We  use  this  acceleration  level  since  there  are 
guidance  and  control  systems  as  well  as  telemetering  units 
hardened  to  operate  at  this  50  Kg’s  level.  The  added 
velocity  at  the  exit  of  the  accelerator  tube  as  a  function  of 
the  tube  length  is  shown  in  Fig.  6.  Thus,  assuming  initial 
velocity  of  1500  m/sec.,  the  exit  velocity  for  a  8  m  long 
accelerator  tube  will  be  3500  m/sec.  and  for  32  m  long 
tube  the  exit  velocity  will  be  5500  m/sec.  It  was  shown  in 
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Ref.  12,  using  an  energy  balance  analysis,  that  it  is 
possible  to  reach  in  the  EPA  a  maximum  velocity  of  about 
6  times  the  detonation  velocity  of  the  mixture  used  in  the 
tube.  Therefore,  using  mixtures  with  detonation  speeds  of 
2,500  m/sec.  and  up  to  3,000  m/sec.  enable  the  utilization 
of  the  EPA  to  reach  about  18,000  m/sec.,  well  bq^ond 
escape  velocity  from  earth.  This  may  be  a  veiy  promising 
method  to  obtain  a  “single  stage  to  orbit”  mission. 

The  diameter  of  the  acceleration  section  is  about  1  m  for  a 
25  cm  diameter  projectile/missile  and  40  cm.  for  a  10  cm. 
diameter  projectile/missile,  or  smaller  or  larger  for  smaller 
or  larger  projectile/missiles. 

B.  Aerodynamic  design  of  the  EPA  projectile. 

After  being  stripped  from  its  s^t  and  upon  entering  the 
acceleration  tube  section  the  projectile  flies  freely  in  the 
gas  mixture  which  is  in  the  tube.  The  projectile  is  aimed  to 
pierce  the  entrance  diaphragm  near  its  center  and  should 
fly  near  the  tube  center.  Therefore,  the  projectile  is 
designed  to  be  aerodynamically  stable,  both  statically  and 
dynamically.  In  order  to  insure  stability  over  the  range  of 
flight  Mach  numbers  the  basic  configuration  needed  for 
external  propulsion,  i.e.  the  cone,  forward  facing  step  « 
cylindricd  with  boattail  aft  section,  is  fitted  with  small 
elongated  fins  at  the  end  of  the  boattail  section.  At  first,  a 
design  with  fins  with  a  span  equal  to  the  maximum  body 
diameter  is  investigated.  This  configuration  with  6  fins  is 
shown  in  Fig.  7.  In  order  to  stabilize  the  flight  trajectory  in 
the  center  of  the  tube  the  possibility  of  imparting  a  slow 
spin  by  canting  the  fins  is  also  tested.  In  Fig.  7  the  fins  are 
canted  to  3®  to  the  model  axis  while  the  tail  with  parallel 
fins  is  also  shown.  The  aerodynamic  coeflBcients  are 
measured  in  the  supersonic  wind  tunnel  of  the 
Aerodynamic  Laboratoiy  in  the  Technion.  The 
aerodynamic  coefficient  are  measured  at  M  =  3.4  and  these 
values  are  extrapolated  to  the  hypersonic  flight  Mach 
number  using  slender  body  analysis  methods.  Preliminary 
wind  tunnel  results  indicate  that  the  aerodynamic  center  of 
the  6  fins  configuration  is  positioned  so  as  to  obtain 
reasonable  stability  margin  for  the  positions  of  the  center 
of  gravity  of  the  projectile.  The  aerodynamic  coefficient 
evaluated  from  the  wind  tunnel  results  are  used  to  evaluate 
the  flight  trajectoiy  by  a  six  degrees  of  freedom  program. 
This  work  is  currently  underway  and  will  be  included  in 
future  report.  In  parallel,  a  firing  test  is  being  planed  using 
the  Ram  Accelerator  facility  at  the  Army  Research 
Laboratories  of  the  US  Army  at  Aberdeen  Proving 
Grounds  in  Maryland,  USA. 

Following  this  investigation  on  the  aerodynamic 
stabilization  and  flight  trajectoiy  it  is  planned  to  study  the 
heat  transfer  rates  to  the  projectile  surfaces  and  determine 
a  practical  method  for  thermal  protection  at  the  hypersonic 
speeds  that  can  be  achieved  by  projectiles  in  the  EPA. 


C.  Evaluation  of  the  combustion  characteristics  on  the 
projectile. 

The  performance  characteristics  of  the  projectile  flight  in 
the  accelerator  tube  is  evaluated  by  the  use  of  the  computer 
code  which  is  based  on  the  solution  of  the  Navier-Stokes 
equations  for  flows  with  chemical  processes.  This 
computer  code  was  developed  for  the  Ram  Accelerator 
calculations  at  ARL  by  Nusca  (Refs.  16,  17).  This  program 
uses  the  Rockwell  Science  Center  USA-RG  (Unified 
Solution  Algorithm  Real  Gas)  code  written  \yy 
Chakravarthy  et  al  (Refs.  18,  19).  This  CFD  code  includes 
the  full  3D  unsteady  Reynolds-Averaged  Navier-Stokes 
(RANS)  equations  including  equations  for  chemical 
kinetics  (finite-rate  and  equilibrium).  These  equations  are 
cast  in  conservation  form  and  converted  to  algebraic 
equations  using  upwind  finite-difference  and  finite-volume 
formulations.  The  equations  are  solved  using  a  second- 
order  TVD  (total  variation  diminishing)  scheme  which  is 
used  to  insure  non-oscillatoiy  numerical  behavior.  The 
flow  field  with  the  combustion  generated  by  the  forward 
facing  step  is  presented  in  Fig.  8  and  for  the  model  with 
the  ring  wing  in  Fig.  9.  Following  these  calcidations,  test 
runs  of  the  forward  facing  step  configurations  and  the  ring 
wing  projectiles  are  being  planned  in  the  ARL  Ram 
Accelerator  fticility  (Refs.  20,  21)  using  subcaliber  models. 

USE  OF  THE  EXTERNAL  PROPULSION 
ACCELERATOR  AS  TEST  FACILITIES 

A.  HYPERSONIC  FLIGHT  RANGE 
The  External  Propulsion  Accelerator  (EPA)  can  be  used  as 
the  launcher  of  aerodynamic  models  to  hypersonic  speeds 
in  a  ground  test  facility.  The  models,  encased  in  a  properly 
designed  sabot,  are  accelerated  in  the  EPA  to  the  desired 
hypersonic  speed  and  fired  into  an  instrumented  ballistic 
range.  The  EPA  is  theoretically  capable  to  launch 
projectiles  to  velocities  of  up  to  6  times  the  detonation 
velocity  of  the  fuel-oxidizer  mixture,  in  the  order  of  12,000 
m/sec  to  18,000m/sec,  so  that  it  is  well  suited  to  be  a 
laimcher  for  the  hypersonic  ballistic  range.  The  free  flight 
range  technology  seems  to  enable  the  best  simulation  of 
atmospheric  free  flight  conditions  in  a  groimd  based 
facility.  The  range  concept  is  most  suitable  since  it  uses 
quiescent  arbitrary  gas  in  its  test  sections  in  which  the 
model  flies  at  the  design  hypersonic  speed.  While  in  other 
test  facilities,  such  as  the  various  types  of  shock  tunnels 
and  gun  tunnels,  it  is  required  to  compress  and  heat  the 
test  gas  to  high  stagnation  temperatures  and  pressures 
followed  by  a  rapid  expansion  in  high  Mach  number 
nozzles.  So  that  in  such  facilities  the  simulation  is 
distorted  by  the  effects  of  dissociation  and  ionization  of  the 
gas  in  the  shock  tunnel  stagnation  chamber  as  well  as 
disturbances  due  to  nozzle  flow  non-uniformities,  effects  of 
boundaiy  layer  disturbances,  etc.  The  test  range  can  be 
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divided  into  sections  separated  diaphragms,  enabling 
simulation  of  flight  at  various  pressure  altitudes  in  a  single 
shot.  Therefore,  such  a  hypersonic  range  with  the  EPA 
laimcher  can  be  an  excellent  tool  for  hypersonic  aircraft 
research  and  development 

The  requirements  for  a  hypersonic  ground  based  flight  test 
range  and  the  feasibility  of  the  use  of  light-gas  gun 
launcher  was  discussed  by  Witcofiki  et  al  (Ref  22).  In  this 
study  the  use  of  an  Electromagnetic  launcher  (EML)  and 
of  the  Ram  Accelerator  (RA)  where  also  considered.  It  was 
concluded  that  with  the  present  state  of  technology  the 
most  suitable  candidate  is  the  two-stage  light-gas  gun 
which  can  be  made  to  launch  models  of  about  20  cm.  span 
and  length  of  up  to  90  cm.  to  speeds  of  up  to  6  km/sec.  It 
was  contemplated,  that  at  a  later  stage,  using  more 
advanced  launcher  technology,  the  models  can  be 
increased  to  30  cm  to  45  cm  in  span  and  the  velocity  range 
will  be  increased  to  10  km/sec  to  15  km/sec. 

THE  TEST  FACILITY 

In  this  application,  the  hypersonic  testing  fecility  is 
basically  similar  to  the  AHAF  concept  presented  in  Ref 
22.  The  &cility  will  be  comprised  from:  1.  The  Accelerator 
section.  2.  Sabot  separation  section.  3.  Ballistic  Range  for 
aerodynamic  models  free  flight  testing.  4.  Decelerator 
section  for  model  retrieval.  An  artist  view  of  such  facility, 
reproduced  from  Ref.  22,  is  shovm  in  Fig.  10. 

1.  The  Accelerator  Section 

The  accelerator  is  comprised  of  two  parts,  an  initial 
accelerator  for  accelerating  the  model  with  its  sabot  to  the 
required  insertion  speed  of  1,400  to  1,800  m/sec.  (  Mach 
number  between  5  to  6.6)  and  the  External  Propulsion 
launcher  tube  which  will  accelerate  the  projectile,  in  which 
the  model  is  inserted,  into  the  test  range.  The  flight  Mach 
Number  in  the  ballistic  range  can  be  then  set  to  the 
required  hypersonic  speed  with  the  possibility  of  reaching 
M  =  30.  The  type  of  initial  accelerator  and  the  length  of 
the  External  Propulsion  Launcher  depend  on  the  allowable 
acceleration  for  the  model.  Simple  aerodynamic  shapes 
used  for  hypersomc  flight,  such  as  blunt  reentry  bodies, 
slender  shuttle  type  configurations  and  wave  riders  with 
hardened  on  board  instrumentation  can  be  built  to  operate 
up  to  50,000  g’s  acceleration  level.  In  this  case,  a  gun 
powder  laimcher  can  be  used  for  the  initial  acceleration 
and  the  acceleration  level  of  the  External  Propulsion 
launcher  can  be  adjusted  to  the  desired  level.  The  gun 
launcher  will  require  lengths  of  5  to  20  meters.  This 
section  ivill  be  followed  by  the  External  Propulsion 
launcher.  The  lengths  of  the  EPA  section  will  depend  on 
the  achievable  and  allowable  acceleration  levels.  Using  the 
level  of  50,000  g’s,  the  required  lengths  of  the  EPA  (Fig. 
6)  wiU  be  about  2  meters  for  M  =  7  testing,  8  meters  for  M 
=  11  testing,  about  18  meters  for  M  =  15  testing  and  for 
Mach  number  30  testing  a  launcher  length  of  about  120 


meters  will  be  required.  Using  the  conservative  evaluation 
for  the  thrust  coefficient,  F/pA  =  3,  the  thrust  generated  on 
a  25  cm  diameter  projectile  is  nearly  300,000  Kg.  for 
initial  pressure  of  50  atm.  The  projectile  weight,  including 
the  model  and  sabot,  is  then  about  58  Kg.  for  the  50,000 
g’s  acceleiatioa  This  weight  can  be  doubled  if  the  initial 
pressure  is  increased  to  100  atm.  enabling  robust  structure 
for  the  model  and  sabot  as  well  as  for  the  onboard 
instrumentation.  The  initial  accelerator  as  well  as  the  EPA 
can  be  “cushioned”  for  lower  acceleration  levels.  The 
length  of  this  launcher  varies  inversely  as  the  acceleration 
levels,  i.e.  lower  acceleration  will  require  the 
correspondingly  longer  launcher  lengths. 

The  External  Propulsion  launcher  diameter  should  be 
about  4  times  the  model  maximum  diameter  or  span. 
Therefore,  a  launcher  of  Im  diameter  will  enable 
launching  of  models  of  25  cm.  span  encased  in  a  projectile 
of  comparable  diameter.  The  launcher  tube  diameter  can 
be  reduced  to  40  cm  for  10  cm  projectile  diameter.  Of 
course  larger  models  can  be  considered  with  increase  in 
the  diameter  of  the  accelerator. 

2.  Blast  Chamber  and  Sabot  Separation  Section. 

The  blast  of  the  combusting  gas  mixture  as  well  as  the 
blast  of  the  gun  for  the  initial  acceleration  follow  the 
projectile  and  should  be  retained  in  a  blast  chamber  at  the 
exit  of  the  launching  tube.  This  chamber  is  followed  by  the 
sabot-model  separating  section,  as  shown  in  Fig.  10. 

Since  the  projectile  shape  required  for  the  initial 
accelerator  and  for  the  External  Propulsion  Accelerator 
should  be  of  a  specified  axi-symmetric  shape  the 
aeroilynamic  model  must  be  encased  in  a  sabot.  The 
outside  shape  of  this  sabot  should  be  in  accordance  with 
the  External  Propulsion  projectile  design  while  supporting 
the  model  to  withstand  the  acceleration  levels.  Upon  exit 
from  the  accelerator  launcher  the  sabot  should  be  separated 
either  mechanically  or  by  special  mechanism  which  will  be 
initiated  at  the  exit  of  the  launcher  and  the  sabot  parts 
should  be  retained  in  this  separation  section  while  the  free 
model  flies  into  the  test  range. 

3.  Ballistic  Range  Section. 

The  free  model  can  now  fly  in  the  instrumented  ballistic 
range  where  its  instantaneous  trajectory  positions  and 
attitudes  can  be  photographed  and  measured.  Further 
ofiboard  instrumentation  may  include  optical 
instrumentation  for  measurements  of  the  vibrational, 
translational  and  electron  temperatures  and  interferometric 
photographs  for  measurements  of  density  profiles.  In 
addition  various  onboard  instrumentation  can  be  installed 
in  the  model,  such  as  wall  temperatures  and  pressures, 
accelerometers  for  model  acceleration,  strain  gages  for 
loads  and  stresses,  and  these  can  be  monitored  by 
telemetering  system.  If  we  consider  models  of  10  cm  to  25 
cm  in  diameter  or  span,  such  models  are  sufficiently  large 
to  include  most  details  of  the  aerodynamic  design. 
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including  deflected  control  surfeces  and/or  small  control 
jets  actuated  by  small  rockets  or  gas  supply. 

The  ballistic  range  can  either  be  an  atmospheric  range  or 
be  enclosed  in  a  large  tube  enabling  controlled  atmosphere 
for  simulation  of  pressure-altitude  effects  and  also 
atmospheric  composition  of  other  planets.  The  enclosed 
test  range  in  a  large  tube  can  be  divided  into  sections 
separated  by  diaphragms,  enabling  simulation  of  flight  at 
various  pressure  altitudes  and/or  various  gases  in  a  single 
shot. 

4.  Deceleration  Section. 

The  energy  imparted  to  the  model  traveling  at  Mach  6  to 
30  is  very  large  and  must  be  absorbed  in  order  to  stop  the 
model.  For  simple  models,  the  models  can  be  expandable 
and  a  conventional  “catcher”  can  be  used.  For  more 
sophisticated  and  expansive  models  there  may  be  some 
possibilities  of  some  energy  absorbing  methods  for 
arresting  the  models,  at  least  in  the  lower  hypersonic  speed 
range. 

B.  THE  HYPERSONIC  SCRAMJET  TEST  FACILITY 
The  combustion-detonation  phenomena  in  a  scramjet 
engine  as  well  as  in  an  Oblique  Detonation  Wave  Engine 
can  be  studied  in  the  External  Propulsion  Accelerator.  This 
facility  provides  the  means  to  investigate  the  combustion- 
detonation  process  independently  from  the  fuel  mixing 
process.  In  this  case,  a  projectile  simulates  the  engine 
geometry  and  the  combustion  of  the  premixed  fiiel/oxidizer 
mixture  is  initiated  in  the  model  flying  in  the  External 
Propulsion  Accelerator.  The  scramjet  combustion  also 
provides  net  thrust  for  the  model  flight  in  the  Accelerator. 
Thus,  this  facility  enables  testing  of  the  combustion 
characteristics  as  function  of  inlet,  combustor  and  nozzle 
geometry,  fuel/oxygen  chemical  parameters  and  flight 
velocity.  The  effects  of  various  fiiel/oxidizer  mixtures  can 
be  studied  by  including  sections  separated  by  diaphragms 
in  the  accelerator  tube.  These  sections  can  be  filled  with 
various  mixtures  to  different  initial  pressures,  so  that 
effects  of  various  pressure  levels  can  also  be  simulated  in 
this  facility. 

The  scramjet  model  will  be  composed  of  a  shallow  sharp 
nosed  cone  acting  as  the  scramjet  inlet,  a  ring  wing 
simulating  the  cowling  of  the  scramjet  around  a  cylindrical 
center  body  and  a  conical  afterbody  simulating  the 
scramjet  exhaust.  Combustion  is  generated  behind  the 
shock  wave  established  in  front  of  the  wedge  shape  lip  of 
the  ring  wing.  The  flow  field  is  illustrated  in  Fig.  9.  The 
scramjet  combustion  in  this  configuration  will  produce 
high  temperature  high  pressure  combustion  products 
resulting  in  significant  net  positive  thrust. 

THE  SCRAMJET  TEST  FACILITY 

In  this  application,  the  hypersonic  scramjet  testing  facility 

will  be  comprised  from:  1.  The  External  Propulsion 


Accelerator  section.  2,  Transparent  test  section.  3. 
Deccelerator  section  for  model  retrieval. 

1.  The  External  Propulsion  Accelerator  Section 
The  tests  are  conducted  in  a  regular  External  Propulsion 
Accelerator  using  a  specially  designed  model  for  the 
scramjet  simulation.  The  model  with  its  sabot  are 
accelerated  by  a  gun  powder  charge  to  the  required 
insertion  speed  of  1,000  up  to  1,900  m/sec.  ( Mach  number 
between  3  to  5).  In  a  simple  configuration  the  projectile 
body  will  simulate  the  engine  centeibody  and  a  properly 
shaped  ring  wing  around  this  centeibody  will  represent  the 
engine  cowling.  The  shapes  of  the  engine  intake  and  the 
shape  of  the  combustion  region  inside  the  engine  can  be 
also  simulated  by  this  model.  Combustion-detonation  will 
be  initiated  by  the  shock  waves  generated  at  the  intake  to 
the  ring  wing  and  it  will  cause  combustion  in  the  space 
between  the  centeibody  and  the  ring-cowling.  The  flow 
will  be  then  expanded  on  the  aft-body  of  the  model, 
simulating  the  engine  exhaust  and  nozzle  flow 
characteristics.  The  operation  of  the  scramjet  will  produce 
net  thrust  for  the  acceleration  of  the  model.  Therefore, 
measurements  of  the  model  acceleration-deceleration  will 
enable  the  evaluation  of  the  thrust  obtained  by  this 
scramjet  engine  as  a  function  of  model  velocity.  More 
sophisticated  measurements  can  be  obtained  using  on¬ 
board  instrumentation  with  telemetering.  Models  with 
hardened  on-board  instrumentation  can  be  built  to  operate 
up  to  50,000  g’s  acceleration  level.  The  acceleration  level 
in  the  projectile  will  depend  on  the  thrust  generated  by  the 
scramjet.  Since  the  gas  mixtures  can  be  varied  in  the 
accelerator,  this  will  be  a  method  to  study  the  effects  of 
various  chemical  fuel  compounds  on  the  combustion  and 
the  thrust  generation.  Such  studies  can  be  made  when  the 
accelerator  tube  is  divided  into  a  number  of  sections 
separated  by  diaphragms.  Each  section  of  the  accelerator 
tube  can  be  filled  with  a  certain  fiiel/oxidizer  mixture  at 
various  initial  pressures.  So  that  in  a  single  experiment  the 
combustion  characteristics  of  flight  in  various  mixtures  at 
specified  initial  pressures  can  be  studied  These  test 
sections  can  be  positioned  at  the  end  of  the  basic 
accelerator  tube  in  which  the  projectile  is  accelerated  to  the 
desired  velocity. 

The  length  of  the  facility  will  be  determined  by  the 
acceleration  of  the  model  and  the  speed  range  for  the  tests. 
For  the  External  Propulsion  launcher  at  50,000  g’s,  the 
required  lengths  of  the  EPA  (Fig.  6)  will  be  about  2  meters 
for  M  =  7  testing,  8  meters  for  M  =  11  testing,  about  18 
meters  for  M  =  15  testing  and  for  Mach  number  30  testing 
a  launcher  length  of  about  120  meters  will  be  required. 
More  sophisticated  models  requiring  lower  values  of 
acceleration  can  be  used  by  proper  modifications  to  the 
launcher.  Lower  initial  accelerations  can  be  achieved  by 
initial  acceleration  of  either  a  Ram  Accelerator  stage 
initiated  by  a  small  light  gas  gun  or  a  small  rocket  or  the 
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initial  acceleration  be  achieved  direcUy  a  rocket  motor 
which  will  separate  at  the  entrance  to  the  External 
Propulsion  launcher.  The  length  of  this  launcher  varies 
inversely  as  the  acceleration  levels,  i.e.  lower  acceleration 
will  require  the  correspondingly  longer  lengths. 

The  External  Propulsion  launcher  diameter  should  be 
about  4  times  the  model  maximum  diameter.  Therefore,  a 
launcher  of  Im  diameter  will  enable  launching  of  models 
of  25  cm.  diameter  and  the  laimcher  tube  diameter  can  be 
reduced  to  40  cm  for  10  cm  models.  Of  course  larger 
models  can  be  considered  with  increase  in  size  of  the 
facility. 

2.  Transparent  Test  Section. 

The  possibility  of  including  a  transparent  section  into  the 
Ram  Accelerator  facility  was  studied  at  ARL  by 
Kruczynski  et  al.  (Ref  15).  Such  a  transparent  section  at 
the  end  of  the  launcher  tube  will  enable  flow  visualization 
photographs  on  the  model,  particularly,  at  the  intake  and 
exhaust  regions. 

3.  Deceleration  Section. 

The  energy  imparted  to  the  model  traveling  at  Mach  6  to 
15  is  very  large  and  must  be  absorbed  in  order  to  stop  the 
model.  For  simple  models,  the  models  can  be  expandable 
and  a  conventional  “catcher”  can  be  used.  For  more 
sophisticated  and  expansive  models  there  may  be  some 
possibilities  of  some  energy  absorbing  methods  for 
arresting  the  models,  at  least  in  the  lower  hypersonic  speed 
range. 

SUMMARY  AND  CONCLUSIONS 

The  two  methods  for  in-tube  chemical  accelerators 
utilizing  premixed  gaseous  detonative  mixtures;  the  Ram 
Accelerator  (RA)  and  the  External  Propulsion  Accelerator 
(ERA)  are  examined  using  an  energy  balance  analysis.  It  is 
shown  that  the  maximum  velocity  of  the  projectile  in  the 
RA  is  limited,  due  to  the  high  drag  caused  by  the  choking, 
to  about  1.3  the  detonation  velocity  while  the  maximiun 
velocity  of  the  projectile  in  the  EPA  can  reach  up  to  6 
times  the  detonation  velocity. 

The  results  of  numerical  simulation  of  the  flow  over 
projectile  configurations  using  a  Navier-Stokes  code  with 
chemical  reactions  show  that  combustion  can  be  generated 
by  a  small  step  on  the  shoulder  of  the  projectile.  These 
calculations  indicate  that  positive  thrust  is  obtained. 

It  is  shown  that  the  projectile  can  be  stabilized 
aerodynamically  by  6  small  fins  which  can  be  canted  to 
impart  roll  to  the  flying  projectile. 

The  applications  of  the  EPA  to  ground  test  facUities  such 
as  a  launcher  for  a  hypersonic  flight  range  and  as  a  facility 
for  combustion  research  in  scramjet  engines  are  discussed. 
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Figure  1.  The  Ram  Accelerator  projectile  in  the  detonative 
mode  (Ref.  1) 


Figure  2.  The  External  Propulsion  Accelerator  projectile 
with  a  forward  facing  step  (Ref.  2) 
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Figure  3a.  Results  of  variable  nitrogen  experiments  plotted 
as  thermally  choked  heat  release  vs.  flight  Mach 
number  (from  Ref.  14) 


Ground-Based  Flight-Test  Range”,  AIAA  Paper  91-0296, 
1991 
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Figure  3b.  Results  of  variable  methane  experiments  plotted 
as  thermally  choked  heat  release  vs.  flight  Mach 
number  (from  Ref.  14) 
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Figure  4.  Velocity-distance  data  comparing  aluminum  and 
titanium  nose  projectiles  (from  Ref.  14) 
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Figure  5.  Schematic  of  the  External  Propulsion  Accelerator 
facility 


Figure  7.  An  External  Propulsion  Accelerator  projectile 
with  a  1  mm  step  and  6  fins  tail 
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Figure  6.  The  added  velocity  as  a  function  of  accelerator 
tube  length  for  50,000  g’s  acceleration 
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Figure  8.  The  temperature  ratio  distribution  around  the  1 
mm  step  model  for  Mach  number  6  flow  with 
chemical  reactions  in  the  EPA. 


Figure  9.  The  temperature  ratio  distribution  around  the 

ring  wing  model  for  Mach  number  6  flow  without 
chemical  reactions  in  the  EPA. 
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Figure  10.  Schematic  view  of  the  hypersonic  flight  range 
facility  using  EPA  launcher 
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INVESTIGATION  OF  THE  COMBUSTION  INDUCED  BY  A  STEP  ON  A 
PROJECTILE  FLYING  AT  HYPERSONIC  SPEED  IN  AN 
EXTERNAL  PROPULSION  ACCELERATOR 

J.  Rom^j-M.J.  Nusca^,  D.L.  Kruczynski^,  M.J.  Lewis^,  A.K.  Gupta®,  and  J.  Sabean® 


Abstract 

The  establishment  of  combustion  by  a  small  forward 
facing  step  on  a  hypervelocity  projectile,  in  the  External 
Propulsion  Accelerator,  is  demonstrated  via  a  computa¬ 
tional  fluid  dynamics  simulation.  Computations  of  the 
flow  over  the  projectile  flying  at  hypersonic  speeds  in  a 
combustible  mixture  are  presented.  The  computations 
indicate  that  a  small  forward-facing  step,  of  the  order 
of  1  mm,  can  initiate  detonation  at  the  step  and  sus¬ 
tain  combustion  over  the  afterbody  and  blunt  base  of 
the  projectile,  flying  at  Mach  numbers  above  5.  The  in¬ 
teraction  between  the  shock  wave  from  the  shallow  cone 
forebody  and  the  detached  shock-detonation  wave  from 
the  forward- facing  step,  positioned  on  the  shoulder  of 
the  projectile  (i.e.,  the  forebody-afterbody  junction)  re¬ 
sults  in  confinement  of  a  combustion  layer  on  the  exter- 
mi  surface  of  the  projectile  afterbody  and  on  the  blunt 
base.  The  high  pressure,  high  temperature  products  of 
the  combustion,  which  are  confined  to  the  external  com¬ 
bustion  layer,  are  then  expanded  over  the  afterbody 
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and  into  the  wake  of  the  projectile  and  can  be  used 
to  generate  thrust.  Thus  this  layer,  confined  by  the 
aerodynamic  interactions,  acts  as  an  external  combus¬ 
tion  chamber.  This  is  an  excellent  demonstration  of 
the  external  propulsion  method.  Some  preliminary 
computed  results  are  presented  for  a  32mm  projectile 
flying  in  a  120mm  accelerator  tube,  at  Mach  num¬ 
bers  of  6  and  10.  The  tube  is  filled  with  a  fuel-rich 
methane/oxygen/nitrogen  mixture  at  an  initid  pressure 
of  50  atm.  The  results  show  that  for  projectiles  with  step 
heights  of  0.5, 1,  and  1.5  mm,  combustion  is  established 
and  maintained. 

Introduction 

The  question  of  establishment  and  stabilization  of  a 
combustion  or  detonation  front  on  a  projectile  flying  at 
hypersonic  speeds  in  a  detonable  mixture  is  of  inter¬ 
est  both  as  a  fundamental  combustion  problem  and  for 
practical  application  in  the  External  Propulsion  Accel¬ 
erator.  It  is  also  relevant  to  the  Oblique  Detonation 
Wave  (ODW)  engine,  the  Ram  Accelerator,  and  various 
scramjet  engines.  The  External  Propulsion  Accelerator 
was  proposed  by  J.  Rom^  and  certain  aspects  of  its  char¬ 
acteristics  were  presented  in  References  2  and  3.  Some 
aspects  of  the  establishment  of  a  combustion  front  ahead 
of  the  forward-facing  step  on  a  hypervelociy  projectile, 
and  on  a  spherically  blunted  body  in  hypersonic  flow  of 
a  detonable  mixture  (including  oscillations)  were  inves¬ 
tigated  by  Tivanov  and  Rom.^»®»® 

The  External  Propulsion  Accelerator  has  been  devel¬ 
oped  following  that  of  the  Ram  Accelerator  by  Hertzberg 
and  colleagues  at  the  University  of  Washington.^  The 
Ram  Accelerator  is  based  on  the  idea  of  using  a  pre¬ 
mixed  fuel/oxidizer/diluent  mixture  in  a  sealed  launch 
barrel,  and  injecting  a  projectile  at  supersonic  speeds 
into  the  mixture.  Thrust  is  generated  by  the  ramjet 
combustion  obtained  between  the  tube  wall  and  the  pro¬ 
jectile,  where  the  projectile  acts  as  a  ramjet  centerbody 
and  the  barrel  is  basically  a  long  cowl.  In  the  super¬ 
sonic  mode  combustion  is  initiated  by  the  interaction  of 
shock  waves,  generated  at  the  projectile  nose,  and  re¬ 
flected  between  the  tube  walls  and  the  projectile  body. 


Unlike  the  Ram  Accelerator,  for  the  External  Propul¬ 
sion  Accelerator,  thrust  is  generated  by  combustion  ini¬ 
tiated  by  a  normal  shock  wave  generated  on  the  projec¬ 
tile  body  by  a  ramp  or  forward-facing  step  (or  by  a  ring 
wing^)  independent  of  the  barrel  walls.  Since  the  pro¬ 
jectile  flies  in  a  premixed  fuel/oxidizer  atmosphere,  the 
problem  of  mixing  fuel  and  oxidizer,  which  plagues  air- 
breathing  propliision  methods  proposed  for  hypersonic 
flight,  is  eliminated.  However,  there-are  remaining  dif- 
flcult  issues  such  as  the  ability  to  initiate  and  stabilize 
the  combustion  front  on  the  projectile  in  such  a  way  that 
positive  thrust  can  be  generated.  Previous  studies^t^  in¬ 
dicated  that,  if  a  relatively  small  forward-facing  step 
on  the  projectile  shoulder  (i.e.  the  forebody-afterbody 
junction)  is  able  to  induce  a  detonation  wave  and  a  sus¬ 
tained  combustion  front,  then  reasonable  thrust  levels 
may  be  achieved. 

Therefore,  the  present  investigation  will  study  the 
ability  of  a  small  forward-facing  step  to  ignite  the 
premixed  fuel/oxidizer  mixture  flowing  at  hypersonic 
speeds  over  the  projectile,  and  study  the  flow  field  and 
combuston  front  generated  by  the  interaction  of  the 
shock  wave  from  the  projectile  nose,  and  the  detonation 
w’ave  generated  by  the  step.  This  interaction  results  in 
the  confinement  of  a  high  pressure,  high  temperature 
region  below  a  contact  surface.  This  region,  between 
the  projectile  body  and  the  contact  surface,  acts  like  the 
combustion  chamber  of  an  engine.  However,  in  this  case, 
the  contact  surface  acts  as  the  “engine  cowling”  with¬ 
out  a  structural  boundary.  This  is  an  example  of  a  true 
external  propulsion  “combustion  chamber”.  The  high 
temperature  and  pressure  combustion  products  are  then 
expanded  over  the  afterbody  and  the  base  of  the  projec¬ 
tile  to  produce  thrust.  The  flow  field,  with  chemical 
reactions,  arc  studied  in  detail  using  the  numerical  so¬ 
lution  of  the  Reynolds  Averaged  Navier-Stokes  (RANS) 
equations,  including  the  equations  for  chemical  species. 

The  present  investigation  utilizes  the  solution  of 
the  RANS  equations  for  flows  with  chemical  processes 
adapted  for  the  R.am  Accelerator  configuration  at  the  US 
Army  Research  Laboratory  (ARL)  by  Nusca^-i^  and  are 
summarized  in  this  paper.  In  the  present  calculations 
we  use  the  same  mixture  of  CH4,  O2,  and  N2  selected  at 
the  ARL  for  the  Ram  Accelerator  firings.  The  chemical 
reaction  rate  constants  for  a  mixture,  used  in  the  cal¬ 
culations  for  the  ARL  Ram  Accelerator,  yielded  good 
results  when  compared  with  experimental  data.^^i  Al¬ 
though  the  combustion  conditions  in  the  case  of  the  Ex¬ 
ternal  Propulsion  Accelerator  are  different  from  those 
in  the  Ram  Accelerator  (i.e.  scale  of  the  projectile  and 
method  of  ignition)  it  has  been  assumed  in  this  paper 
that  by  using  the  same  mixture  with  the  same  kinetics 
model,  good  results  can  be  achieved  for  the  External 


Propulsion  Accelerator  as  well.  Actually,  more  sutiable 
mixtures  for  the  External  Propulsion  mode  of  operation 
are  available  and  should  be  used,  but  in  the  present  in¬ 
vestigation  we  chose  to  use  this  mixture  for  which  the 
chemical  reaction  model  was  validated  for  the  Ram  Ac¬ 
celerator.  This  paper  presents  calculations  of  the  effects 
of  step  height  and  Mach  number  for  several  projectile 
geometries,  on  the  combustion  process  and  thrust  gen¬ 
erated  on  the  projectile.  These  calculations  are  also  used 
to  study  the  flow  ahead  of  the  forward-facing  step,  in¬ 
cluding  the  formation  of  a  combustion-detonation  front. 

Reacting  Flow  Model 

Computational  fluid  dynamics  (CFD)  flow  simula¬ 
tions  for  the  ram  accelerator  projectile  were  performed 
at  U.S.  Army  Research  Laboratory  (ARL)  using  the 
Rockwell  Science  Center  USA-RG  (Unified  Solution 
Algorithm  Real  Gas)  code. ^2-14  been 

used  successfully  at  the  ARL  for  simulation  of  a  full- 
bore  120mm  ram  accelerator  projectile  by  Nusca^^^ 
This  CFD  code  solves  the  full,  3D,  unsteady  Reynolds- 
Averaged  Navier-Stokes  (RANS)  equations  including 
equations  for  chemical  kinetics  (finite-rate  and  equilib¬ 
rium).  These  partial  differential  equations  are  cast  in 
conservation  form  and  converted  to  algebradc  equations 
using  an  upwind  finite-volume  formulation.  Solution 
takes  place  on  a  mesh  of  nodes  distributed  in  a  zonal 
fashion  around  the  projectile  and  throughout  the  flow 
field  such  that  sharp  geometric  corners  and  other  details 
are  accurately  represented.  The  conservation  law  form 
of  the  equations  assures  that  the  end  states  of  regions 
of  discontinuity  (shocks,  detonations,  deflagrations)  are 
physicaJly  correct  even  when  smeared  over  a  few  compu¬ 
tational  cells.  The  Total  Variation  Diminishing  (TVD) 
technique  is  employed  to  discretize  inertia  terms  of  the 
conservation  equations,  while  the  viscous  terms  are  eval¬ 
uated  using  an  unbiased  stencil.  Flux  computations 
across  cell  boundaries  are  based  on  Roe’s  scheme  for 
hyperbolic  equations.^®  Spatial  accuracy  of  third-order 
can  be  maintained  in  regions  of  the  flow  field  with  con¬ 
tinuous  variation  while  slope  limiting,  used  near  large 
flow  gradients,  reduces  the  accuracy  locally  to  avoid  spu¬ 
rious  oscillations. 

The  RANS  equations  for  2D/axisymmetric  reacting 
flow  {N  species  mixture)  are  written  in  the  follow¬ 
ing  conservation  form^^  for  dependent  variables  of  en- 
ergy  (e),  density  (p),  axial  and  radial  momentum  flux 
(pu.pv),  and  species  flux  (per). 
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The  shear  stress  terms  are  given  by, 
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In  these  equations  (Ti  and  u,*  are  the  mass  fraction  and 
chemical  production  term  for  the  i-th  species.  The 
species  viscosity  {p)  and  thermal  conducivity  («)  are 
referenced  to  /io,  /c©  and  T©  using  Sutherland’s  law, 

Pi  _(  T  + 

Pc.  \ToJ  T  +  S, 
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where  T©  and  S  vary  with  species.^*  The  mixture  vis- 
cosity  and  thermal  conductivity  are  determined  using 
Wilke’s  law^^  denoting  /  as  /i  or  /c, 
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where  X,-  and  Mi  are  the  mole  fraction  and  molecular 
weight  of  the  i-th  species,  respectively.  Pick’s  law  is  used 
to  relate  the  mixture  diffusivity  to  the  mixture  viscosity 
through  the  Schmidt  number  Sc  =  pirnUpD),  The  spe¬ 
cific  heat,  enthalpy  and  Gibbs  free  energy  of  each  species 
(per  mass)  are  given  by  the  following  fourth  order  tem¬ 
perature  (T)  polynomial  curve  fits,^®  {^^Ji  is  the  heat 
of  formation,  and  Sff,-  is  the  specific  gas  constant  for  the 
i-th  species). 
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The  mixture  enthalpy,  total  energy  per  unit  volume,  and 
ratio  of  specific  heats  are  given  by, 
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A  co-volume  equation  of  state  is  used  to  relate  pressure 
(p)  to  temperature  (T). 
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where,  V  is  the  volume,  b  is  the  covolume  (see  Ref.  19), 
n  is  the  number  of  moles  of  gas,  and  S  is  the  universal 
gas  constant. 

The  system  of  Navier-Stokes  equations  (Eq.  1)  is 
valid. for  the  lammzir  flow  of  a  viscous,  Newtonian  fluid.* 
In  reality,  the  flow  will  remain  laminar  up  to  a  certain 
critical  value  of  the  Reynolds  number,  pULjfji,  where 
U  and  L  are  representative  values^  of  the  velocity  and 
length.  Above  this  critical  value  the  flow  becomes  tur¬ 
bulent  and  is  characterized  by  the  appearance  of  fluc¬ 
tuations  in  all  variables  />,  T,  etc.)  around  mean 
values.  These  fluctuations  are  statistical  in  nature  and 
cannot  be  described  in  a  deterministic  way.  The  cur¬ 
rent  model  uses  the  Reynolds-averaging  approximation 
along  with  algebraic  turbulence  models  that  relate  fluid 
viscosity  (/i)  to  other  variables.  See  Reference  10  for 
further  details. 

Hydrocarbon  mixtures  such  as  3CH4  -h  2O2  +  ION2, 
pressurized  to  50-100  atm,  are  commonly  used  for  ram 
accelerator  testing  at  ARL.  A  single  test  can  incorpo¬ 
rate  several  different  gas  mixtures  in  different  sections 
of  the  accelerator  tube,  separated  by  diaphragms.  Fuel 
rich  mixtures  arc  used  with  fuel  equivalence  ratios  as 
high  as  3.  For  these  conditions  a  good  understand¬ 
ing  of  CH4/O2  chemical  kinetics,  especially  for  P  >  10 
atm,  is  not  available.^^  Accurate  numerical  simulation 
of  hydrocarbon-based  reacting  flow  is  very  demanding 
in  terms  of  computational  resources  since  the  number 
of  intermediate  species  and  the  number  of  kinetic  steps 
are  prohibitively  large  (i.e.,  >  15).  For  high  pressure 
systems  the  computational  investment  may  not  be  jus¬ 
tified  due  to  uncertainty  in  measured  reaction  kinetics. 
Thus,  global  reaction  mechanisms  (neglecting  interme¬ 
diate  steps)  have  been  used, 

CH4  +  2O2  CO  +  2H2  +  I.5O2  (17) 

CO  -l-  .5O2  ^  CO2 

2H2  +  O2  ^  2H2O 

The  reaction  rate  is  defined  using  the  Law  of  Mass  Ac¬ 
tion  and  an  Arrhenius  expression  for  C,  the  specific  re¬ 
action  rate  constant. 

N 
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where  AT°  is  the  collision  frequency,  the  exponential 
term  is  the  Boltzmann  factor,  Eq  is  the  activation  en¬ 
ergy,  and  Ui  represents  the  stoichiometric  coefficient. 
Exponents  a  thru  /  are  chosen  to  fit  results  from  flame 
experiments  using  large  kinetic  mechanisms. 


An  Arrhenius  reaction  rate  may  not  be  correct  for  the 
high  pressure  flows,  and  much  of  the  necessary  reaction 
rate  data  for  hydrocarbon  fuels  has  not  been  measured 
above  10  atm.^°  In  addition,  the  3-step  kinetics  mecha¬ 
nism  described  ajbove  is  based  on  the  low  pressure  work 
of  Westbrook^^  where  the  exponents  determined  by 
matching  experimentally  measured  flame  speed.  For  the 
ARL  full-bore  ram  accelerator,  where  mixture  ignition 
is  largely  determined  by  shock-induced  (i.e.  reflected 
shocks  from  the  accelerator  tube  waJl)  and  viscous  heat¬ 
ing,  the  exponents  a  thru  /  were  altered  to  the  degree 
necessary  in  order  to  match  accelerator  tube  wall  pres¬ 
sures  measured  at  the  ARL  (See  References  8-11).  It 
was  found  that  ram  accelerator  performance  was  deter¬ 
mined  to  a  large  degree  by  the  reaction  induction  length 
(i.e.  ignition  delay  time)  which  is  a  function  of  scale 
and  mixture  (as  well  as  pressure).  For  the  preliminary 
CFD  computations  described  in  the  present  paper,  no 
attempt  has  been  made  to  alter  the  kinetic  steps  or  the 
reaction  rate  expression  for  the  current  sub-bore  projec¬ 
tile  with  ignition  via  a  step  on  the  projectile.  Experi¬ 
mental  confirmation  of  mixture  ignition  (i.e.  location, 
extent,  and  resulting  tube  wall  pressure  history)  is  re¬ 
quired  in  order  to  validate  the  kinetics  model  used  in 
the  CFD  code.  However,  the  code  has  been  extensively 
validated  for  high  speed  compressible  gas  dynamics.^^ 

An  additional  issue  that  must  be  resolved  is  tur¬ 
bulence  modeling.  No  attempt  has  been  made  in  the 
present  work  (or  that  described  by  Nusca  for  the  ARL 
projectile)  to  test  the  dependence  of  results  on  vari¬ 
ous  models  of  turbulence.  Nusca^^  has  documented 
the  effects  of  viscosity  for  the  full-bore  projectile  and 
has  achieved  good  results  for  that  shape  using  algebraic 
turbulence  models.  No  attempt  has  been  made  in  the 
present  applications  to  alter  those  models. 

Results 

Figure  1  shows  the  computational  mesh  for  the  pro¬ 
jectile  with  a  step  located  at  the  forebody-midbody  junc¬ 
tion,  1mm  in  height.  The  forebody  is  conical,  85.1  mm 
in  length  with  cone  half-angle  10  degs.  The  midbody 
is  cylindrical,  30  mm  in  length  and  32  mm  in  diameter. 
The  afterbody  is  conical,  98.6  mm  in  length  and  with  a 
half-angle  of  3.5  degs.  The  grid  was  generated  in  zones 
so  that  sharp  body  junctions  are  accurately  represented. 
Approximately  11000  grid  cells  cover  the  projectile  and 
600  grid  cells  reside  in  the  wake  (for  the  half-plane). 

A  gas  mixture  of  3CH4-I-2O2+10N2  was  used  for  the 
CFD  simulations  described  in  this  section.  The  mixture 
has  a  frozen  speed  of  sound  a  =  361  m/s,  ratio  of  spe¬ 
cific  heats  7  =  1.376,  and  a  Chapman- Jouget  detonation 
velocity  Vcj  =  1450  m/s.  The  initial  pressure  and  tem- 
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perature  of  the  mixture  (before  projectile  injection)  was 
50  atm  and  300  K,  respectively. 

Figure  2  shows  the  computed,  frozen  flow  Mach  num¬ 
ber  contours  in  grey-scale  (bright  white  represents  zero 
value,  i.e.  the  projectile,  while  dark-  black  represents 
a'  value  of  6)  for  the  projectile  geometry  displayed  in 
Figure  1  (step  height  of  1  mm)  and  a  freestream  Mach 
number  of  6-(V"  =  2175  m/s,  V/Vcj  =  ^5).  Important 
gas  dynamic  features  are  labeled.  A- bow  shock  extends 
from  the  projectile  nose,  while  the  step  produces  a  nor¬ 
mal  shock  and  a  small  stagnation  region.  The  merging 
of  these  shocks  produces  a  curved  “transmitted  shock” 
that  reflects  from  the  accelerator  tube  w’all,  downstream 
of  the  projectile.  Behind  the  transmitted  shock,  a  con¬ 
tact  surface  is  formed  (not  fully  visible  in  Mach  number 
contours)  that  separates  flows  of  different  density,  tem¬ 
perature  and  velocity  (i.e.  flow  having  passed  thru  dif¬ 
ferent  shock  waves).  The  recirculation  zone  behind  the 
projectile  (wake)  is  clearly  subsonic.  The  results  for  a 
freestream  Mach  number  of  10  {V=  3610  m/s,  V/Vcj  = 
2.5),  displayed  in  Figure  3,  show  a  dowstream  bending  of 
the  bow  and  transmitted  shocks  as  well  as  a  strengthen¬ 
ing  (i.e.  producing  lower  veJues  of  the  dowmstream  Mach 
number  flow,  greyer  color)  of  the  transmitted  shock  and 
step-generated  normal  shock. 

Figures  4  and  5  show  water  (H2O  is  a  major  product 
of  the  combustion)  mass  fraction  contours  in  grey-scale 
(bright  white  representing  the  absence  of  water  and  dark 
black  representing  a  mass  fraction  of  0.07  and  above) 
for  a  freestream  Mach  number  of  6.  Figure  5  is  closer 
view  of  the  step  region.  Major  combustion  activity  is 
predicted,  using  the  current  kinetics  mechanism  (Eq.  17 
and  18),  in  the  step-generated  stagnation  region  (behind 
the  normal  shock)  and  continuing  around  the  step  and 
then  over  the  projectile  surface.  This  region  of  strong 
combustion  and  combustion  products  effectively  resides 
between  the  contact  surface  and  the  wall.  The  projectile 
wake  entrains  combustion  products.  Note  in  Figure  5 
that  the  production  of  H2O  is  reduced  somewhat  (i.e. 
lighter  grey  colors)  after  expansion  around  the  step,  but 
is  not  extinguished.  Figure  5  also  shows  a  close  view  of 
the  same  step  height  but  at  a  freestream  Mach  number  of 
10.  The  H2O  mass  fraction  distribution  if  more  uniform 
and  extends  farther  into  the  flow’field  radially  away  from 
the  projectile. 

A  close  examination  of  the  flow  structures  is  instruc¬ 
tive.  The  bow  shock  wave  from  the  nose  does  not  initi¬ 
ate  chemical  reaction.  The  detached  normal  shock  wave 
ahead  of  the  step  initiates  a  very  strong  detonation  w’ave 
which  intersects  the  bow  wave.  At  the  intersection  point 
a  transmitted  shock  wave  is  directed  towards  the  tube 
w’all.  The  strength  of  this  transmited  wave  depends  on 


the  flow  Mach  number  and  step  height.  The  transmited 
shock  weakens  as  its  distance  from  the  projectile  in¬ 
creases,  until  reaching  the  tube  wall  and  reflecting.  The 
pressure  behind  this  reflected  shock  wave  is  much  lower 
in  comparison  to  the  pressure  behind  the  reflected  shock 
waves  obtained  in  the  Ram  Accelerator.  At -the  inter¬ 
section  point  between  the  bow  shock  wave  and  the  det¬ 
onation  wave,  ahead  of  the  step,  a  contact-combustion 
surface  is  also  generated.  This  contact  surface  begins  at 
the  intersection  point,  which  is  a  few  step  heights  above 
the  surface,  and  is  propagated  from  there  nearly  paral¬ 
lel  to  the  projectile  surface  towards  the  projectile  base 
and  further  downstream,  above  the  wake.  The  temper¬ 
ature  of  the  flow  with  combustion  products,  behind  the 
step,  is  about  5  to  8  times  the  freestream  temperature, 
and  the  pressure  is  about  10  to  15  times  the  freestream 
pressure.  This  contact  surface  separates  the  hot  flow 
of  the  combustion  products,  processed  by  the  detona¬ 
tion  ahead  of  the  step,  and  the  outer  flow  that  passed 
through  the  transmited  shock  wave.  Therefore,  the  con¬ 
tact  surface  acts  as  an  aerodynamic  cowling  while  the 
region  between  the  contact  surface  and  the  projectile 
body  surface  acts  in  a  similar  way  as  the  combustion 
chamber  in  a  conventional  ramjet  engine.  However,  the 
“combustion  chamber”  is  confined  by  the  aerodynamic 
contact  surface  without  the  structural  cowling  of  a  con¬ 
ventional  engine.  The  high  pressure,  high  temperature 
combustion  products  in  this  aerodynamic  “combustion 
chamber”  can  then  be  expanded  over  the  projectile  after¬ 
body  to  produce  thrust.  This  is  the  first  demonstration, 
by  a  numerical  simulation,  of  an  aerodynamically  con¬ 
fined  combustion  region  external  to  the  projectile,  which 
can  be  utilized  for  thrust  generation. 

Figure  6  shows  water  mass  fraction  contours  in  the 
vicinity  of  the  0.5  and  1.5  mm  steps  for  a  freestream 
Mach  number  of  6.  The  H2O  mass  fraction  distributions 
are  similar  to  those  for  the  1  mm  step  height  but  with 
the  combustion  region  remaining  close  to  the  projectile 
body  for  the  0.5  mm  step  and  extending  farther  into  the 
flowfleld  for  the  1.5  mm  step.  The  larger  step  causes 
a  more  extensive  stagnation  region  in  front  of  the  step 
(behind  the  normal  shock). 

As  the  step  height  is  increased,  the  conditions  and  the 
position  of  the  contact  surface  are  effected,  resulting  in 
increased  mass  of  combustion  products  at  higher  pres¬ 
sures  and  temperatures,  but  with  an  increase  in  drag 
due  the  increased  step  height.  It  is  interesting  to  note 
that,  at  a  certain  Mach  number,  the  height  of  the  com¬ 
bustion  layer  above  the  projectile,  normalized  by  step 
height,  is  nearly  constant  for  the  three  step  heights  (0.5, 
1,  1.5  mm)  considered.  This  ratio  is  found  to  be  about 
5.5±0.2  at  Mach  6.  The  observation  that  combustion 
can  be  initiated  and  sustained  for  the  0.5mm  step  is  in 
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Figure  4:  Reacting  HoO  Contours  (bright  white  =  0,  dark  black  >  0.07),  A/qo  =  6,  1  mm  Step. 


Figure  5:  Reacting  HjO  Contours  (bright  white  =  0,  dark  black  >  0.07),  Mqo  =  6  and  10,  1  mm  St 


Figure  6:  Reacting  H2O  Contours  (bright  white  =  0.  dark  black  >  0.07),  =  6,  0.5  and  1.5  mm  Steps 


general  agreement  with  the  approximate  analysis  of  Ref¬ 
erence  4,  where  the  minimum  step  height  for  establish¬ 
ing  a  combustion-detonation  front  ahead  of  a  step  for 
stoichiometric  CH4,  O2  mixture  is  evaluated.  This  step 
height  is  about  0.7  mm  at  Mach  number  6  and  pres¬ 
sure  of  one  atmosphere.  The  agreement  between  the 
prediction  of  the  approximate  analysis^  and  the  present 
numerical  simulation  is  good.  * 

The  computed  pressure  distributions  on  the  projec¬ 
tile  surface  were  integrated  to  form  a  thrust  coefficient, 
F/PA.  where  P  is  the  initial  mixture  pressure  (50  atm) 
and  A  is  the  cross  sectional  area  of  the  accelerator  tube 
(diameter  of  120mm).  For  a  freest  ream  Mach  number 
of  6  and  the  1  mm  step,  the  non-reacting  flow  yielded 
F/PA  =  -0.4  (i.e.  drag),  while  the  reacting  flow  yielded 
F/PA  =  3.1.  Projectiles  with  step  heights  of  0.5  mm 
and  1.5  mm,  and  non-reacting  flow  yielded  F/PA  =  -0.3 
and  -0.9,  respectively.  With  chemical  reactions,  these 
cases  yielded  F/PA  =  3.4  and  3.06,  respectively.  The 
0.5  mm  step  yields  a  smaller  drag  that  must  be  cancelled 
by  the  thrust  produced  by  combustion  and  thus  yields  a 
larger  thrust  coefficient.  The  1.5  mm  step  yields  a  larger 
drag  that  must  be  cancelled  by  the  thrust  produced  by 
combustion.  Nevertheless  the  more  extensive  combus¬ 
tion  for  this  step  height  yields  a  slightly  smaller  thrust 
coefficient. 

The  authors  feel  confident  about  the  fidelity  of  the 
CFD  code  with  regards  to  gas  dynamic  phenomena  dis¬ 
played  in  this  section.  Grid  resolution  studies  should  be 
conducted  to  ensure  accurate  representation  of  the  wake 
region.  The  accuracy  of  the  assumed  chemical  kinetics 
mechanism  can  be  argued,  especially  at  high  pressure 
and  fuel  rich  conditions.  Further  research  and  validation 
is  warranted.  The  presented  results  should  therefore  be 
considered  preliminary. 

Conclusions 

The  present  numerical  simulation  of  the  effects  of  a 
forward-facing  step  positioned  on  the  body  of  a  hyper¬ 
velocity  projectile,  clearly  indicates  that  the  step  will 
produce  chemical  reactions.  The  present  calculations 
also  show  that  these  reactions  do  not  extinguish  after 
flow  expansion  around  the  step,  and  are  sustained  over 
the  afterbody  and  the  wake  of  the  projectile,  and  can 
be  used  to  produce  thrust.  However,  since  the  projectile 
shape  and  combustion  scsdes  are  different  in  the  present 
External  Propulsion  case  from  those  in  the  Ram  Accel¬ 
erator  case,  these  results  must  be  validated  by  an  ex¬ 
perimental  firing  of  a  projectile  with  a  step.  Such  an 
experimental  program  is  planned  by  the  third  author  at 
the  ARL,  and  should  be  reported  in  future  publications. 


It  is  shown  that  the  interaction  between  the  shock- 
wave  from  the  shallow*  nose  cone  and  the  detached  shock- 
detonation  wave  from  the  forward-facing  step,  which  is 
positioned  on  the  shoulder  of  the  projectile,  results  in 
confinement  of  a  combustion  layer  on  the  external  sur¬ 
face  of  the  projectile  afterbody  and  behind  the  blunt 
base.  The  high  pressure,  high  temperature  products 
of  the  combustion,  which  are  confined  to  the  external 
combustion  layer,  are  then  expanded  into  the  rear  and 
the  base  flow  of  the  projectile  and  can  be  used  to  gen¬ 
erate  thrust.  Thus,  this  layer,  confined  by  the  aero¬ 
dynamic  interactions,  acts  as  an  external  combustion 
chamber.  This  is  an  excellant  demonstration  of  the  ex¬ 
ternal  propulsion  method. 
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Stability  of  Hypersonic  Reacting  Stagnation  Flow  of  a 
Detoi^ble  Gas  Mixture  by  Dynamical  Systems  Analysis 

GENADI  nVANOV*  and  JOSEF  ROM^ 

Technion — Israel  InstituU  of  Technology,  Haifa  32000,  Israel 


The  stability  characteristics  of  the  reacting  hypersonic  flow  of  the  fucl/oaddizer  mixture  in  the  stagnation 
region  of  a  blunt  body  are  studied.  The  conditions  for  oscillations  of  the  combustion  front  are  assumed  to  be 
determined  mainly  by  the  flow  conditions  at  the  stagnation  region.  The  density  at  the  stagnation  region  is 
assumed  to  be  constant  at  hypersonic  flow  conditions.  By  assuming  a  simplified  flow  model,  the  time 
dependent  flow  equations,  including  the  heat  addition  due  to  the  chemical  reactions,  are  reduced  to  a 
second-order  nonlinear  differential  equation  for  the  instantaneous  temperature.  The  solutions  are  analyzed 
assuming  a  one-step  chemical  reaction  with  zero-order  and  first-order  processes  using  dynamical  systems 
methods.  These  methods  are  used  to  determine  the  stability  boundaries  in  terms  of  the  flow  and  chemical 
reason  parameters.  It  is  shown  that  the  zero-order  reaction  has  nonperiodic  solutions  that  may  lead  to 
explosion  whereas  the  first-order  and  higher-order  reactions  may  have  periodic  solutions  indicating 
oscillations.  The  zero-order  analysis  also  reaffirms  the  requirements  for  a  minimum  blunt  body  for  the 
establishment  of  a  detonation  Gn  agreement  with  classical  detonation  theory)  and  the  first-order  analysis 
indicates  a  minimum  body  size  for  establishment  of  oscillations.  The  oscOlation  frequencies  are  calculated 
using  the  small  perturbation  approximation  for  the  temperature  oscillations.  These  frequencies  are  compared 
with  results  from  published  data  on  spheres  and  hemisphere  c^indrical  bodies  fired  into  hydrogen'^oaq^gen 
and  acetylene  oxygen  mixtures.  Very  good  agreement  is  found  between  the  measured  and  calculated  results. 
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p  average  density 
E{T  -  r*) 

INTRODUCTION 

Under  certain  flow  conditions,  osciOations  of 
the  combustion  front  occur  on  blunt  bodies 
flred  into  detonable  mixtures  at  hypersonic 
speeds  [1,  2].  The  conditions  for  the  stabili^  of 
the  chemically  reacting  hypersonic  flow,  the 
establishment  of  a  steady  detonation  or  oscilla- 
toiy  reaction,  are  very  important  because  blunt 
bodies  are  used  for  the  generation  and  stabi¬ 
lization  of  reaction  fronts  in  supersonic  and 
hypersonic  propulsion  systems  such  as  the 
Oblique  Detonation  Wave  Engine  (ODWE) 
[3,  4]  and  the  External  Propulsion  Accelerator 
[5,  6].  Studies  of  the  development  of  combus¬ 
tion  and  detonation  fronts  on  blunt  bodies  in 
supersonic  and  hypersonic  flows  were  pre¬ 
sented  in  Refs.  7-11.  Although  the  oscillation 
of  the  combustion  front  is  discussed  in  these 
studies,  the  problem  is  not  analyzed  in  detail. 
A  wave  interaction  model  that  explains  the 
mechanism  for  the  establishment  of  oscilla¬ 
tions  near  the  stagnation  region  of  a  blunt 
body  was  presented  by  McVey  and  Toong  [12]. 
In  the  last  few  years  investigations  utilizing 
numerical  method  to  calculate  the  details  of 
the  time  dependent  flow  structure  near  the 
stagnation  region  of  a  blunt  body  at  hypersonic 
speeds  including  the  effects  of  chemical  reac¬ 
tions  were  presented  [13-15].  These  studies 
used  numerical  computations  to  map  the  in¬ 
stantaneous  wave  structure  near  the  stagnation 
region  of  a  blunt  body  and  to  infer  the  oscilla¬ 
tion  frequencies  for  these  conditions.  A  nu¬ 
merical  parametric  study  of  the  effects  of  vary¬ 
ing  body  diameter  is  presented  in  Ref.  16, 
showing  that  the  oscillations  are  not  only  a 
function  of  the  projectile  velocity  in  relation  to 
the  Chapman-Jouguet  detonation  velocity,  but 
also  depend  on  the  projectile  diameter  as  well 
as  on  the  induction  time,  reaction  rate  con¬ 
stant,  activation  energy  and  heat  release. 

It  can  be  seen  both  in  the  experimental  flow 
photographs  [1,  2]  as  well  as  in  the  results  of 
the  numerical  calculations  [13-16]  that  the  de¬ 


tached  shock  wave  remains  steady  and  is  only 
slightly  affected  by  the  oscillating  combustion 
front  in  the  stagnation  region.  Therefore,  the 
oscillating  combustion  is  bounded  on  one  side 
by  the  solid  body  surface  and  on  the  other  side 
by  the  steady  structure  of  the  detached  shock 
wave.  In  effect,  the  combustion  front  is  oscil¬ 
lating  like  a  membrane  in  a  region  bounded  by 
the  solid  body  surface  on  one  side  and  the 
stationary  detached  shock  wave  surface  on  the 
other  side.  This  results  in  a  “whistle”  effect  of 
the  vibrating  combustion  front.  This  oscillation 
model  suggests  the  application  of  the  dynami¬ 
cal  system  analysis  for  the  study  of  combustion 
oscillations  in  hypersonic  stagnation  flow  and 
provides  the  motivation  and  is  the  justification 
for  the  present  analysis. 

The  oscillations  are  observed  to  occur  at 
projectile  velocities  close  to  the  detonation  ve¬ 
locity.  When  the  projectile  velocity  is  increased 
considerably  above  the  detonation  speed  the 
combustion  front  merges  into  the  shock  to 
form  a  detonation  wave  and  no  oscillations  are 
expected. 

In  order  to  investigate  the  stability  of  the 
hypersonic  flow  of  the  fuel/oxidizer  mixture  in 
the  stagnation  region  of  a  blunt-nosed  body, 
the  time  dependent  flow  equations  including 
the  effects  of  chemical  reactions  are  studied. 
The  approximations  of  inviscid  flow  and  con¬ 
stant  density  in  the  stagnation  region  are  used. 
These  approximations  are  suitable  for  hyper¬ 
sonic  flows  at  high  Reynolds  numbers.  For 
example,  in  the  External  Propulsion  Accelera¬ 
tor  the  initial  pressures  are  in  the  range  of  20 
to  200  atm.  resulting  in  very  high  Reynolds 
numbers  and  very  thin  boundary  layers.  The 
fact  that  the  density  in  the  stagnation  region  of 
a  blunt  body  in  hypersonic  flow  is  roughly 
constant  is  discussed  by  Hayes  and  Probstein 
[17]  and  more  recently  in  Ref.  18. 

Studies  of  the  critical  conditions  in  chemi¬ 
cally  reacting  systems  by  the  examination  of 
the  equations  for  the  instantaneous  tempera¬ 
ture  (from  the  energy  equation)  and  the  species 
consumption  rate  (from  the  chemical  reaction 
equation)  are  presented  in  Refs.  19-22.  It  is 
shown  that  a  nonlinear  second  order  differen¬ 
tial  equation  describes  the  instantaneous  tem¬ 
perature  field.  This  nonlinear  differential 
equation  can  be  studied  using  the  theory  of 
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nonlinear  oscillators  presented  by  Andronov, 
et  al.  [23]  and  the  developments  of  the  mathe¬ 
matical  theory  for  these  equations  is  presented 
in  Refs.  24-26.  This  dynamical  ^tems  analy¬ 
sis  can  be  used  to  determine  the  conditions  for 
stability  of  the  solutions  and  their  periodic 
behavior. 


gas  equation  p  =  pRT,  in  Eq.  1,  gives 
dT  (  dT  3T\ 


=  -(2- 

1  dT 

dT 

y)  u— 

-1-  u — 

\ 

dy 

.r  ( 

E  \ 

r'K' 

'  rt)' 

THE  EQUATIONS  FOR  THE 
STAGNATION  REGION  FLOW  WITH 
CHEMICAL  REACTIONS 

The  time-dependent  energy  equation  for  the 
flow  of  the  fuel/oxidizer  mixture  including  the 
heat  addition  due  to  the  chemical  reactions 
(neglecting  effects  of  heat  conduction  and  mass 
diffusion)  is 


^vP\  —  +  “T"  + 


dp  dp  IE 

u - h  u - 1-  OBIA,  exp - 

dx  dy  RT 


where  the  chemical  reaction  rate  is  expressed 
by  the  Arrhenius  rate  equation 


In  order  to  accentuate  the  dynamics  of  the 
reactive  flow  system,  we  assume  a  simplified 
model  for  the  flow  in  the  stagnation  region.  It 
is  assumed  that  near  the  stagnation  point  of 
the  blimt  body  the  temperature  and  species 
density  vary  linearly,  so  the  temperature  gradi¬ 
ents  and  the  species  density  gradients  in  the  x 
and  the  y  directions  can  be  expressed  by  the 
approximate  linear  relations.  For  these  rela¬ 
tions  we  shall  assume  an  effective  temperature, 
T*,  and  an  effective  mass  fraction  of  the  main 
reaction  component,  /?*.  In  this  case  these 
gradients  are: 


dT 

T-T* 

dT 

T-T* 

dx  ~ 

R,  ’ 

H 

= 

B-B* 

B-B* 

dx  ~ 

R.  ' 

- 

dy  ~ 

H 

P”A„exp\  - 


The  initial  and  boundary  conditions  on  the 
temperatures  are:  at  r  =  0,  r(0)  =  Tq  and  on 
the  downstream  side  of  the  shock  wave  T 
(shock  wave)  =  T,  while  the  projectile  body 
temperature  is  T^.  Assuming  a  one-step  reac¬ 
tion,  the  main  component  of  reaction  varies 
during  the  reaaion  from  the  initial  value  /3  = 
Po  =  0)  to  the  final  value  p  =  Usually 
Bo  =  1  and  for  complete  reaction  Bead  0. 

For  hypersonic  flow  it  is  justified  to  assume, 
as  a  first  approximation,  that  the  density  in  the 
stagnation  region  behind  the  detached  shock 
wave  is  constant  [17,  18],  represented  by  an 
average  density  p-  In  this  case,  using  the  ideal 


Introducing  these  values  into  the  equations 
for  the  temperature  and  the  species  concentra¬ 
tions,  Eqs.  2  and  3,  gives 


dT  QA„r 

—  = - =— exp 

dt  c^p 


,  /  u  u\ 

-(.2- y)  —  -  —  (T-T*)  (5a) 

,  ri  I 


dB  /  E 

—  =  —A.  B"  exp - 


The  standoff  distance  of  the  detached  shock 
wave,  can  be  estimated  from  the  solution 
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for  the  hypersonic  flow  over  a  blunt  flat-nosed 
body  [17], 


R,= 


H 


1  -b 


where  b  =  — 
2 


(6) 


4  e  (  i, 

*”3;  2  ri;  ' 


and  ^  =  — . 

Pt 

From  Ref.  17  the  velocity  components  u  and  v 
near  the  stagnation  point  are 


u  = 


'  4-t) 


(8) 


and 


l-3e 


sinh 


eR, 


(9) 


--  - 

The  quantity;^^:,  ;;is  the  x  coordinate  of  the 
point  being  invesfigated,  and  the  quantity  x  is 
the  coordinate  at  entry  into  the  shock  layer  of 
the  streamline  passing  through  the  point  being 
investigated  [17],  Near  the  stagnation  point,  let 
y  =  KqH  and  x^  =  where  the  parame¬ 
ters  kq  and  /<i  range  between  0  and  1.  Using 
these  expressions  for  u,  v,  x,  and  y  (Eqs.  6-9), 
it  is  now  possible  to  evaluate 


u 

X 

where 


u 

H 


U.F 
H  ’ 


-  (1  -  w/T 


/Cl  + 


,  8£ 

~  T 


+  e 


3e 


1  -3e 


sinh 


/cqVI  -  3/ 


c(]  -b) 


(10) 


Introducing  the  following  dimensionless 
variables: 


18=-. 


B  = 


P  ^ 


a  = 


(2  -  y)Fu^t* 
~~H 


ixcJT*  exp— 
r*  = - ^ 


Q^nP 


X"-l  ’ 


and  the  Damkohler  number, 
t*u.F  a 


Dam  = 


H  2-y’ 


and  using  these  variables  in  Eqs.  5a  and  5b 
gives 

de  (  e  \ 

^  I  -  (2  -  r)Dam  6, 


1  +  jx6 


e 


—  =  -Bp''  exp- 
dr  ^1  +  /j.d 


(lla) 

-  Dam( p  -  p*) 

(llb) 


The  phase  plane  equation  is 

-(2  -  7)Dam0 


de 

dp 


1  +fj.e 


—Bp’'  exp 


6 


1  +/10 


-Dam(^  -  p*) 

(12) 

Equations  11a  and  11b  can  be  combined  to 
eliminate  p  and  provide  a  single  equation 
for  6.  At  the  stagnation  region  with  uniform 
reactive  mixture,  we  can  assume  that  ep/dx 
and  ep/  dy  c  1,  therefore,  the  value  of 
( p  —  p*)  and  then  the  combined  dy¬ 
namic  reaction  equation  is 


e+  ad - 


6(.6  -I-  0:0) 
(1  +  fidf 


+  nB(6  +  cxd)^  exp 


e 


1  +  1x6 


i/fi 


=  0, 


(13) 
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where  the  dot  represents  differentiation  with 
respect  to  t  and  the  n  represents  the  order  of 
the  chemical  reaction. 

The  dynamic  reaction  equation,  Eq.  D,  is  a 
second-order  nonlinear  differential  equation  of 
the  form 

e+f{e,e)e^g{6)  =  Q.  (i4) 

The  terms  fid,  6)  and  gid)  can  be  viewed 
as  representing  “damping”  and  a  “restoring 
force”  effects,  respectively.  It  should  be  noted 
that  these  terms  are  largely  determined  by  the 
body  size  parameter  a  =  (2  —  y)Fuj*/H 
(which  is  proportional  to  the  Damkohler  num¬ 
ber)  and  by  the  mixture  chemical  and  physical 
properties  represented  by  B  and  /x.  This  anal¬ 
ogy  is  demonstrated  by  the  “energy”  formula¬ 
tion  presented  in  Appendix  A. 

This  formulation  of  the  dynamic  reaction 
equation  developed  here  for  the  hypersonic 
stagnation  Sow  is  analogous  to  the  mechanical 
and  electrical  vibration  equations.  Therefore, 
the  large  mathematical  literature  and  the 
methods  developed  for  the  analysis  of  Dynami¬ 
cal  systems  can  be  applied  to  this  reacting 
system.  The  solution  of  this  equation  and  the 
singularities  which  appear  in  these  solutions 
can  be  used  to  define  the  stability  boundaries 
of  the  reactions,  the  possibility  of  divergence 
leading  to  explosion  or  existence  of  periodic 
solutions  indicating  oscillations.  In  this  analysis 
we  shall  use  the  solutions  of  similar  differen¬ 
tial  equations  presented  by  Andronov  et  al. 
[23],  and  in  Refs.  24-26  and  30-34. 

Some  indications  of  the  properties  of  these 
equations  can  be  obtained  from  numerical  so¬ 
lutions.  The  calculation  of  the  temperature 
parameter,  d,  variation  as  a  function  of  time 
parameter,  t,  for  the  case  of  a  =  0.1  and  for 
chemical  reactions  of  the  order  n  =  0, 1,  2  and 
autocatalytic  are  presented  in  Fig.  1.  The  case 
of  zero-order  reaction  shows  fast  increasing 
temperature  leading  to  explosion.  The  effect  of 
increasing  values  of  a  on  the  temperature  rise 
is  shown  in  Fig.  2.  For  small  a  the  tempera¬ 
ture  rises  fast  leading  to  ejq)losion.  As  a  in¬ 
creases  the  temperature  rise  decreases.  Since 
this  parameter  a  is  inversely  proportional  to 
the  body  size  H,  it  is  seen  that  above  a  certain 
body  size  the  temperature  of  reaction  diverges 


Fig.  1.  The  variation  of  9  as  a  function  of  r  for  various 
chemical  reactions,  with  n  «  0, 1,  2  and  autocatalytic. 


and  for  small  body  size  the  temperature  size  is 
limited.  This  indicates  that  there  is  a  minimum 
body  size  required  to  generate  explosive  reac¬ 
tion,  which  is  in  agreement  with  both  experi¬ 
mental  observations  [27]  and  analysis  [18]. 

Solutions  for  the  Zero-Order  Chemical 
Reaction,  n  =  0 

In  the  case  of  zero-order  reaction  Eqs.  11a  and 
11b  are 

dd  (  d  \ 

(15b) 

For  this  case  there  is  decoupling  between  the 
equations,  so  that  the  temperature  variation 
can  be  obtained  from  Eq.  15a  as  shown  in  Figs. 


Fig.  2.  The  variation  of  6  as  a  function  of  r  for  zero-order 
reaction  for  various  values  of  a. 
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1  and  2.  It  is  of  interest  to  examine  the  station¬ 
ary  states  in  the  B-a  plane,  using  the  relation 

e 

exp*: - :: 

1  -f-  fi6  ,  .  .  .  _ 

a  - - ,  obtained  from  Eq.  15a.  The 

& 

diagram  of  the  stationaiy  states  is  shown  in 
Fig.  3  for  cases  where  n  <  1/4.  There  are  two 
extreme  points,  A  and  B,  corresponding  to 
values  of  (6*,  af )  for  A  and  (02 » a?)  for  B,  as 
noted  on  the  graph,  where 


K2  = 


(1  —  2/i)  ±  \/l  -  4/1. 

V 


a 


1,2 


exp 


K2  ] 
I  +  Mer.2/ 
K2 


and 


The  low-temperature  reaction  following  the 
curve  up  to  point  A,  where  0  <  0f ,  represents 
the  slow  low  temperature  combustion  process. 
The  region  of  operation  between  A  and  B, 
where  0f  <  0  <  02 ,  represents  a  metastable 
state  and  the  high  temperature  region  of  the 
curve  above  point  B,  where  0  >  02 ,  may  be 
interpreted  as  an  eiqilosive  process.  In  fact, 


this  is  a  bifurcation  0(a)  where  af,  aj  are 
the  bifurcation  values  when  the  “jump”  condi¬ 
tions  are  realized.  This  is  a  duck  trajeaoiy 
with  unstable  state  within  0f,  02,  and  stable 
states  when  a  <  a^;  0  <  0*,  and  when  a  > 
a*;  8  >  02,  where  slow  combustion  or  explo¬ 
sive  regimes,  respectively,  exist 
These  extreme  points  can  be  used  to  define 
the  minimum  body  size.  Since  a  is  inversely 
proportional  to  the  body  radius  the  maximum 
value  of  aj  defines  a  minimum  value  of 
Therefore,  to  obtain  stable  detonation  we  must 
have 


(2  -  y)Fu^t* 


«2 


(16) 


This  result  is  in  good  agreement  with  the 
minimum  body  size  estimated  in  Ref.  18. 


Solutions  of  the  First-Order  Reaction 
Equations,  n  =  1 

Using  the  formulation  of  the  two  first-order 
equations  for  the  time-dependent  temperature 
and  species  concentration  parameters,  Eqs.  11a 
and  11b,  these  become  for  the  n  =  1  case: 

de  I  e  \ 

(17a) 

dp  I  S  \ 

"  Dam(^  - 

(17b) 


In  the  stagnation  region,  for  the  case  of 
uniform  combustion,  we  can  assume  Dam(  — 
P*)  =  0.  Then,  the  diagram  of  the  null-cline 
curves  in  the  phase  plane  0  —  jS  is  shown  in 
Fig.  4  for  various  values  of  a.  These  curves 
have  the  duck  shape  with  two  extreme  points, 
and  an  example  of  these  points  on  one  curve  at 
(0**,  p**)  and  at  (0**,  P**)  is  indicated  on 
Fig.  4.  These  curves  indicate  three  regions  of 
stability,  0  lower  than  0f*7^region  of  stable 
low  temperature;  0  between  0**  and  02 ’'-A-a 
metastable  state  and  for  0  above  0* *7texplo- 
sive  regime. 

The  shape  of  the  null-cline  curve  is  used  as  a 
guide  for  the  construction  of  the  phase  dia- 
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gram  for  a  range  of  values  for  a,  as  shown  in 
Fig.  5.  Starting  from  the  point  5  =  0;  ^  =  1,  in 
the  cases  of  large  a  (small  If)  the  solutions 
follow  the  stable  low  temperature  reaction  re¬ 
gion.  For  moderate  values  of  o,  the  solutions 
follow  the  metastable  region  (a^  in  Fig.  5).  For 
small  a  ( ttj  in  Fig.  5),  the  temperature  jumps 
to  the  explosive  reaction  region. 

The  second-order  equation,  Eq.  14,  is  used 
for  the  investigation  of  the  possibility  of  peri¬ 
odic  solutions.  In  the  case  of  zero-order  reac¬ 
tion,  n  =  0  and  g(0)  =  0,  there  can  be  no 
periodic  solutions.  For  the  first-order  reaction 
case,  n  =  1  in  Eq.  14,  the  expressions  for 
/(d,  6)  and  g(d)  are,  respectively, 


f(e,e)  =  a 


ad  -  6 
(1  + 


+  B  exp 


6 


1  +  1x6 
(18a) 


g(e)  =Ba0exp|j-j^j.  (185) 

Since  6g(0)  >  0  for  ^  >  0,  the  point  [6{t  =  0), 
=  0)]  is  a  singular  point.  TTiis  singular  point 


is  inside  the  closed  trajectories  of  the  energy 
curves,  as  shown  in  Ref.  25,  where  the  energy 
curves  are 


B)  =  GiB)  +  —6^  =  const, 
where 

rg(B)de^G(B),  (19) 

■'0 

Since  g(-B)i^  ^g(0X  these  paths  are  not 
sjTnmetrical  trajectories  near  the  initial  point 
[d{t  =  0),  Bit  =  0)].  The  conditions  defined  by 
Eq.  19  are  sufficient  to  obtain  oscillations  of  B 
and  B  or  to  obtain  monotonic  decay  of  the 
function  Bit\  so  that  lim^(r),^o  = 
lim  =  0.  It  is  expected  that  the  peri¬ 

odic  solutions  of  Eq.  14,  for  /z  =  1,  will  be 
cycles  located  around  the  origin  of  coordinates. 
The  existence  of  a  unique  solution  of  thk 
equation  and  its  periodic  behavior  are  dis¬ 
cussed  in  Refs.  25,  26,  and  30-34.  Since  these 
analyses  present  general  theorems,  their  appli¬ 
cation  to  the  present  case  requires  specific 
considerations  which  are  presented  in  Ap- 
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pendix  B.  When  a  =  0,  the  trajectories  are 
open,  indicating  nonperiodic  solution  and  pos¬ 
sible  divergence  to  explosion  as  in  the  zero- 
order  case.  As  a  increases  the  trajectories 
become  closed  indicating  a  limit  cycle  with 
decreasing  amplitude,  as  shown  in  Fig.  6  for 
various  initial  conditions.  As  discussed  in 
Ref.  24,  the  existence  of  a  limit  cycle  indicates 
that  periodic  solutions  are  possible.  These  pe¬ 
riodic  solutions  can  be  studied  from  the  char¬ 
acteristics  of  the  singular  points  of  the  system. 
The  oscillatory  solutions  and  their  frequencies 
can  be  also  studied  using  perturbation  analysis. 

EVALUATION  OF  OSCILLATION 
FREQUENCIES  USING  PERTURBATION 
ANALYSIS  FOR  THE  FIRST-ORDER 
REACTION,  n  =  1 

We  apply  small  perturbation  analysis  to  the 
instantaneous  temperature  and  reacting  species 
equations  fc>r  the  first-order  reactions,  Eqs.  17a 
and  17b,  respectively,  by  assuming  small  distur¬ 
bances  imposed  on  the  average  flow  parame¬ 
ters,  6  =  6+  a',  p  =  p  +  p’. 

In  the  stagnation  flow,  for  first-order  reac¬ 
tion,  n  =  1,  and  Dam( p  —  p*)  =  0,  0  and  p 
are  evaluated  using  Eqs.  17a  and  17b  to  obtain 

-  (  ^  \ 

P  ■; - =  -  (2  —  ■y)Dam  6  =  0,  (20a) 

\l  +  IJ.6) 


equations  are  obtained: 


\  1 

1  +  fid 


-  a  \6’ 


Bp  exp 


1  + 


This  linearization  requires  P6'  ,  fi6'  -s:  1, 

and  'p 

(  6’  ]  6' 

*^\1+M®/  l  +  fid 

The  stability  of  the  solutions  of  these  equa¬ 
tions  and  their  oscillation  frequencies  can  be 
studied  by  the  solution  of  the  characteristic 
equation  for  the  eigenvalues, 

A^-l-  Ib - ^lexpf — A 

1  1  -f  /i6  j  i  1  / 


Ba  exp  - =  =  0. 

I  1  +  Mel 


By  introducing  the  perturbation  values  into 
Eqs.  17a  and  17b,  the  following  linearized 


I- 


Fig.  6,  Phase  porffait  for  first-order  reaction  for  various 
initial  conditions  and  various  a. 


It  is  well  known  that  the  stability  and  fre¬ 
quencies  of  oscillations  are  determined  by  the 
complex  roots  of  this  characteristics  equation. 
For  the  hydrogen-oxygen  mixture,  flowing  at 
5000  m/s,  where  =  0.2,  B  =  0.5,  y  =  1.4, 
t*  =  3.5  X  10"’,  and  using  Kq,  Kj  =  0.5, 
the  eigenvalues  for  this  case,  for  a  — 
(0.01,0.2,1,2),  are  plotted  in  Fig.  7.  It  is  seen 
that  small  a  <  1,  there  are  positive  real  parts 
indicating  regimes  of  diverging  temperatures. 
Here  we  have  the  possibility  of  explosion.  The 
range  of  values  of  the  imaginary  parts  of  the 
eigenvalues  for  this  mixture  indicate  that 
the  frequencies  of  the  oscillations,  /  = 
Im(A)/27rr*,  are  between  0.05  and  5.00  MHz, 
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n»(  X) 


Fig.  7.  Variation  of  the  real  and  imaginaiy  parts  of  the 
eigenvalues  for  a  hydrogen-oxygen  mixture. 


Fig.  8.  Neutral  stability  curves,  Re(A)  -  0,  for 
hydrogen-oxygen  mixture. 


which  agrees  well  with  frequencies  measured 
in  experiments. 

The  complex  eigenvalues  vary  as  a  function 
of  the  temperature  parameter  for  a  known 
mixture.  The  oscillations  occur  when  the  imag¬ 
inary  part  is  nonzero.  The  range  where  the 
imaginaiy  part  is  not  zero,  Le.,  the  range  of 
oscillations,  is  obtained  from  the  solution  of 


5  - 


exp 


e 

l  +  fxO  j 


|2 


+  a 


I 

—  4Ba  exp 

I 


1  +  fid  j 


^0. 


(23) 


The  boundary  within  which  oscillation  oc- 
ctirs  is  obtained  when  5  =  0.  The  neutral  sta¬ 
bility  boundaries  between  stable  and  unstable 
regimes,  defined  by  the  curves  of  Re(A)  =  0, 
are  calculated  for  this  hydrogen-oxygen 
mixture,  as  shown  in  Fig.  8,  for  various  values 
of  a.  The  variation  of  the  firequency  of  oscilla¬ 
tions,  as  determined  by  the  values  of  the  imagi¬ 
nary  part  of  the  eigenvalues,  as  a  function  of  6 


is  presented  in  Fig.  9.  For  small  values  of  a 
the  frequencies  are  low,  although  as  a  in¬ 
creases  the  values  of  Im(A)  increase.  There 
seems  to  be  a  peak  of  frequency  for  certain 
values  of  6,  indicating  an  effect  like  “reso¬ 
nance.”  For  increasing  a  the  peak  in  the  val¬ 
ues  of  Im(A)  Jncreases  and  as  is  increased 
the  value  of  6  at  the  peak  of  the  Im(A)  also 
shifts  to  higher  values  of  6.  The  variation  of 
this  peak  value  of  Im(A)  at  increasing  a  and  6 
is  also  indicated  on  Fig.  9. 


DISCUSSION  OF  THE  DYNAMICAL 
SYSTEMS  ANALYSIS  FOR  THE 
PERTURBATION  EQUATIONS 

The  perturbation  equations  presented  in  Eqs. 
21a  and  21b,  can  be  expressed  in  the  “classic” 
form 


e'  =ae’  +  bP', 

13'  =  -cG’  -dl3', 
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where  a,  b,  c,  and  d  are 


a  = 


/3exp 


6 


1  +  /xe, 


b  =  exp 


Bp  exp 


1  +  IJ.6 

— i 

1  +  fjid  I 

e 


c  = 


1  +  ^jl6 


1  +  1x6 


—  a 


d  = 


B  exp 


6  ' 
1  +  1x6 1 


Following  the  analysis  of  Ref.  23  we  define 


A  =  —ad  +  be 


as  the  trace  and  the  determinant  of  the  matrix 
of  coefficients. 

Then 


cr  = 


1  +  fx6 


■=  -B 


exp 


a 


■=]- a  (24) 
1  +  A®; 


A  =  aB  expl 


6 


1  +  t^S, 


(25) 


In  Fig.  10,  the  well  known  diagram  of  a-  vs. 
A  is  presented,  showing  the  regions  of  stability 
and  instability  as  well  as  the  ^e  of  singulari¬ 
ties  in  each  region.  The  stationary  states  are 
evaluated  by  Eqs.  20a  and  20b  and  the  nature 
of  singularities  near  these  stationary  states  can 
be  evaluated  by  analyzing  the  perturbation 
equations  using  the  parameters  a-  and  A  as 
defined  in  Eqs.  24  and  25,  respectively.  It  was 
already  shown  that  for  the  first-order  reaction 
case,  n  =  1,  there  js  o^y  one  singular  point  at 
the  origin,  where  6  =  p  =  0  and  the  eigenval¬ 
ues  are  Aj  =  -a  and  X2  =  indicating  a 
stable  node  type  singularity.  This  result  can  be 
shown  to  apply  also  to  higher  order  reaction 
cases  where  n  ^  1. 

We  examine  the  condition  cr  =  0.  The  real 
part  of  the  eigenvalue,  A,  =  0,  but  the  imagi¬ 
nary  part  A;  0,  indicates  neutral  oscillations 
corresponding  to  a  singularity  of  a  center  type. 
The  corresponding  body  size  which  is  evalu¬ 
ated  from  the  corresponding  value  of  a  can  be 


Fig.  10.  Relation  between  states  of  equilibrium  and  types 
of  singularities  of  the  characteristics  equation. 


a-  =  a  —  d, 
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calculated  from  Eq.  24,  when  cr  =  0,  so  that 
this  critical  value  of  is 

(2-r)f«.a-M«)«p(Y^) 

(26) 

This  is  the  body  radius  required  for  obtain¬ 
ing  neutral  oscillations  of  the  reactions  in  the 
stagnation  point  value  for  first-order  reactions. 
This  critic^  body  size  is  different  from  the 
minimum  body  radius  required  for  establish¬ 
ment  of  detonation  evaluated  from  the  appro¬ 
priate  criterion  for  the  zero-order  reaction, 
presented  in  Eq.  16. 


COMPARISON  WITH  EXPERIMENTAL 
DATA 

Minimum  Body  Size 

The  analysis  for  the  zero-order  reaction  case 
resulted  in  the  evaluation  of  a  minimum  body 
radius  for  establishment  of  detonation  in  the 
stagnation  point,  presented  in  Eq.  16.  The  min¬ 
imum  radius  can  be  calculated  using  the  proper 
parameters  for  various  mixtures,  such  as  hy¬ 
drogen-oxygen,  methane-oxygen  and  acety¬ 
lene-oxygen.  For  these  mixtures  the  calculated 
body  radius  is  about  1-3  mm  at  flow  Mach 
number  of  5  down  to  0.1-03  mm  at  Mach 
number  of  10.  These  results  are  in  good  agree¬ 
ment  with  the  previous  estimates  of  Ref.  18. 
Experimental  evidence  showing  the  require¬ 
ment  for  a  minimum  size  cylindrical  rod  for 
establishment  of  detonation  is  presented  in 
Ref.  27.  The  experiments  are  carried  out  in  an 
expansion  tube  filled  with  methane— oxygen- 
nitrogen  mixture  accelerated  to  velocity  of 
about  2300  m/s  at  initial  pressure  of  03  atm 
For  these  conditions  detonation  was  not 
achieved  for  rod  diameters  lower  than  3  mm 
This  value  is  within  the  range  of  estimated 
values  by  the  present  analysis,  Eq.  16.  The 
minimum  body  size,  obtained  here  by  the  dy¬ 


namical  tystems  analysis,  is  in  agreement  with 
previous  results  obtained  using  classical  homo¬ 
geneous  ignition  theory  [35-37]. 

An  indication  that  there  is  a  minimum  body 
diameter  for  oscillations  to  occur,  as  indicated 
by  the  critical  body  size  in  Eq.  26,  can  be  seen 
from  the  numeric^  study  on  the  effect  of  body 
size  variation  presented  in  Ref.  16.  It  was 
shown  there  that  whereas  oscillations  were  ob¬ 
served  on  a  15  mm  sphere,  such  oscillations 
were  not  obtained  for  Ae  same  flow  conditions 
on  a  23  mm  sphere. 

Frequencies  of  Oscillations 

There  are  few  ejqperimental  measurements  of 
the  oscillations  obtained  on  spheres  and  spher¬ 
ical  nosed  cylindrical  projectiles  flying  at  high 
supersonic  speeds  in  hydrogen/oxygen  diluted 
with  nitrogen  or  with  argon  and  acetylene/ 
oxygen  mixtures  [1,  2,  12,  28,  29].  The  test 
conditions  and  the  measured  frequencies  re¬ 
ported  in  these  dted  references  are  presented 
in  Table  1. 

It  was  shown  here  that  for  reactions  of  order 
one  and  higher,  oscillations  can  be  obtained 
when  the  phase  trajectories  are  closed  indicat¬ 
ing  a  limit  cycle.  The  period  of  oscillations  is 
estimated  from  the  relation  for  the  limit  cycle, 

exp 


Introducing  the  values  for  the  hydrogen- 
oxygen  system,  the  estimated  frequencies  by 
this  relation  are  in  the  range  of  0.1  to  13  MHz. 
This  range  of  estimated  frequencies  covers  very 
well  the  experimental  values  presented  in  Table 
1. 

The  oscillation  frequencies  are  evaluated 
from  the  perturbation  analysis  by  calculating 
the  values  of  the  imaginary  part  of  the  eigen¬ 
values  for  the  first-order  reaction,  n  =  1,  using 
Eq.  22.  The  oscOlation  frequency  is  then. 


^  Im(  A) 


2c„^r*  exp|-i 


exp 


1  -f 


1  + 
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TABLE  1 

Comparison  of  Measured  and  Calculated  Frequencies  of  Oscillations  on  Blunt  Bodies  Flying  at  Hypersonic  Speeds 

in  Detonative  Mixtures 


/  / 


Mixture 

M. 

4. 

(m/s) 

(m/s) 

Pm 

(atm) 

D 

(m/s) 

d-/D 

H 

(mm) 

(MHz) 

Exp. 

(MHz) 

Caic. 

Ref. 

2H2  +  Oj  +  3.76N2 

4.18 

403 

1685 

0.421 

2055 

a82 

7.5 

0.15 

0.12 

2 

2H2  +  O2  +  3.76Nj 

5.04 

403 

2029 

0.421 

2055 

a99 

7.5 

1.04 

1.15 

2 

H2  +  Air 

4,65 

409 

1900 

0.55 

1950 

a97 

4JS 

0.85 

1  Q 

0.9 

28 

H2  +  Air 

4.71 

409 

1925 

0.55 

1950 

0.99 

4.5 

1.0 

1.15 

130 

28 

H2  +  Air 

4.16 

409 

1700 

0.55 

1950 

a87 

4.5 

03 

033 

28 

H2  +  Air(30%H2) 

4.8 

409 

1963 

0.25 

1950 

1.01 

10 

*0.1 

ai2 

29 

O2  +  C2H2(15:1) 

5.16 

328 

1690 

0.197 

1692 

1.00 

6J5 

*0.06 

0.07 

12 

H2  +  02  +  Ar  (1:2:2) 

4.54 

370 

1680 

a263 

1910 

ass 

635 

*036 

0.29 

1 

These  calculated  frequencies  are  also  in¬ 
cluded  in  Table  1,  for  comparison  with  the 
experimental  data.  The  values  of  the  parame¬ 
ters  used  in  these  calculations  are  presented  in 
Table  2.  In  all  cases  we  chose  the  values  of 
kq  =  0.5  and  =  0.75,  as  these  seem  to  fit 
somewhat  better  the  results,  particularly  those 
of  Ref.  28,  instead  of  kq  =  0.5  and  =  0.5 
which  were  chosen  previously.  There  is  a  re¬ 
markable  agreement  between  the  meastired 
and  calculated  frequencies  for  all  the  reported 
ejqjeriments,  the  differences  range  from  4.3% 
to  20%,  as  indicated  in  Table  2.  These  differ¬ 


ences  are  relatively  small  considering  the  rough 
approximations  used  in  the  analysis. 

DISCUSSION  AND  CONCLUSIONS 

The  present  analysis  shows  that  the  chemical 
reaction  process  in  the  stagnation  region  of  a 
blunt-nosed  body  in  hypersonic  flow  of  a  react¬ 
ing  mixture  behaves  as  a  dynamical  system 
analogous  to  the  vibrating  mechanical  and 
electrical  systems.  In  the  stagnation  region  the 
combustion  front  oscillates  between  the 
nonoscOlating  detached  shock  wave  and  the 


TABLE 2 


The  Physical  Parameters  Used  in  the  Calculations  of  the  Frequency  of  the  Oscillations  in  Table  1* 


Mixture 

Q 

(MJ/kg) 

r* 

(KD 

B 

E/R  X  10^ 

X  10®  X  p 

e 

/Exp 

(MH£I 

/Calc. 

(MHz) 

%  error 

2H2  +  O2  +  3.76N2 

32 

1967 

0.205 

0.55 

0.96 

1.25 

131 

0.15 

o.ii 

20 

2H2  +  O2  +  3.76N2 

32 

2200 

0.23 

0.69 

0.96 

1.25 

0.74 

1.04 

1.15 

10.6 

H2  +  Air 

1.98 

2085 

0.208 

0.62 

1.0 

32 

03 

0.85 

0.9 

6 

Hj  +  Air 

1.98 

1800 

0.18 

0.46 

1.0 

32 

134 

1.15 

1.20 

4.3 

H2  +  Air 

1.98 

1650 

0.165 

038 

1.0 

32 

339 

03 

033 

10 

+  Air  (30%  Hj) 

1,98 

1830 

0.183 

0.47 

1.0 

32 

1.4 

*  0.1 

0.12 

20 

O2  +  C2H2  (15:1) 

1.49 

2060 

0.165 

0.16 

1.25 

(2630  X  p)^^) 

0.41 

*0.06 

0.07 

17 

H2  +  02  +  Ar(l:2:2) 

1.90 

2050 

0.214 

0.81 

0.96 

1.25 

1.2 

*  0.26 

0.29 

11.5 

is  the  pressure  equal  to  the  pressure  behind  the  normal  detached  shock  wave,  p. 
r*  is  evaluated  as  follows:  T”  —  OJ  CTc-j  +  where  Tc-/  is  the  Qiapman-Jouguet  detonation  temperature  and  Tg 
is  the  Eckeret  effective  temperature,  -  OJ  {?„  +  T,)  +  0.22  (T,  -  T,)  where  T,  -  p./Rp,  and  T  -  +  0.5r(-y  - 

1)^^]  with  r  -  0.85. 

6  is  evaluated  by  the  relation 


£  y  +  1 

(J-T’),  where  T  -  Tc- , 


PT-i 


^  ^  for  the  Acytden^^^gen  mixture  the  value  is  (2.62  X  10*  X  p),  where  p  is  the  average  density  equal  to  the 


density  behind  thej 


normal  detached  shock  wave,  p. 
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solid  body  surface.  The  instantaneous  temper¬ 
ature  variation  in  the  hypersonic  flow  with 
chemical  reactions  which  may  lead  to  detona¬ 
tion  or  to  oscillatory  combustion  is  formulated 
as  a  second  order  differential  equation  with 
“damping”  and  “restoring  force”  terms.  These 
“damping”  and  “restoring  force”  terms  are 
dependent  on  the  flow  characteristics  and  on 
the  chemical  and  physical  parameters  of  the 
fuel/oxidizer  mixture.  This  formulation  is  de¬ 
rived  from  the  energy  equation  for  the  hyper¬ 
sonic  flow  at  the  stagnation  region  of  a  blunt- 
nosed  body  and  the  associated  species  concen¬ 
tration  equation  for  the  reacting  mixture.  In 
this  formulation  these  equations  can  be  pre¬ 
sented  either  as  two  first-order  nonlinear  dif¬ 
ferential  equations  for  the  temperature  and 
species  concentration  or  as  a  second-order 
nonlinear  differential  equation  for  the  instan¬ 
taneous  temperature. 

In  this  development  some  rough  approxima¬ 
tions  are  required,  as  presented  and  discussed 
in  the  paper.  These  include:  (1)  the  assumption 
of  “constant  density”  in  the  stagnation  region; 

(2)  the  assumption  of  linear  variation  of  the 
flow  parameters  near  the  stagnation  point 
which  enabled  the  evaluation  of  the  local 

and  d/dy  derivatives  near  the  stagnation  point 
in  terms  of  average  values  (Eqs.  4a  and  4b); 

(3)  representing  the  velocity  components  in  the 
stagnation  region  u,  v  by  Eqs.  7-9,  evaluated 
for  blunt  body  in  hypersonic  flow  [17];  (4) 
choosing  average  values  of  the  parameters 

=  0.5  and  kj  =  0.75  for  the  representative 
y/H  and  x-^/R,  positions  in  the  stagnation 
region,  respectively,  and  (5)  assuming  uniform 
reactive  mixture  in  the  stagnation  region,  so 
that  (/3-/3*)->’Oin  this  region. 

The  result  is  a  formulation  of  the  dynamical 
systems  model  for  the  combustion  front  which 
can  oscillate  behind  the  steady  detached  shock 
wave  and  ahead  of  the  blimt  nosed  body  in 
hypersonic  flow.  The  full  analogy  between  this 
combustion  process  and  mechanical  and  elec¬ 
trical  dynamic  systems  is  a  very  intriguing  sub¬ 
ject  for  studies.  It  is  shown  in  this  investigation 
that  it  is  possible  to  identify  “damping”  and 
“restoring  force”  in  terms  of  the  flow  and 
chemical  and  physical  parameters.  This  is  made 
even  clearer  using  the  energy  analogy  formula¬ 
tion  presented  in  Appendix  A. 


The  analysis  of  the  zero-order  reaction  indi¬ 
cates  the  existence  of  a  solution  leading  to 
explosion  when  the  body  is  large  enough.  This 
is  used  for  the  determination  of  the  minimum 
blunt  body  size  required  to  establish  combus¬ 
tion-detonation  process.  This  result,  obtained 
here  by  the  dynamical  systems  analysis,  is  in 
agreement  with  results  obtained  previously  us¬ 
ing  classical  homogeneous  ignition  theory 
[35-37]  and  with  experimental  observations 
[27].  It  is  also  shown  that  for  first-order  reac¬ 
tions  (and  higher  orders)  the  possibility  of  peri¬ 
odic  solutions  exists.  Using  small  perturbation 
anafysis,  a  characteristic  equation  for  thie 
eigenvalues  is  obtained  from  which  the  oscilla¬ 
tion  frequencies  are  evaluated.  The  calculated 
frequencies  are  compared  with  measured  fre¬ 
quencies  on  spheres  and  hemispherical  cylin¬ 
drical  bodies.  The  differences  vary  between  4% 
and  20%,  which  is  a  very  good  agreement  in 
this  case.  The  fact  that  such  good  agreement  is 
obtained,  in  spite  of  the  very  rough  approxima¬ 
tions  used,  is  an  indication  that  the  present 
formulation  of  the  combustion  process  using 
the  dynamical  systems  analysis  is  reasonable 
for  the  stability  analysis  of  hypersonic  flow  of 
reacting  mixtures  in  the  stagnation  region  of 
blunt-nosed  bodies. 
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APPENDIX  A:  ENERGY  ANALOGY 
FORMULATION 


An  additional  view  of  the  analogy  between  the 
dynamical  systems  in  the  hypersonic  reacting 
stagnation  flow  and  the  dynamical  systems  in 
mechanical  and  electrical  vibrations  can  be 
derived  by  examining  the  energy  formulation 
of  the  dynamic  equation.  The  energy  formula¬ 
tion  can  be  obtained  by  examining  £q.  13  for 
the  zero-order  reaction  case,  n  =  0,  by  first 
differentiating  this  equation  with  respect  to  r 
and  obtaining  the  relation 
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exp 
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1  +  fie 


(1  +  tier 


(Al) 


Using  the  relation  6=6 


ide 


and  integrating 


Eq.  Al  between  the  limits  of  Sq 
convenienc^can  assume  /u.  =  0  (since  our  in¬ 
terest  is  mainly  to  obtain  the  energy  terms), 
then 
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exp(20o) 
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This  relation  has  the  form  of  an  energy  repre¬ 
sentation. 


KCd)  +  V(a,e)=£(cx,eo,6o),  (A3) 

where  AT  =  6^/2  represents  the  “kinetic  en¬ 
ergy,” 

V  =  a0exp(0)  -  [exp(2e)/2]  -  a^^/2, 
represents  “potential  energy,”  and 

E  =  6g/2  a6g  exp(eo) 

-[exp(20o)/2]  -  crV/2 

=  const^represents  “total  energy.”  /\  'O  J' ) 
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Therefore,  Eq.  A3  has  the  form  of  the  en¬ 
ergy  conservation  of  a  dynamical  ^tem,  and 
can  be  used  as  a  phase  trajectory  equation 
obtained  from  “energy”  conservation  relations 
when  E  =  const  The  “total  energy”  constant 
E  is  determined  by  the  values  of  which  can 
be  viewed  as  determining  the  “energy  reserve” 
for  the  system  while  the  parameter  a  repre¬ 
sents  the  “stiffiness”  of  the  system.  This  me¬ 
chanical  analogy  shows  that  for  this  zero-order 
reaction  no  periodic  solutions  exist,  in  agree¬ 
ment  with  our  analysis. 

Examining  the  second-order  differential 
equation  for  the  instantaneous  temperature 
parameter  in  the  zero-order  reaction  case,  Al, 
shows  that  the  “damping,”  which  is  repre¬ 
sented  by  the  term  multiplying^  Eq.  Al,  is 
approximately  equal  to  [e^fl)  —  a]  assuming 
fj,  =  0  for  this  case.  So  when  a  is  increased, 
i.e.,  body  size  is  decreased,  then,  the  “damping” 
is  also  decreased. 

APPENDIX  B:  MATHEMATICAL 
CONDITIONS  FOR  SOLUTIONS  OF  THE 
DYNAMICAL  REACTION  EQUATION 

Conditions  for  Periodic  Solutions 

The  dynamical  reaction  equation,  Eq.  14,  is 

e+f(e,e)e  +  g(.e)  =o.  (bi) 

I 

•'  /Using  the  analysis  of  Eq.  34,  examining  the 
solutions  of  this  equation  for  the  following 
conditions; 

1.  f{do,  0o)  <  0,  where  6^  =  d(t  =  0) 

and  Bo  =  B(t  =  0)  (B2) 

2.  lg(6)l  +  lBlf(B,e)  2:s>0JorlBl>  Bg 

(B3) 

3.  f(6,B)  +f(B,-6)  >0 

for  |0|  +  |e|  ^  const.  (B4) 

In  these  cases,  according  to  Ref.  34,  the 
dynamical  reaction  equation,  Eq.  Bl,  has  at 
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This  requirement  is  still  very  restrictive  in 
comparison  to  those  presented  in  relations  B5b 
or  B6,  which  lead  to  the  conditions 

a<l-B  if  00  =  0  fromEq.BSb, 
a<l  M  B  E\  IX  e  from  Eq.  B6 

(where  f  is  an  infinitesimal  value)  while  from 
the  condition  of  Eq.  B8,  it  can  be  shown  that 


cannot  be  treated  with  the  present  one-step 
reaction  analysis. 

Conditions  on  the  “Damping”  and  “Restoring 
Force”  Parameters 

Examining  the  dynamic  reaction  equation,  Eq. 
Bl,  neutral  free  oscillations  occur  when 

/(0,  0)  =  0;  then 


a  <  1  4- 
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/I-2 


0"“ 
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-f-  B  exp 
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So  that  in  order  to  obtain  periodic  solutions,  it 
is  required  that  a  <  1  4-  In  such  cases  when 
a  is  kept  constant  and  for  various  initial  tem¬ 
peratures,  00,  the  solution  is  a  family  of  dosed 
trajectories  as  outlined  in  the  phase  portrait 
shown  in  Fig.  11.  Solutions  are  shown  when 
starting  from  the  stationary  state  0,  which  is  a 
stable  node  evaluated  from  Eq.  17a,  and  for 
various  other  initial  conditions.  An  example  of 
a  phase  portrait  for  various  values  of  a  and 
various  initial  conditions  is  shown  in  Fig.  6. 
Some  of  these  curves  are  not  completely  dosed; 
this  is  an  indication  that  “jump”  conditions 
may  occur  due  to  effect  of  energy  source.  The 
study  of  these  effects  requires  analysis  using  at 
least  a  two-step  chemicd  reaction  model  and 


Fig.  11.  Phase  portrait  for  first-order  reaction  for  a  - 
const. 


=  0.  (BIO) 

Therefore,  in  this  case,  when  the  chemical  and 
physical  parameters  fulfill  the  relation 
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p  exp 

'  1 
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1  1  +  M0/ 
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the  conditions  for  free  oscillations  exist 
The  term  /(0,  0)  in  Eq.  Bl  can  be  inter¬ 
preted  as  being  analogous  to  a  “damping”  term 
for  the  dynamic  system.  So  that  increasing  or 
decreasing  /(0, 0)  will  correspondingly  affect 
the  oscillations  frequendes  and  amplitudes. 
The  variations  of  fid,  d)  for  obtaining  the 
conditions  for  the  limit  cycle  to  exist  are  now 
exanmed.  Let  /(0,  0)  vary  in  the  range  0  ^ 
/(0,  0)  ^  fid,  dX  where  /  has  a  finite  value. 
The  charaaeristics  of  nonlinear  differential 
equations  indicate  that,  if  fiO,  0)  <  0  and  since 
/  is  larger  than  fid,  0),  then  the  periodical 
solutions  that  e^t  for  fid,  0)  exist  also  at  the 
conditions  of  fid,  0).  Therefore,  osdllations 
(when  th^  occur)  are  notftes&oye^by  the 
increase  of  the  “damping”  term  m  the  equa¬ 
tion  but  are  affected  when  the  stationary  (equi¬ 
librium)  state  becomes  unstable.  The  condition 
0g(0)  >  0  can  be  interpreted  as  being  propor¬ 
tional  to  the  “restoring  force,”  in  an  analogy  to 
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a  mechanical  dynamic  system,  acting  to  return 
the  system  towards  its  equilibrium  state,  to¬ 
wards  the  point  5=0  =  0.  The  term  g(5)  in 
Eq.  B1  represents  the  “restoring  force”  in  the 
dynamic  equation,  therefore,  g(ff)  ^  0.  since 
g(-6)  -giO)  this  restoring  force  is  unhar¬ 
monic. 

Oscillatory  solutions  are  obtained  when 
fid,  9)  changes  sign  as  6  increases  as  the  pro¬ 
cess  develops  with  time.  Also,  starting  from  the 
initial  conditions  5o  =  =  0)  2nd  6^  = 

5(r  =  0),  the  function  +  e,  +  s) 
should  decrease  to  /  <  0.  This  is  illustrated  in 
the  phase  diagram  presented  in  Fig.  12  where 
the  variation  of  f{6,d)  and  5  as  a  function  of 
6  for  various  values  of  a  are  shown.  It  is  seen 
in  this  figure  that  the  curves  for  fid,  6)  first 
decrease  to  negative  values  and  as  d  increases 
this  function  increases  to  positive  values.  Phys¬ 
ically  this  process  can  be  explained  as  indicat¬ 
ing,  first,  as  /  <  0,  the  system  absorbs  the 
external  energy  so  that  the  oscillation  energy  is 
increasing,  then,  as  /  >  0,  the  oscillation  en¬ 
ergy  decreases.  This  pattern  is  seen  for  the  / 
curves  for  a  values  of  0.99  and  1.1.  However, 
for  smaller  values  of  a,  for  example  the  a  — 
0.8  curve,  most  of  the  /  values  are  negative, 
almost  over  the  complete  range  of  tempera¬ 
tures,  indicating  that  the  system  in  this  case 
absorbed  large  amount  of  energy  into  its  oscil¬ 
latory  system,  eventually,  at  sm^er  a  the  os¬ 
cillations  transform  into  fast  combustion- 


Fig-  IZ  fi$,  6)  and  6  variation  with  6  for  first-order 
reaction  for  various  initial  conditions. 


detonation.  On  the  other  hand,  when  a  is 
inaeased  /  is  shifted  to  positive  values,  as 
seen  for  a  =  1.6,  where  the  complete  curve 
has  positive  values  indicating  a  low  tempera¬ 
ture  combustion  without  oscillations. 
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CALCULATIONS  OF  COMBUSTION  IN  A  SCRAMJET  ENGINE  MODEL  USING  THE  EXTERNAL  PROPULSION 

ACCELERATOR  AS  A  TCST  FACE-irY 
J.  Rom^ ,  M.J.  Nusca^ ,  M.J.  Lewis^ ,  A.K.  Gupta^  and  J.  Sabcan^ 


EXTENDED  ABSTTIACT 


ABSTRACT 

The  combustion-detonation  phenomena  in  a  scramjet 
engine  as  well  as  in  an  Oblique  Detonation  Wave  Engine 
can  be  studied  in  the  External  Propulsion  Accelerator. 
This  facility  provides  the  means  to  investigate  the 
combustion-detonation  process  independently  from  the 
fuel  mixing  process.  In  this  case,  a  projectile  simulates  the 
engine  geometry  and  the  combustion  of  the  premixed 
fuel/oxidizer  mixture  is  initiated  in  the  model  flying  in  the 
External  Propulsion  Accelerator.  The  scramjet  combustion 
also  provides  net  thrust  for  the  model  flight  in  the 
Accelerator.  Thus,  this  facility  enables  testing  of  the 
combustion  characteristics  as  function  of  inlet,  combustor 
and  nozzle  geometry,  fuel/oxygen  chemical  parameters 
and  flight  velocity.  The  effects  of  various  fuei/oxidizer 
mixtures  can  be  studied  including  sections  separated  by 
diaphragms  in  the  accelerator  tube.  These  sections  can  be 
filled  with  various  mixtures  to  different  initial  pressures, 
so  that  effects  of  various  pressure  levels  can  also  be 
simulated  in  this  facility. 

INTHODUCnON 

The  External  Propulsion  Accelerator  was  proposed  by  J. 
Rom  (Ref.  1)  and  certain  preliminaiy  aspect  of  its 
characteristics  were  presented  in  Refe.  2,  and  3.  Some 
analytical  investigations,  based  on  many  simplifying 
assumptions,  for  the  establishment  of  combustion  front 
ahead  of  the  forward  facing  step  and  on  spherical  nosed 
blunt  bodies  in  hypersonic  flows  of  detonable  mixtures. 


including  the  studies  of  the  oscillations  that  may  occur 
were  investigated  by  Tivanov  and  Rom  (Ref.  4,  5  and  6). 
A  more  detailed  stucfy  of  the  combustion-detonation 
characteristics  using  the  solutions  of  the  full  Navier- 
Stokes  equations  with  chemical  processes  on  projectiles  in 
the  External  propulsion  Accelerator  is  presented  in  Ref.  7. 
The  External  Propulsion  Accelerator  has  been  developed 
following  the  development  of  the  Ram  Accelerator  by  A. 
Hertzberg  (Ref  8)  at  the  University  of  Washington.  The 
Ram  Accelerator  method  utilizes  the  idea  of  using  a 
premixed  fuel/oxidizer  mixture  in  the  launcher  barrel  and 
injecting  a  projectile  at  supersonic  speeds  into  this 
mixture.  Thrust  is  generated  by  the  ramjet  combustion 
obtained  between  the  tube  wall  and  the  projectile,  where 
the  projectile  acts  as  a  ramjet  centerbody  and  the  barrel  is 
basically  a  long  cowling,  and  combustion  is  initiated  and 
stabilized  by  the  interaction  of  the  shock  waves  from  the 
projectile  nose  with  the  barrel  walls.  However,  as  distinct 
from  the  Ram  Accelerator,  the  External  Propulsion 
Accelerator  use  thrust  that  is  generated  purely  by  the 
interactions  between  the  shock  vmcs  from  the  projectile 
nose  with  combustion  generated  by  normal  shock  waves 
generated  on  the  projectile  bocfy  a  ramp  or  a  forward 
facing  step,  or  a  ring  wing,  completely  independent  of 
the  barrel  walls.  Since  the  projectile  flies  in  a  premixed 
fuel/oxidizer  atmosphere,  the  difficult  problem  of  mixing 
the  fuel  and  oxidizer,  which  plagues  all  airbreathing 
propulsion  methods  proposed  for  hypersonic  flight,  is 
eliininated.  However,  some  important  and  difficult 
combustion  issues,  including  the  ability  to  initiate  and 
stabilize  the  combustion  front  in  the  engine  in  such  a  way 
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that  positive  thrust  can  be  generated,  can  be  investigated 
for  a  uniformly  premixed  mixture  at  various  ^ght 
conditions.  Therefore,  with  proper  simulation  of  the 
intake  geometry  and  the  flow  through  the  engine,  the 
conditions  for  the  initiation  and  the  stabilization  of  the 
combustion  in  the  engine  and  the  corresponding  thrust 
generation  can  be  studied. 

In  this  paper  we  present  the  calculation  of  the  flow  on  a 
scramjet  model  simulating  a  sciamjet  engine  geometiy 
flying  in  the  External  Propulsion  Accelerator  njcing  the 
numerical  programs  developed  for  the  calculations  of  the 
flow  in  the  External  Propidsion  Accelerator.  The  recent 
study,  presented  in  Ref  7,  utilizes  the  solution  of  the 
Navier-Stokes  equations  for  flows  with  chemical  processes 
which  was  adapted  for  the  Ram  Accelerator  calculations 
at  ARL  by  Nusca  (Refs.  9,  10).  Calculations  of  the  flow 
with  combustion  on  the  projectile  designs  proposed  for  the 
External  Propulsion  Accelerator  are  presented  in  Ref  7. 
This  calculation  uses  the  Rockwell  Science  Center  USA- 
RG  (Unifled  Solution  Algorithm  Real  Gas)  code  written 
by  Chakravarthy  et  al  (Refe.  11,  12).  This  CFD  code 
includes  the  full  3D  unsteady  Reynolds-Averaged  Navier- 
Stokes  (RANS)  equations  induding  equations  for 
chemical  kinetics  (finite-rate  and  equilibrium).  These 
equations  are  cast  in  conservation  form  and  converted  to 
algebraic  equations  using  upwind  finifft^Hifn^rgnce  and 
fimte-volume  formulations.  The  equations  are  solved 
using  a  second-order  TVD  (total  variation  diminishing) 
scheme  which  is  used  to  insure  non-oscillatoiy  numerical 
behavior.  Following  these  calculations,  a  scramjet  model 
configuration  will  be  proposed  for  tests  in  the  ARL  Ram 
Accelerator  facility  (R^.  13,  14)  using  subcaliber 
models. 

The  scramjet  model  will  be  composed  of  a  shallow  sharp 
nosed  cone  acting  as  the  scramjet  inlet,  a  ring  wing 
simulating  the  cowling  of  the  scramjet  around  a 
cylindrical  center  body  and  a  conical  afterboefy  simulating 
the  scramjet  exhaust  Combustion  is  generated  behind  the 
shock  wave  established  in  front  of  the  wedge  shape  lip  of 
the  ring  wing.  The  scramjet  combustion  in  this 
configuration  will  produce  high  tenq)erature  high  pressure 
combustion  products  resulting  in  significant  net  positive 
thrust 

HYPERSONIC  SCRAMJET  TEST  FACILITY  BASED 
ON  THE  EXTERNAL  PROPULSION  ACCELERATOR 

In  this  application,  the  hypersonic  scramjet  testing  facility 
will  be  comprised  fi’om:  1.  The  External  Propulsion 
Accelerator  sectiort  2.  Transparent  test  section.  3. 
Deccelerator  section  for  model  retrieval. 

1.  The  External  Propulsion  Accelerator  Section 

The  tests  are  conducted  in  a  regular  External  Propulsion 

Accelerator  using  a  specially  designed  model  for  the 


scramjet  simulation.  The  model  with  its  sabot  are 
accelerated  by  a  gun  powder  charge  to  the  required 
insertion  speed  of  1,000  up  to  1,900  m/sec.  (  Mach 
number  between  3  to  5).  In  a  simple  configuration  the 
projectile  body  will  simulate  the  engine  centerbody  and  a 
properly  shaped  ring  wing  around  this  centeibo^  will 
present  the  engine  cowling.  The  shapes  of  the  engine 
intake  and  the  shape  of  the  combustion  region  inside  the 
engine  can  be  also  simulated  by  this  model.  Combustion- 
detonation  will  be  initiated  by  the  shock  waves  generated 
at  the  intake  to  the  ring  wing  and  it  will  cause  combustion 
in  the  space  between  the  centerbody  and  the  ring-cowling. 

The  flow  will  be  then  e^qpanded  on  the  aft-body  of  the 
model,  simulating  the  engine  exhaust  and  nozzle  flow 
characteristics.  The  operation  of  the  scramjet  wfll  produce 
net  thrust  for  the  acceleration  of  the  model.  Therefore, 
measurements  of  the  model  acceleration-deceleration  will 
enable  the  evaluation  of  the  thrust  obtained  this 
scramjet  engine  as  a  function  of  model  velocity.  More 
sophisticated  measurements  can  be  obtained  using  on¬ 
board  instrumentation  with  telemetering.  Models  with 
hardened  on-board  instrumentation  can  be  built  to  operate 
up  to  50,000  g’s  acceleration  level.  The  acceleration  level 
in  the  projectile  will  depend  on  the  thrust  generated  by  the 
scramjet  Since  the  gas  mixtures  can  be  varied  in  the 
accelerator,  this  will  be  a  method  to  stuefy  the  effects  of 
various  chemical  fuel  compounds  on  the  combustion  and 
the  thrust  generatiorL  Such  studies  can  be  made  when  the 
accelerator  tube  is  divided  into  a  number  of  sections 
separated  by  diaphragms.  Each  section  of  the  accelerator 
tube  can  be  filled  with  a  certain  fuel/oxidizer  mixture  at 
various  initial  pressures.  So  that  in  a  single  experiment 
the  combustion  characteristics  of  flight  in  various 
mixtures  at  specified  initial  pressures  can  be  studied 
These  test  sections  can  be  positioned  at  the  end  of  the 
basic  accelerator  tube  in  which  the  projectile  is  accelerated 
to  the  desired  velocity. 

The  length  of  the  facility  will  be  determined  by  the 
acceleration  of  the  model  and  the  speed  range  for  the  tests. 

The  External  Propulsion  launcher  for  50,000  g’s  will  be 
about  10  meters  for  M  =  7  testing,  40  meters  for  M  =  10 
testing  and  about  150  meters  for  M  =  15  testing.  More 
sophisticated  models  requiring  lower  values  of 
acceleration  can  be  used  by  proper  modifications  to  the 
launcher.  Lower  initial  accelerations  can  be  achieved  by 
initial  acceleration  of  either  a  Ram  Accelerator  stage 
initiated  by  a  small  light  gas  gun  or  a  small  rocket  or  the 
initial  acceleration  be  achieved  directly  by  a  rocket  motor  I 

which  will  separate  at  the  entrance  to  the  External  » 

Propulsion  launcher.  The  length  of  this  launcher  varies  \ 
inversely  as  the  acceleration  levels,  i.e.  lower  acceleration  j 

will  require  the  correspondingly  longer  lengths.  I 

The  External  Propulsion  launcher  diameter  should  be  ? 
about  4  times  the  model  maximum  diameter.  Therefore,  a 
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launcher  of  Im  diameter  will  enable  launching  of  models 
of  25  cm.  diameter  and  the  launcher  tube  diamrtftr  can  be 
reduced  to  40  cm  for  10  cm  models.  Of  course  larger 
models  can  be  considered  with  increase  in  size  of  the 
iacility. 

2.  Transparent  Test  Section. 

The  possibili^  of  including  a  transparent  section  into  the 
Ram  Accelerator  fecility  was  studied  at  ARL  by 
Kruczynski  et  al.  (Ref.  15).  Such  a  transparent  section  at 
the  end  of  the  launcher  tube  will  enable  flow  visualization 
photographs  on  the  model,  particularly,  at  the  intairp  and 
exhaust  regions. 

3.  Deceleration  Section. 

The  energy  imparted  to  the  model  traveling  at  Mach  6  to 
15  is  very  large  and  must  be  absorbed  in  order  to  stop  the 
model.  For  simple  models,  the  models  can  be  expandable 
and  a  conventional  “catcher”  can  be  used.  For  more 
sophisticated  and  e^qransive  models  there  may  be  some 
possibilities  of  some  energy  absorbing  methods  for 
arresting  the  models,  at  least  in  the  lower  hypersonic 
speed  range. 
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Abstract 

The  use  of  the  External  propulsion 
Accelerator  (EPA)  for  launching  models  of 
hypersonic  aerodynamic  configurations  into  an 
instrumented  ballistic  range  is  discussed.  The 
aerodynamic  model  is  encased  inside  an 
axisymmetric  projectile  designed  to  be  accelerated 
to  high  speed  in  the  EPA.  Accelerator  lengths 
required  to  achieve  hypersonic  speeds  are 
estimated  to  vary  from  10  meters  for  Mach  7,  40 
meters  for  Mach  10,  150  meters  for  Mach  15,  and 
700  meters  for  Mach  30,  assuming  a  limit  of 
50,000  g's  acceleration.  For  a  model  span  of  10  cm 
to  25  cm,  the  launch  tube  diameters  are  40  cm  and 
100  cm,  respectively.  Using  this  EPA  launcher  will 
enable  exact  simulation  of  hypersonic  flight  in 
ground  facilities  where  both  the  gas  composition 
and  pressure  can  be  controlled  in  the  ballistic 
range. 

Introduction 

One  of  the  practical  applications  of  the 
External  Propulsion  Accelerator  (EPA)  is  it’s  use 
as  the  launcher  of  aerodynamic  models  to 
hypersonic  speeds  in  a  ground  test  facility.  The 
models,  encased  in  a  properly  designed  sabot,  are 
accelerated  in  the  EPA  to  the  desired  hypersonic 
speed  and  fired  into  an  instrumented  ballistic 
range.  In  the  conventional  ranges  the  model 
velocity  is  limited  to  the  launch  velocities  of  the 
chemical  explosives,  which  is  less  than  2000 
m/sec.,  so  that  these  ballistic  ranges  cannot 
simulate  flight  above  Mach  number  6  in  room 
temperature  air. 
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The  EPA  is  theoretically  capable  to 
launch  projectiles  to  velocities  of  up  to  6  times  the 
detonation  velocity  of  the  fuel-oxidizer  mixture,  in 
the  order  of  12,0()0  m/sec,  so  that  it  is  well  suited 
to  be  a  launcher  for  the  hypersonic  ballistic  range. 
The  requirements  for  a  hypersonic  ground  based 
flight  test  range  and  the  feasibility  of  the  use  of 
light-gas  gun  launcher  was  discussed  by  Witcofski 
et  al  (Ref.  1).  In  this  study  the  use  of  an 
Electromagnetic  launcher  (EML)  and  of  the  Ram 
Accelerator  (RA)  where  also  considered.  It  was 
concluded  that  with  the  present  state  of  technology 
the  most  suitable  candidate  is  the  two-stage  light- 
gas  gun  which  can  be  made  to  launch  models  of 
about  20  cm.  span  and  length  of  up  to  90  cm.  to 
speeds  of  up  to  6  km/sec.  It  was  contemplated,  that 
at  a  later  stage,  using  more  advanced  launcher 
technology,  the  models  can  be  increased  to  30  cm 
-  45  cm  in  span  and  the  velocity  range  will  be 
increased  to  10  km/sec  to  15  km/sec. 

The  free  flight  range  technology  seems  to 
enable  the  best  simulation  of  atmospheric  free 
flight  conditions  in  a  ground  based  facility.  The 
range  concept  is  most  suitable  since  it  uses 
quiescent  arbitrary  gas  in  its  test  sections  in  which 
the  model  flies  at  the  design  hypersonic  speed. 
While  in  other  test  facilities,  such  as  the  various 
types  of  shock  tunnels  and  gun  tunnels,  it  is 
required  to  compress  and  heat  the  test  gas  to  high 
stagnation  temperatures  and  pressures  followed  by 
a  rapid  expansion  in  high  Mach  number  nozzles. 
This  process  causes  disturbances  due  to 
dissociation  and  ionization  of  the  gas  in  the 
stagnation  region  as  well  as  disturbances  due  to 
nozzle  expansion  non-uniformities,  boundary  layer 
effects,  etc.  which  are  encountered  in  such 
facilities. 

The  External  Propulsion  Accelerator  was 
proposed  by  J.  Rom  (Ref.  2)  and  certain 
preliminary  aspect  of  its  characteristics  were 
presented  in  Refs.  3,  and  4.  Some  analytical 
investigations,  based  on  many  simplifying 
assumptions,  for  the  establishment  of  combustion 
front  ahead  of  the  forward  facing  step  and  on 
spherical  nosed  blunt  bodies  in  hypersonic  flows  of 
detonable  mixtures,  including  the  studies  of  the 
oscillations  that  may  occur  were  investigated  by 
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Figure  1.  External  propulsion  accelerator  facility. 


Tivanov  and  Rom  (Ref.  5,  6  and  7).  A  more 
detailed  study  of  the  combustion-detonation 
characteristics  using  the  solutions  of  the  full 
Navier-  Stokes  equations  with  chemical  processes 
on  projectiles  in  the  External  propulsion 
Accelerator  is  presented  in  Ref.  8. 

The  External  Propulsion  Accelerator  has 
been  developed  following  the  development  of  the 
Ram  Accelerator  by  A.  Hertzberg  and  his  co¬ 
workers  (Ref.  9)  at  the  University  of  Washington. 
The  Ram  Accelerator  method  utilizes  the  idea  of 
using  a  premixed  fuel/oxidizer  mixture  in  the 
launcher  barrel  and  injecting  a  projectile  at 
supersonic  speeds  into  this  mixture.  Thrust  is 
generated  by  the  ramjet  combustion  obtained 
between  the  tube  wall  and  the  projectile,  where 
the  projectile  acts  as  a  ramjet  centerbody  and  the 
barrel  is  basically  a  long  cowling,  and  combustion 
is  initiated  and  stabilized  by  the  interaction  of  the 
shock  waves  from  the  projectile  nose  with  the 
barrel  walls. 

In  contrast  to  the  Ram  Accelerator,  in  the 
External  Propulsion  Accelerator  the  thrust  is 
generated  purely  by  the  interactions  between  the 
shock  waves  from  the  projectile  nose  with 
combustion  generated  by  normal  shock  waves 
generated  on  the  projectile  body  by  a  ramp  or  a 
forward  facing  step  or  by  a  ring  wing,  completely 
independent  of  the  barrel  walls.  A  schematic 
diagram  of  this  facility  is  shown  in  Fig.  1 . 

Since  the  projectile  flies  in  a  premixed 
fuel/oxidizer  atmosphere,  the  difficult  problem  of 
mixing  the  fuel  and  oxidizer,  which  plagues  all 
airbreathing  propulsion  methods  proposed  for  the 
remaining  difficult  issues  of  the  ability  to  initiate 
and  stabilize  the  combustion  front  on  the  projectile 
in  such  a  way  that  positive  thrust  can  be  generated. 
The  previous  studies  (Refs.  3,  4  and  8)  indicated 
that  a  relatively  small  forward  facing  step  on  the 
projectile  shoulder  will  be  able  to  induce  a 
detonation  wave  and  induce  large  combustion 
around  the  projectile  aft-body  and  base  resulting  in 
reasonably  high  thrust  levels. 


The  recent  study,  presented  in  Ref.  8, 
utilizes  the  solution  of  the  Navier-Stokes  equations 
for  flows  with  chemical  processes  which  was 
adapted  for  the  Ram  Accelerator  calculations  at 
ARL  by  Nusca  (Refs.  10,  11).  Calculations  of  the 
flow  with  combustion  on  the  projectile  designs 
proposed  for  the  External  Propulsion  Accelerator 
are  presented.  This  calculation  uses  the  Rockwell 
Science  Center  USA-RG  (Unified  Solution 
Algorithm  Real  Gas)  code  written  by  Chakravarthy 
et  al  (Refs.  12,  13).  This  CFD  code  includes  the 
full  3D  unsteady  Reynolds-Averaged  Navier-Stokes 
(RANS)  equations  including  equations  for 
chemical  kinetics  (finite-rate  and  equilibrium). 

These  equations  are  cast  in  conservation 
form  and  converted  to  algebraic  equations  using 
upwind  finite-difference  and  finite-volume 
formulations.  The  equations  are  solved  using  a 
second-order  TVD  (total  variation  diminishing) 
scheme  which  is  used  to  insure  non-oscillatory 
numerical  behavior.  Following  these  calculations, 
the  forward  facing  step  configuration  is  tested  in 
the  ARL  Ram  Accelerator  facility  (Refs.  14,  15) 
using  subcaliber  models. 

The  projectile  external  shape  is  composed 
of  a  shallow  sharp  nosed  cone  followed  by  a 
forward  facing  step,  a  short  length  of  a  cylindrical 
body  and  a  conical  afterbody  with  a  blunt  base. 
Combustion  is  generated  by  the  detonation  wave 
established  in  front  of  the  step,  using  steps  of  about 
1  mm  height.  The  combustion  results  in  a  region  of 
high  pressure,  high  temperature  gas  engulfing  the 
back  parts  of  the  projectile  resulting  in  significant 
net  positive  thrust.  This  region  of  combustion  is 
confined  by  the  hypersonic  flow  structure 
generated  by  the  interactions  of  the  shock  wave 
ahead  of  the  step  which  is  followed  by  the 
combustion  -  detonation  front.  This  is  an  excellent 
demonstration  of  external  propulsion  in  which 
supersonic  combustion  is  confined  by  the  external 
flow  without  any  structural  boundaries.  The  step 
projectile  geometry  and  non-reacting  flow  structure 
is  illustrated  for  a  1  mm  step  projectile  in  Fig.  2. 
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a.)  b.) 


Figure  2.  External  propulsion  hypersonic 
projectile  with  forward-facing  step,  a.) 
projectile  fired  at  U.S..  Army  research 
Laboratory  120  mm  accelerator  b.) 
computational  solution  for  pressure  contours 
over  this  projectile  (non-reacting  case.) 

The  calculations  in  Ref.  8  indicated  that 
the  thrust  levels  which  can  be  achieved  in  the  EPA 
are  on  the  order  F/pA  =  3.  This  thrust  level  can  be 
increased  with  optimized  configurations  of  the 
projectile  and  more  energetic  mixtures.  The 
axisymmetric  projectile  configuration  needed  for 
the  thrust  generation  in  the  launcher  can  be  used 
as  the  cover  for  the  hypersonic  model  which  can 
be  installed  inside  the  projectile.  Upon  exit  of  the 
hypervelocity  projectile  from  the  launcher  tube  the 
projectile  can  be  opened  an  separated  from  the 
model  which  is  then  free  to  fly  at  the  launching 
hypersonic  speed  in  the  instrumented  range. 


Hypersonic  Test  Facility 

In  this  application,  the  hypersonic  testing 
facility  will  be  comprised  from:  1.  The  Accelerator 
section.  2.  Sabot  separation  section.  3.  Ballistic 
Range  for  aerodynamic  models  free  flight  testing. 
4.  Decelerator  section  for  model  retrieval.  These 
parts  are  described  below: 


Figure  3.  External  propulsion  accelerator 
projectile  types,  a.)  forward-facing  step  with 
shock-induced  combustion  b.)  ring  projectile 
with  concentric  cowl,  which  can  operate  as  a 
scramjet. 

1.  The  Accelerator  Section 

The  accelerator  is  comprised  of  two  parts, 
an  initial  accelerator  for  accelerating  the  model 
with  its  sabot  to  the  required  insertion  speed  of 
1,000  to  1,400  m/sec.  (  M  between  3  to  4)  and  the 
External  Propulsion  launcher  tube  which  will 
accelerate  the  model  with  its  sabot  into  the  test 
range.  The  flight  Mach  Number  in  the  ballistic 
range  can  be  then  set  to  the  required  hypersonic 
speed  with  the  possibility  of  reaching  M  =  30. 

The  type  of  initial  accelerator  and  the 
length  of  the  External  Propulsion  Launcher  depend 
on  the  allowable  acceleration  for  the  model. 
Simple  aerodynamic  shapes  used  for  hypersonic 
flight,  such  as  blunt  reentry  bodies,  slender  shuttle 
type  configurations  and  wave  riders  with  hardened 
on  board  instrumentation  can  be  built  to  operate  up 
to  50  Kg's  acceleration  level.  In  this  case,  a  gun 
powder  launcher  can  be  used  for  the  initial 
acceleration  and  the  acceleration  level  of  the 
External  Propulsion  launcher  can  be  adjusted  to 
the  desired  level. 

The  gun  launcher  will  require  lengths  of 
up  to  5  meters  and  the  External  Propulsion 
launcher  for  50  Kg's  will  be  about  10  meters  for  M 
=  7  testing,  40  meters  for  M  =  10  testing  and  about 
150  meters  for  M  =  15  testing.  For  Mach  number 
30  testing  a  launcher  length  of  about  700  meters 
will  be  required.  For  more  sophisticated  models 
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Figure  4.  Proposed  scheme  for  using  an 
external  propulsion  projectile  for  sabot 
launching  of  a  hypersonic  freeflight  model. 

requiring  lower  values  of  acceleration,  the  initial 
accelerator  can  be  either  a  Ram  Accelerator  stage 
initiated  by  a  small  light  gas  gun  or  a  small  rocket 
or  the  initial  acceleration  can  be  achieved  directly 
by  a  rocket  motor  which  will  separate  at  the 
entrance  to  the  External  Propulsion  launcher.  The 
length  of  this  launcher  varies  inversely  with  the 
acceleration  levels,  i.e.  lower  acceleration  will 
require  the  correspondingly  longer  lengths. 

The  External  Propulsion  launcher  diameter 
should  be  about  4  times  the  model  maximum 
diameter  or  span.  Therefore,  a  launcher  of  Im 
diameter  will  enable  launching  of  models  of  25 
cm.  span  or  diameter  and  the  launcher  tube 
diameter  can  be  reduced  to  40  cm  for  10  cm 
models.  Of  course  larger  models  can  be  considered 
with  increase  in  size  of  the  facility. 

2.  Sabot  Separation  Section. 

Since  the  projectile  shapes  required  for 
the  initial  accelerator  and  for  the  External 
propulsion  Accelerator  should  be  of  a  specified 
axisymmetric  shape,  the  aerodynamic  model  must 
be  encased  in  a  sabot,  as  illustrated  in  Fig.  4.  The 
outside  shape  of  this  sabot  should  be  in  accordance 
with  the  external  propulsion  projectile  design  while 
supporting  the  model  to  withstand  the  acceleration 
levels.  Upon  exit  from  the  accelerator  launcher  the 
sabot  should  be  separated  either  mechanically  or 
by  special  mechanism  which  will  be  initiated  at 
the  exit  and  the  sabot  parts  should  be  retained  in 
this  separation  section  while  the  free  model  flies 
into  the  test  range. 


3.  Ballistic  Range  Section. 


The  free  model  can  now  fly  in  the 
instrumented  ballistic  range  where  its 
instantaneous  trajectory  positions  and  attitudes  can 
be  photographed  and  measured  as  well  as  any  on 
board  instrumentation  can  be  monitored  by 
telemetering  system.  If  we  consider  models  of  10 
cm  to  25  cm  in  diameter  or  span,  such  models  are 
sufficiently  large  to  include  most  details  of  the 
aerodynamic  design,  including  deflected  control 
surfaces  and/or  small  control  jets  actuated  by 
small  rockets  or  gas  supply. 

The  ballistic  range  can  either  be  an 
atmospheric  range  or  be  enclosed  in  a  large  tube 
enabling  controlled  atmosphere  for  simulation  of 
pressure-altitude  effects  and  also  atmospheric 
composition  effects  studies. 


4.  Deceleration  Section. 

The  energy  imparted  to  the  model 
traveling  at  Mach  6  to  30  is  very  large  and  must  be 
absorbed  in  order  to  stop  the  model.  For  simple 
models,  the  models  can  be  expandable  and  a 
conventional  "catcher”  can  be  used.  For  more 
sophisticated  and  expansive  models  there  may  be 
some  possibilities  of  some  energy  absorbing 
methods  for  arresting  the  models,  at  least  in  the 
lower  hypersonic  speed  range. 


Conclusions 

The  concept  of  using  the  EPA  for 
launching  large  aerodynamic  models  to  hypersonic 
speeds  has  been  introduced.  It  has  been  shown 
that  with  acceleration  levels  which  can  be 
achieved  in  the  EPA,  about  50,000  g's,  the  length 
of  launcher  needed  for  hypersonic  speeds  would  be 
reasonable,  about  150  meters  for  Mach  15  flight. 

A  ballistic  range  for  testing  aerodynamic 
models  of  span  25  cm  could  be  constructed  using 
an  accelerator  tube  of  diameter  100  cm.  The 
model  would  be  encased  in  an  EPA  projectile- 
sabot  and  separated  before  entering  the  free  flight 
ballistic  range.  This  would  enable  exact  simulation 
of  hypersonic  flight 
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ABSTRACT 

The  initiation  of  detonation  and  the  subsequent  acceleration  of  the  projectile  in  the  in-tube  chemical 
launchers  for  accelerating  projectiles  to  l^rpervelocity  are  investigated  using  an  energy  balance  analysis. 
The  question  of  initiation  of  detonation  the  hypervelocity  projectile  is  evaluated  by  the  use  of  the  blast 

wave  analogy.  In  this  analysis  the  blast  energy  required  for  the  initiation  of  detonation  is  equated  with 
the  work  done  by  the  drag  force  of  the  traveling  projectile.  This  analysis  can  be  applied  to  the  two  types 
of  in-tube  chemical  laimchers,  the  Ram  Accelerator  and  the  External  Propulsion  Accelerator.  It  is  shown 
that  in  the  External  Propulsion  Accelerator,  where  the  detonation  is  initiated  by  the  hypersonic  flow  of  a 
detonable  gas  mixture  flowing  into  a  forward  facing  step,  the  drag  of  this  blunt  step  is  more  then 
sufficient  to  initiate  detonation  at  projectile  velocity  well  below  the  Chapman-Jouguet  detonation 
velocity.  Furthermore,  the  energy  balance  analysis  is  also  applied  to  evaluate  the  maximum  velocity 
which  can  be  achieved  in  the  in-mbe  accelerators  when  the  available  chemical  reaction  energy  is  utilized. 
It  is  shown  that  in  the  Ram  Accelerator,  due  to  the  high  drag  caused  by  the  choking  of  the  flow  between 
the  projectile  and  the  tube  wall  as  the  projectile  Mach  number  is  inaeasing,  the  maximum  projectile 
velocity  is  limited  to  about  1.2  - 1.3  the  Chapman-Jouguet  detonation  velocity.  It  is  observed  that  the 
experimental  measurements  of  the  velocity  to  “unstart”  in  the  Ram  Accelerator  are  at  most  at  1.15  -  1.2 
times  the  Chapman  Jouguet  detonation  velocity,  which  is  in  agreement  with  the  results  of  this  analysis. 

In  the  External  Propulsion  Accelerator,  the  total  drag  coefficient  of  the  projectile  is  much  lower  and 
remains  almost  constant  at  the  hypersonic  Mach  numbers,  while  the  useful  chemical  reaction  energy 
increases  with  increasing  flight  Mach  number.  Therefore,  the  maximum  projectile  velocity  in  the 
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External  Propulsion  Accelerator  can  reach  the  level  of  6  times  the  Ch^man-Jougoet  detonation  velocity, 
which  can  be  well  above  escape  velocity  into  space. 

NOMENCLATURE 
B  (Q/c,)(E/RP.) 

Cp  specific  heat  coefiBcient  at  constant  pressure 

Cv  specific  heat  coefficient  at  constant  volume 
Cd  drag  coefficient 

d  diameter 

D  detonation  speed 

Dam  Damkohler  number 
e  natural  log  base 

E  activation  energy 

H  step  height 

M  Mach  number 

p  pressure 

Pr  Prandtl  number 

Q  heat  of  chemical  reactions 

R  gas  constant 

Rt  radius  of  nose  of  the  body 
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R.  distance  between  nonnal  shock  and  step  face 

Re  R^olds  number 

t  time 

T  Temperature 

u,  V  velocity  components  in  the  x,  y  direction 
V  effective  velocity 

V  velocity  vector 

X,  y,  z  Cartesian  coordinates 

Z  coefficient  of  chemical  reaction  in  the  Arrhenius  rate  equation 
p  mass  fraction  of  species;  ratio  of  combustion  tube  to  projectile  diameters 
y  specific  heat  ratio 

Tic  Carnot  efBciency 

X  detonation  cell  size 

A  He B/y  v Xr 

0  E(T-T*)/RT*^ 

p  density 

T  induction  time 

T,  measure  of  time  of  flow  between  shock  wave  and  step  face 

Xr  measure  of  time  of  chemical  reaction 
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INTRODUCTION 


The  in-tube  chemical  launchers  for  accelerating  projectiles  to  hypervelocity  utilize,  in  the 
superdetonative  mode  of  operation,  the  possibilities  of  generating  continuous  thrust  initiating 

detonation  in  the  premixed  fiiel/oxidizer  mixture  by  a  shock  wave.  The  first  method  proposed  for  an  in¬ 
tube  chemical  launcher  was  the  Ram  Accelerator,  originated  and  developed  by  A.  Hertzberg  and  his 
colleagues  at  the  University  of  Washington  (Ref.  1).  The  concept  of  the  Ram  Accelerator,  operating  in 
the  superdetonative  mode,  is  based  on  utilization  of  the  scramjet  cycle,  where  the  projectile  acts  as  a  fi%e 
centerbody  and  the  tube  as  an  extended  cowling.  The  projectile  is  fired  at  high  initial  speed  into  the  tube 
which  is  filled  with  the  premixed  fiiel/oxidizer  mixture,  where  the  detonation  velocity  in  this  mixture  is 
lower  than  the  projectile  velocity.  The  sharp  nosed  projectile  diameter  is  sUghUy  less  than  the  tube 
diameter  (typically  70%  to  80% )  therefore  the  nose  shock  wave  is  reflected  from  the  tube  wall  into  the 
projectile  centerbody.  Under  proper  conditions  this  reflected  shock  wave  initiates  a  detonation  process  so 
that  when  the  products  of  the  chemical  reactions  are  expanded  on  the  rear  part  of  the  projectile,  thrust  is 
generated.  Another  method  for  operating  the  chemical  in-tube  accelerator  is  based  on  the  utilization  of 
the  external  propulsion  cycle,  proposed  ly  Rom  (Ref.  2,3).  In  deference  from  the  Ram  Accelerator,  in 
the  External  Propulsion  Accelerator  the  projectile  diameter  is  much  smaller  than  the  tube  diameter  so 
that  there  is  no  interaction  between  the  flow  over  the  projectile  and  the  tube  wall  over  the  complete 
length  of  the  projectile.  In  this  case  the  detonation  is  established  by  aerodynamic  means  on  the  projectile, 
such  as  a  forward  facing  step  on  the  projectile  shoulder  (or  a  blunt  leading  edge  of  a  ring  wing  positioned 
on  the  rear  part  of  the  projectile).  In  tlie  case  of  the  projectile  with  the  forward  facing  step,  a  normal 
shock  wave  is  estabUshed  ahead  of  the  forward  facing  step  in  supersonic  and  hypersonic  flows.  When  a 
detonable  gas  mixture  is  flowing  at  these  high  speeds  into  the  forward  facing  step,  this  normal  shock 
wave  compresses  and  heats  up  the  mixture  to  very  high  pressures  and  temperatures.  For  high  enough 
projectile  velocity,  the  heating  and  compression  of  the  combustible  mixture  ty  this  shock  wave  can  result 
in  initiation  of  a  detonation  front  in  the  region  between  the  shock  wave  and  the  face  of  the  step.  Due  to 
the  interaction  between  the  shock  wave  from  the  projectile  nose  with  the  detached  shock-detonation  wave 
established  in  front  of  the  step,  the  high  pressme  high  temperature  chemical  reaction  products  are 
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confined  1^  the  contact  sur&ce  generated  at  the  wave  intersection  point.  Therefore,  the  aeixxlynaniically 
confined  region  filled  with  the  chemical  reaction  products  is  in  effect  an  “external  combustion  chamber” 
which  extends  fiirther  on  the  rear  part  and  in  the  base  region  of  the  projectile  and  is  used  to  produce 
thrust  on  the  projectile. 

It  is  obvious  that  the  projectile  configurations  and  therefore  the  initiation  of  detonation  and  the 
combustion  process  are  very  different  in  the  cases  of  the  Ram  Accelerator  (RA)  and  in  the  External 
Propulsion  Accelerator  (EPA).  In  the  EPA,  using  the  projectile  configuration  with  a  forward  &cing  step, 
the  detonation  is  initiated  by  the  detached  (nearly  normal)  shock  wave  which  is  established  at  hypersonic 
speed  in  fiont  of  the  step.  An  investigation  of  the  required  step  height  and  of  the  position  of  the  reaction 
fi-ont  ahead  of  the  step  was  presented  by  Tivanov  and  Rom  (Ref  4).  Furthermore,  there  are  numerical 
solutions  for  the  flow  over  a  projectile  with  a  step,  presented  Rom  et  al  (Ref  5)  and  an  example  of 
such  a  calculation  is  also  included  in  a  paper  by  Nusca  (Ref  6).  It  was  shown  by  Tivanov  and  Rom  (Ref 
4)  that  step  heights  of  less  than  1  mm  may  enable  initiation  of  detonation  for  hydrocarbon  fiiels,  such  as 
methane  and  acetylene,  as  well  as  hydrogen,  where  all  of  these  fuels  are  premixed  with  oxygen.  Since 
these  studies  are  based  on  rough  approximations,  it  is  important  to  investigate  further,  based  on  the 
present  knowledge  of  detonation  processes,  the  initiation  of  detonation  and  the  type  of  combustion 
process  that  we  may  expect  around  the  forward  facing  step  and  the  over  the  projectile  in  the  EPA 

In  the  present  paper  the  initiation  of  detonation  and  the  subsequent  acceleration  of  the  projectile  in  the 
in-tube  chemical  launchers  for  accelerating  projectiles  to  hypervelocity  are  investigated  using  an  energy 
balance  analysis.  The  question  of  initiation  of  detonation  by  the  hypervelocity  projectile  is  evaluated  by 
the  use  of  the  blast  wave  analogy.  In  this  analysis  the  blast  energy  required  for  the  initiation  of 
detonation  is  equated  with  the  work  done  Ity  the  drag  force  of  the  traveling  projectile.  In  the  case  of  the 
Ram  Accelerator,  due  to  the  large  drag  caused  by  the  choking  of  the  flow  between  the  projectile  and  the 
tube  wall  as  the  projectile  Mach  number  increases,  the  drag  work  is  sufficient  to  induce  detonation  after 
the  projectile  is  accelerated  to  above  a  certain  velocity.  It  is  shown  that  in  the  External  Propulsion 
Accelerator,  where  the  detonation  is  initiated  by  the  hypersonic  flow  of  a  detonable  gas  mixture  flowing 
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into  a  forward  feeing  step,  the  drag  of  this  blunt  step  is  more  then  sufheient  to  initiate  at 

projectile  velocity  well  below  the  Chapman-Jouguet  detonation  velocity. 

The  energy  balance  analysis  is  also  applied  to  evaluate  the  maxinuun  velocity  which  can  be  achieved  in 
the  in-tube  accelerators.  It  is  shown  that  in  the  Ram  Accelerator,  due  to  the  high  drag  caused  by  the 
choking  of  the  flow  between  the  projectile  and  the  tube  wall,  the  maximum  projectile  velocity  is  limitfvi 
to  about  1.2  - 1.3  the  Chapman-Jouguet  detonation  velocity.  It  is  observed  that  the  experimental 
measurements  of  the  velocity  to  “unstart”  in  the  Ram  Accelerator  are  at  most  1.15  - 1.2  times  the 
Chapman  Jouguet  detonation  Mach  number,  which  are  in  agreement  with  the  results  of  this  analysis. 

When  this  analysis  is  applied  to  the  External  Propulsion  Accelerator,  it  is  noted  that  the  total  drag 
coefficient  of  the  projectile  is  much  lower  and  remains  almost  constant  at  the  hypersonic  Mach  numbers. 
In  this  case  the  available  chemical  reaction  energy  increases  with  increasing  flight  Mach  number.  It  is 
shown  that  the  maximum  projectile  velocity  in  the  External  Propulsion  Accelerator  can  be  in  the  order  of 
6  times  the  Chapman-Jouguet  detonation  velocity,  which  is  well  above  escape  velocity  into  space. 


INITIATION  OF  DETONATION  BY  A  BLUNT  BODY 

A  discussion  on  the  initiation  of  detonation  on  a  blunt  nosed  projectile  is  presented  by  Lee  (Ref.  7,8).  It  is 
proposed  there  to  estimate  the  requirement  for  initiation  of  detonation  by  a  hypervelocity  projectile  by  the 
use  of  the  blast  wave  analogy,  thus  equating  the  blast  energy  required  for  the  initiation  of  detonation  with 
the  work  done  by  the  drag  force  of  the  traveling  projectile.  Various  results  of  the  measurement  of  the 
energy  for  direct  initiation  of  detonation  are  discussed  by  Lee  (Refs.  7,8).  The  analysis  is  then  based  on 
the  relation,  first  stated  by  Zeldovich,  that  this  critical  energy  is  dependent  on  the  third  power  of  the 
induction  time  (or  induction  length)  for  spherical  blast  waves.  Using  strong-blast  theory,  Lee  derived  the 
following  expression  for  the  critical  energy,  E„ ,  for  spherical  detonation 

=  \.692np^D^{25DT^f  (1) 
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and  for  cylindrical  detonation 


E,  =  \2527tp,D\2DT,f  (2) 

where  D  is  the  detonation  velocity  for  the  Chapman-Jouguet  conditions  and  %,  is  the  induction  time.  In 
Ref.  8  Lee  states  that  “if  the  induction  time  evaluated  at  the  shock  tenq>eratuie  corresponding  to  a 
Chapman-Jouguet  detonation  is  used  for  x,  than  the  value  of  E,  evaluated  from  Eq.  1  is  about  three  orders 
of  magnitude  smaller  than  the  experimental  values”.  The  e^iimental  values,  to  which  Lee  refers,  are 
determined  from  the  measurements  of  the  average  cell  size,  X,  and  from  the  evaluation  of  the  critical  tube 
diameter  determined  for  detonation  waves  traveling  into  a  stationary  detonable  mixture.  Lee  then 
assumed  that  the  characteristic  chemical  length  is  the  cell  size,  so  that  the  characteristic  time  is 
proportional  to  X/D. 

Using  semi-empirical  considerations,  Lee  suggests  the  following  relations  for  the  energy  for  direct 
initiation,  for  the  spherical  and  cylindrical  detonation  cases  for  fiiel-air  mixtures. 

For  spherical  detonation  (Ref  7), 

E,=252p,M^^X^  (3) 

For  cylindrical  detonation,  the  critical  energy  per  xmit  length  is  (Ref  7), 

E^^\4.5p^MlX^  (4) 

where  X,  the  cell  size,  is  related  to  the  induction  time.  The  cell  size  was  measured  for  detonations 
traveling  into  a  stationary  gas  mixture  at  relatively  low  pressures,  po  in  the  range  of  10  torr.  to  200  torr. 
and  up  to  atmospheric  pressure. 

The  large  discrepancy  between  the  critical  energy  which  is  evaluated  by  applying  the  strong-blast  theory 
as  presented  in  Eqs,  1  and  2,  using  the  induction  time,  % ,  evaluated  at  the  temperature  corresponding  to 
the  Chapman-Jouguet  detonation  and  the  critical  energy  evaluated  using  the  experimentally  measured 
cell  size,  Eqs.  3  and  4,  indicates  a  strong  deficiency  in  our  ability  to  evaluate  correctly  this  energy  for 
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initiation  of  detonation  on  a  traveling  projectile.  It  seems  reasonable  to  expect  that  in  the  case  of  the  flow 
of  a  detonable  mixture  at  hypersonic  speeds  in  and  about  the  stagnation  region  of  a  blunt  body  the 
characteristic  time  for  detonation  should  be  the  induction  time  evaluated  from  the  chemical  reaction  rates 
at  the  temperature  and  pressure  corresponding  to  the  Chapman-Jouguet  detonation.  This  discrepancy, 
which  was  pointed  out  by  Lee  in  Ref.  7,  should  be  investigated.  Particularly  since  there  are  number  of 
numerical  calculations  Yungster  et  al  (Ret  9),  Nusca  (Ret  6 )  and  Li  et  al  (Ret  10)  and  supported  Ity 
some  RA  experiments  which  indicate  that  the  initiation  of  detonation  is  obtained  and  is  computed  using 
the  induction  time  evaluated  from  the  chemical  reaction  rates.  In  this  respect,  it  will  be  extremely 
important  to  measure  directly  the  initiation  of  detonation  on  a  flying  projectile  in  the  EPA  where  the 
detonation  will  be  initiated  by  a  forward  facing  step. 

When  we  consider  the  initiation  of  detonation  ahead  of  a  blunt  body  in  hypersonic  flow,  it  was  also 
pointed  out  by  Zeldovich  that  detonation  will  occur  when  the  induction  time  is  shorter  than  the  flow  time 
from  the  detached  shock  wave  to  the  free  of  the  body.  The  induction  time  is  determined  by  the  chemical 
reaction  kinetics  of  the  mixture  at  the  conditions  behind  the  detached  shock  wave  while  the  flow  time  is  a 
function  of  the  detachment  distance  which  varies  with  the  radius  of  curvature  of  the  blunt  nose  of  the 
botfy.  Analysis  of  certain  aspects  of  the  hypersonic  flow  over  blunt  bodies  with  chemical  reactions  were 
presented  by  Sichel  and  Galloway  (Ref.  1 1) ,  Galloway  and  Sichel  (Ret  12)  and  Chemyi  (Ref  13).  It  was 
shown  in  these  studies  that  the  establishment  of  combustion  in  the  stagnation  region  depends  on  the  ratio 
of  the  ignition  delay  distance  to  the  shock  standoff  distance.  Sichel  and  Galloway  (Ref.  1 1)  define  a 
parameter  relating  the  ignition  delay  time,  the  body  radius  of  curvature  and  the  flow  velocity.  When  this 
ignition  delay  parameter,  k  =  u«  t  /  IL  (where  u«  is  the  free  stream  velocity,  x  is  the  ignition  delay  time 
and  Rt  is  the  spherical  nose  radius),  is  less  than  1,  then  combustion  will  occur  between  the  detached 
shock  wave  and  the  body  surface. 

Studies  of  hypersonic  flows  with  oblique  detonation  are  presented  in  Refs.  14  and  15.  A  study  of  the 
hypersonic  flow  of  a  premixed  fiiel/oxidizer  mixture  over  a  wedge  and  the  characteristics  of  the  oblique 
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detonation  wave  is  presented  Powers  and  Stewart  (Ref.  14)  and  studies  of  the  oblique  detonation  wave 
by  Li  et  al  in  Ref  15. 


INITIATION  OF  DETONATION  BY  A  FORWARD  FACING  STEP 

Since  the  forward  facing  step  has  an  infinite  radius  of  curvature  at  the  stagnation  region,  the  initiation  of 
detonation  will  occur  when  the  temperature  behind  the  detached  shock  wave  will  be  above  the  minimum 
auto-ignition  temperature  of  the  mixture.  Therefore,  above  a  minimum  projectile  velocity  the  temperature 
behind  the  detached  shock  wave  will  be  higher  then  the  auto-ignition  temperature,  and  a  detonation 
process  will  be  initiated.  Above  the  step,  the  shock  wave  will  curve,  reducing  the  radius  of  curvature  as 
the  distance  above  the  step  increases.  As  discussed  by  Lee  in  Ref  7,  when  the  shock  wave  curvature  is 
reduced  below  a  critical  value  and  hence  the  flow  divergence  is  too  severe,  then  the  reactions  may  be 
quenched  even  when  the  initial  shock  temperature,  in  the  stagnation  point,  exceeds  the  auto-ignition 
limit  Therefore,  in  order  to  insure  stabilization  of  the  reaction  firont  behind  the  detached  shock  wave  a 
minimum  step  height  is  needed.  This  qualitative  argument  is  in  agreement  with  the  results  of  the 
analysis  of  Tivanov  and  Rom  (Ref  4)  where  the  minimum  step  height  concept  was  introduced  and 
evaluated  by  an  approximate  analysis.  In  the  analysis  in  Ref  4,  it  was  shown  that  when  B  t,/  Tr>  1  /  e 
then  detonation  can  occur.  Here,  the  time  of  reaction,  Xr  /  B,  and  the  time  of  the  flow  between  the 
detached  shock  wave  and  the  body,  r, ,  are  estimated  by  the  following  relations: 

T,‘  =  (Z/  p)  exp(-E/RT.)  t,  =  H/  v 

where  H  is  the  step  height  and  v  is  an  average  velocity  in  the  stagnation  region  which  can  be  evaluated 
from  the  solution  of  the  hypersonic  blimt  body  stagnation  flow  field  (Ref  16).  A  minimum  step  height 
needed  for  estabUshing  a  detonation  fi'ont  is  then, 

Hix,  v/eB 
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This  required  step  height  can  be  defined  in  terms  of  a  Damkohler  number.  Dam  =  H  /  Xr  v ,  which  must 
be  equal  or  greater  then  1/  eB  for  detonation  to  be  initiated  and  stabilized,  and  is  determined  by  the 
mixture  properties.  This  minimum  step  height  was  evaluated  in  Ref  4  and  is  presented  for  mixtures  of 
hydrogen-o?^gen,  acetylene-oxygen  and  methane  oxygen  as  a  function  of  the  flow  Mach  number  in  Fig. 
la  and  lb.  It  is  possible  to  define  a  parameter  A  =  HeB/Tr  v,  which  indicates  the  possibility  of 
detonation.  If  A  >  1  then  detonation  can  develop  while  for  A  <  1  only  deflagration  can  occur.  This 
analysis  is  in  agreement  with  the  phenomenological  analysis  of  Sichel  and  Galloway  (Ref  5).  Their 
criterion  for  initiation  of  the  combustion  fi'ont  on  a  blunt  body  represented  by  the  ignition  delay 
parameter  k  =  u«  x  /  Rb  (where  x  is  the  ignition  delay  time  and  Rb  is  the  spherical  nose  radius)  is 
inversely  proportional  to  the  parameter  A. 

THE  ENERGY  REQUIRED  TO  INITIATE  DETONATION 

In  order  to  define  the  conditions  for  the  establishment  of  detonation  on  a  projectile,  consider  the  flow  of 
the  fiiel/oxidizer  mixture  at  hypersonic  velocity  over  the  conical  nose  of  the  projectile,  as  illustrated  in 
Fig.  2.  An  oblique  shock  wave  is  formed  on  the  cone,  attached  to  the  cone  tip,  compressing  and  heating 
the  gas  mixture.  The  nose  cone  angle  is  selected  so  a  to  keep  the  temperature  behind  the  conical  shock 
wave  low  enough  so  as  to  prevent  significant  chemical  reaction  from  occurring  on  the  projectile  nose  due 
to  this  shock  wave  over  the  range  of  the  design  flight  Mach  numbers.  A  forward  facing  step  is  positioned 
at  the  shoulder  of  the  cone  causing  a  normal  shock  wave  ahead  of  the  step  which  compresses  and  heats 
up  the  mixture  so  that  exothermic  reactions  are  sustained  between  the  shock  wave  and  the  face  of  the 
step.  As  the  temperature  and  pressure  behind  the  detached  normal  shock  wave  increase  the  reaction  front 
moves  closer  to  the  shock  wave  until  a  detonation  front  is  established  in  front  of  the  step.  As  the  flows 
nears  the  outer  comer  of  the  step,  the  expansion  fan  from  the  comer  will  weaken  and  curve  the  initially 
normal  detached  shock  wave.  The  reduced  radius  of  curvature  above  the  step  comer  will  reach  a  value 
where  the  reaction  front  will  separate  from  the  curved  shock  wave.  However  for  flows  of  high  enough 
Mach  numbers  the  reaction  will  be  attached  to  the  detached  shock  wave  at  least  till  it’s  intersection  point 
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with  the  conical  shock  wave  fiom  the  nose.  So  that  all  the  incoming  flow  into  the  projectile  nose  iq)  to 
the  intersection  point  between  the  detached  shock  wave  and  the  conical  shodc  wave  will  pass  through  the 
reaction  front  initiated  by  the  detached  shodc  wave.  Furthermore,  at  increasing  flow  Mach  numbers  the 
transmitted  shock  wave  b^rond  the  intersection  may  be  still  strong  enough  to  sustain  the  reaction  front 
further  into  the  flow  above  the  projectile  until  the  local  radius  of  curvature  of  the  shodc  wave  will  be 
reduced  so  as  to  quench  the  reaction.  In  this  case  the  layer  of  the  high  temperature  high  pressure  reaction 
products  confined  the  contact  sheet  acts  as  an  “external  combustion  chamber”  for  the  projectile  and  its 

etiergy  can  be  used  to  generate  thrust  on  the  rear  part  of  the  projectile. 

A.  Evaluation  of  the  Minimum  Mach  Number  for  Initiation  of  Detonation  on  the  Hypervelodty  Projectile 

The  minimum  Mach  number  for  the  initiation  of  detonation  on  a  projectile  with  a  forward  facing  step  in 
the  EPA  can  be  evaluated  by  equating  the  critical  energy  required  for  direct  initiation  of  cylindrical 
detonation,  given  Eq.  2  with  the  work  done  by  the  drag  force  of  the  hyperveloci^  projectile  of 
diameter  d  with  a  step  H  (where  H/d  «  1),  thus 

E,  =  \2527rp,D\2Dr,f  =  05  p„V^{2ndH)Cj,^  (5) 


Thus  the  minimum  projectile  speed  in  this  case  is  expressed,  as  follows: 


HfDzy  1 
-^  =  2.504—1 


D 


d\  H  J 


(6) 


Since  typical  H/d  values  can  be  1/30  and  smaller ,  Dr,  /H  should  be  less  than  l(probably  in  the  range  of 
1/3  to  1/5)  for  detonation  to  occur  and  Cd*,,  for  a  forward  facing  step  is  about  1.  Therefore,  introducing 
these  values  into  Eq.  6,  the  minimmn  projectile  velocity  that  produces  energy  to  balance  the  critical 
initiation  energy  for  detonation  is  evaluated  to  be  much  below  the  detonation  velocity.  Since  the 
projectile  in  the  EPA  is  expected  to  fly  well  above  the  detonation  speed  then  there  is  no  problem  of 
securing  the  energy  required  to  initiate  the  detonation. 
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The  energy  balance,  presented  in  Eq.  5,  can  be  ^plied  also  to  the  projectile  in  the  RA.  In  this  case,  the 
drag  coefiBcient  is  equated  with  that  of  the  drag  of  the  projectile  including  the  effects  of  the  int^rarHon 
between  the  projectile  and  the  tube  wall  and  the  reference  area  is  the  cross  section  of  the  projectile  body. 
As  the  projectile  speed  increases  the  drag  force  on  the  projectile  increases  so  as  to  provide  the  energy  for 
initiation  of  the  detonation.  At  the  lower  initial  speeds  the  required  drag  is  provided  by  the  perforated 
piston  oibutrator  which  is  used  in  the  RA  for  initiation  of  the  process. 

B.  Evaluation  of  the  Maximum  Velocity  by  the  Energy  Balance  on  the  Ifypervelocily  Projectile  in  the 
EPA 


As  discussed  before,  the  work  done  the  drag  force  acting  on  the  hypervelocity  projectile  is  balanced  by 
the  energy  released  by  the  detonating  gas  mixture.  Following  the  discussion  by  Lee  in  Ref  7,  it  is 
possible  to  evaluate  the  maximum  velocity  that  can  be  reached  by  the  projectile  when  all  the  combustion 
energy  of  the  reacting  mixture  is  utilized.  This  energy  balance  can  be  e}q)ressed  by  the  following  relation. 


(7) 


where  Q  is  the  chemical  energy  per  unit  mass  and  qc  is  the  Carnot  cycle  efficiency.  Now, 

and  then  using  Eq.  7  it  can  be  shown  that 


cAr^-i) 


(9) 


and  following  Lee  (Ref  7),  we  can  assume  a  value  of  2/3  for  t|c  (assuming  high  combustion  temperature 
of  3000**  K  and  low  temperature  of  1000°  K).  In  Eqs.  7  and  9,  (3  is  defined  as  the  ratio  of  the  outside 
diameter  of  the  external  combustion  zone  to  the  projectile  diameter.  In  the  EPA  the  thickness  of  the 
external  combustion  layer  varies  as  a  function  of  the  flow  Mach  number.  It  was  found  firom  the  numerical 
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calculations,  which  are  rqx)rted  in  part  in  Re&.  5,6  that  for  the  32  nun  diameter  projectile  with  a  1  mm 
step  the  combustion  layer  thickness  varies  from  3  step  heights  at  Mach  number  5  to  5.S  step  heights  at 
Mach  number  6  and  to  17  step  heights  at  Mach  number  10.  The  value  of  Cd  for  the  cone  with  the  1  mm 
step  is  about  0.26-0.28  for  Mach  numbers  above  5  (for  cone-cylinder  without  the  step  the  drag  coefficient 
is  about  0.07-0.08).  Then,  using  Eq.  9,  the  maximum  projectile  velocity  can  be  about  1.9  times  the 
detonation  velocity  at  Mach  number  5,  about  2.1  times  the  detonation  velocity  at  Mach  number  6  and 
about  3.3  times  the  detonation  velocity  at  Mach  number  10.  At  higher  flight  Mach  numbers  the 
combustion  may  extend  to  the  tube  wall,  then  for  the  120  mm  tube  with  the  32  mm  projectile  the 
maximum  velocity  is  about  6  times  the  detonation  velocity,  which  is  well  above  the  escape  velocity  to 
orbit  into  space.  Therefore,  when  the  projectile  is  injected  at  initial  velocity  between  Mach  5  to  6,  only 
about  21%  to  28%  of  the  energy  of  combustion  is  used  to  overcome  the  drag,  so  that  about  72%  to  79% 
of  the  combustion  energy  is  available  for  accelerating  the  projectile.  Actually  the  available  combustion 
energy  increases  as  the  Mach  number  increases,  reaching  to  about  90%  at  Mach  number  10  and  will  be 
even  higher  as  the  Mach  number  is  increased  further.  In  this  case  of  accelerating  the  projectile  in  the 
EPA  we  can  conclude  that  only  a  small  fraction  of  the  available  combustion  energy  which  can  be 
released  in  the  detonation  process  is  needed  to  overcome  the  drag  of  the  hypervelocity  projectile  and  this 
available  energy  can  be  used  for  accelerating  the  projectile  to  higher  hypervelocity.  Of  course  as  the 
projectile  velocity  increases  the  problems  of  heating  and  ablation  of  the  projectile  surfaces  and  the 
aerodynamic  and  acceleration  loads  on  the  projectile  structure  become  more  critical. 

C.  Evaluation  of  the  Maximum  Velocity  by  the  Energy  Balance  on  the  Hypervelocity  Projectile  in  the  RA 

It  is  of  interest  to  evaluate  the  maximum  velocity  of  the  projectile  in  the  RA  by  this  ener^  balance 
analysis.  Particularly,  since  there  are  some  experimental  results  on  the  operational  limits  in  the  RA.  Thus 
using  a  similar  analysis,  we  can  apply  the  energy  balance  and  evaluate  the  maximum  velocity  applying 
Eq.  9  to  the  case  of  the  projectile  in  the  RA.  Such  evaluation  was  first  done  by  Lee  in  Ref.  7.  There  Lee 
assumes  a  value  of  the  drag  coefficient  for  the  projectile  flying  in  the  RA  tube  to  be  of  order  1,  Cd  =  0(1), 
which  is  equal  to  the  drag  coefficient  of  a  blunt  body  at  hypersonic  speeds.  So,  using  Eq.  9  with  y  =  1.4 
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and  P  =  1.5,  he  estimates  the  value  of  the  maximum  projectile  velocity  to  be  about  1.3  times  the 
detonation  velocity.  He  then  concludes  that  when  the  projectile  velocity  is  equal  to  the  detonation  velocity 
about  75%  of  the  available  energy  is  used  to  overcome  the  drag  so  that  only  a  small  fraction  of  the  energy 
is  available  for  additional  acceleration  of  the  projectile. 

The  projectiles  used  in  the  RA  of  the  University  of  Washington  and  ARL,  as  presented  in  Refs.  1  and  6, 
has  a  drag  coefficient  of  about  0.2  >0.25  at  Mach  numbers  above  5  in  fee  ffigbt  without  the  tube  wall 
interference.  However,  the  multiple  shock  wave  reflections  from  the  tube  wall  to  the  projectile  center- 
section  and  the  blunt  leading  edges  of  the  fins  which  may  initiate  detonation  increase  the  drag  coefficient 
of  the  projectile  considerably  in  its  flight  in  the  accelerator  tube.  As  the  projectile  velocity  increases  the 
strength  of  the  reflected  oblique  shock  waves  increase  and  the  total  pressure  loss  in  the  shock 
interference  region  on  one  hand  provides  energy  for  the  initiation  of  the  detonation  but  on  the  other  hand 
due  to  the  extremefy  high  drag  consumes  an  increasing  portion  of  the  available  reaction  energy  until  the 
energy  required  to  overcome  the  drag  force  wiU  be  larger  than  the  available  energy  from  the  chemical 
reaction  and  we  will  fece  the  conditions  of  “unstart”.  At  these  conditions  the  assumption  of  a  drag 
coefficient  of  about  1,  which  is  the  value  for  a  blunt  Ixxty,  seems  reasonable. 

An  experimental  investigation  as  to  the  limit  of  operation  in  the  RA  was  conducted  at  the  University  of 
Washington  (Ref  17).  The  RA  tube  diameter  is  38  mm  while  the  diameter  of  the  projectile  body  is  29 
mm.  Therefore,  using  Cd  =  0.8-0.9  in  Eq.  9,  and  the  appropriate  values  of  the  other  parameters,  i.e.  y  = 
1.4  and  P  =  1.31,  we  obtain  that  the  maximum  value  of  the  projectile  speed  in  the  University  of 
Washington  RA  facility  is  about  1.2-1. 3  times  the  detonation  velocity.  The  experiments  were  conducted 
with  groups  of  mixtures.  The  first  group  was  based  on  a  methane-oxygen  diluted  with  varying  amounts 
of  nitrogen,  that  is  2.8  Oh  +  2  O2  +  X  Nj .  The  second  group  was  based  on  a  hydrogen-oxygen  diluted 
with  methane,  that  is  2  H:  +  2  Oj  +  X  CHi .  The  five  firmed  projectile  was  fired  into  these  mixtures  and 
its  flight  characteristics  measured.  It  was  foimd  that  for  the  very  energetic  mixtures,  where  the  dilutant 
was  very  small,  the  projectile  unstarted  at  relatively  low  velocity,  even  below  the  Chapman-Jouguet 
velocity.  This  can  be  explained  by  the  fact  that  for  these  very  energetic  mixtures  the  energy  release  by  the 
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detonation  generated  the  reflected  shock  wave  is  so  large  that  it  chokes  the  flow,  inn-Kaging  the 
effective  drag  of  the  projectile  so  that  the  projectile  is  “stopped”  -  then  the  “unstart”  in  the  terminology  of 
Ref.  17  occurs  since  the  stopped  projectile  cants  hits  the  wall  and  is  destroyed.  However,  for  moderate 
dilutants  the  projectile  accelerated  from  an  initial  Mach  number  of  3.8  to  above  the  Chapman-Jouguet 
detonation  velocity  reaching  Mach  numbers  up  to  about  5.8.  In  these  experiments,  as  the  projectile 
velocity  increased,  at  a  certain  speed  the  projectile  experienced  the  phenomenon  of  “unstart”,  which  was 
generally  followed  by  disintegration  of  the  projectile.  These  measurements  of  the  projectile  velocity  to 
“unstart”,  which  are  reported  in  Ref.  17,  are  presented  in  Figs.  3a  and  3b,  for  the  nitrogen  diluted 
mixtures  and  for  the  methane  diluted  mixtures,  respectively.  Observing  the  results  of  the  mftagirpmpnt  g 
of  the  Mach  number  at  “unstart”  conditions,  it  turns  out  that  the  ratio  of  M^  Mo  is  constant  for  all  the 
runs  with  the  methane-oxygen  with  the  nitrogen  dilutants,  M— ./  Ma  =  1.2,  and  for  the  hydrogen-oxygen 
with  methane  dilutants  the  value  was  also  constant,  Ma  =  1. 15  (as  shown  in  the  tables  attached  to 
Figs.  3a  and  3b).  This  is  in  remarkable  close  agreement  with  the  energy  balance  analysis  results  (when 
we  identify  that  the  maximum  velocity  is  just  above  the  velocity  to  unstart),  which  was  calculated  to  be, 

Motw/  Ma  =  1.2-1.3.  We  find  that  there  is  an  excellent  qualitative  agreement,  the  fact  that  M. _ J  Ma  = 

constant,  and  also  reasonable  quantitative  agreement,  less  than  10%  deviation  in  the  value  of  the 
maximum  Mach  number,  considering  the  rough  evaluation  of  the  drag  coefScient  for  the  projectile  with 
the  tube  walls  interference. 

In  the  discussion  in  Ref.  17,  it  is  suggested  that  the  “unstart”  phenomenon  is  caused  by  projectile 
disintegration  due  to  configuration  distortion  by  aerodynamic  heating  and  mechanical  stresses.  It  is 
assumed  that  this  distortion  is  a  result  of  the  aerodynamic  heating  and  pressure  and  acceleration  loads 
which  cause  the  projectile  to  cant  and  be  destroyed,  as  shown  in  Ref.  18.  A  test  of  this  assumption  was 
performed  by  firing  a  titanium  nosed  projectile  instead  of  the  aluminum  projectiles  used  in  all  previous 
tests.  The  titanium  projectile  was  able  to  withstand  the  choking  effects  at  the  maximum  velocity  without 
significant  distorfions  and  actually  exited  from  the  accelerator  tube.  A  flight  trajectory  comparison 
between  the  aliuninum  and  titanium  projectiles  was  presented  in  Ref  17  and  is  shown  in  Fig.  4.  It  is 
interesting  to  note,  that  this  titanium  projectile  reached  a  maximum  velocity  after  traveling  about  13  m  in 
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the  tube  and  for  the  last  2-3  meters  it  continued  at  constant  velocity  of 2300  m/sec.  The  aluminum 
projectile  had  an  unstart  at  a  velocity  of 2200  m/sec  after  traveling  about  10  m  in  the  tube.  For  this 
aluminum  projectile  the  Mach  number  at  unstart  was  1.13  times  the  Chapman-Jouguet  Mach  number 
while  the  mavimiiin  Mach  number  of  the  titanium  projectile  is  1.2  times  Mo ,  which  is  in  close 
agreement  with  the  estimated  value  from  the  energy  balance  analysis. 


SUMMARY  AND  CONCLUSIONS 

It  is  shown  that  the  energy  balance  analysis  of  the  projectile  flight  in  the  in-mbe  chemical  accelerators 
can  be  used  to  evaluate  the  velocity  required  to  initiate  the  detonation  process  on  the  rear  part  of  the 
projectile.  Applying  this  analysis  to  the  projectile  designed  for  the  External  Propulsion  Accelerator, 
where  the  detonation  is  initiated  on  the  blunt  forward  facing  stq),  it  is  shown  that  even  at  velocities 
below  the  Chapman-Jouguet  detonation  velocity  the  drag  work  of  the  step  on  the  projectile  is  sufficient  to 
initiate  the  detonation. 

It  is  shown  that  using  the  energy  balance  analysis  it  is  possible  to  evaluate  the  maximum  velocity  that 
can  be  achieved  in  the  in-tiibe  chemical  accelerators  by  balancing  the  drag  work  of  the  projectile  with  the 
utilization  of  the  total  available  chemical  reaction  energy. 

In  the  Ram  Accelerator  case,  the  high  drag  experienced  by  the  projectile  due  to  the  choking  of  the  flow  at 
increasing  Mach  numbers  results  in  limiting  the  projectile  velocity  to  about  1.2  - 1.3  times  the  Chapman- 
Jouguet  detonation  velocity.  This  result  is  supported  by  the  measurements  of  the  “unstart”  phenomenon 
in  the  RA  at  the  University  of  Washington  (Ref.  17). 

Using  this  analysis,  it  is  shown  that  the  maximum  projectile  velocity  in  the  External  Propulsion 
Accelerator  can  be  few  times  the  Chapman-Jouguet  detonation  velocity  and  may  reach  6  times  this 
velocity,  which  is  well  above  the  escape  velocity  into  space. 
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Fig.  la.  The  calculated  nrinimum  step  height  for  hydrogen  -  oxygen  mixture 
(from  Ref.  4). 


lb.  The  calculate  minimum  step  height  for  the  acetylene  -  oxygen  and  for  the 
methane  -  oxygen  mixtures  (from  Ref  4). 


Heat  Release  Q 


2.  The  flow  Structure  on  the  conical  nose  with  the  forward  facing 
step. 


MIXTURE:  2.8  CH-  +  2  0,  -i*  X  N, 

CJ  Ddonalion 


Mach  Number 


3  a.  Results  of  variable  nitrogen  experiments  plotted  as  thermally 
choked  heat  release  vs.  flight  Mach  number  (from  Ref  17). 


Velocily  (rrVsec) 


CJ  Detonation 


MIXTURE:  2  H,  ♦  2  O,  +  X  Ca, 


3b.  Results  of  variable  methane  experiments  plotted  as  thermally 
choked  heat  release  vs.  fli^t  Mach  number  (from  Ref  17). 


Position  (m) 


4.  Velocity-distance  data  comparing  aluminum  and  titanium  nose 
projectiles  (from  Ref.  17). 


